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Schistura megalodon species nova, а new river loach from the Irra- 


waddy basin in Dehong, Yunnan, China (Teleostei: Cypriniformes: 


Nemacheilidae) 


Marco Endruweit" 


601, Qingshan Road, Licang, Qingdao, Shandong, Qingdao 266001, China 


Abstract: A new species of river loach, Schistura megalodon sp. nov., is described from the Irrawaddy basin in Yingjiang County, 


Dehong Autonomous Prefecture, Yunnan Province, China. The following combination of diagnostic characters serve to distinguish it 


from all other congeners in the given zoogeographical region: a large processus dentiformes in the upper jaw, a short pre-anus length 
of 65.4%-66.3% of SL, long paired fins (pectoral: 20.8%-24.2% of SL; pelvic: 17.0%-20.6% of SL), a wide body оҒ9.7%-11.3% 
of SL at anal fin origin, an incomplete lateral line, the absence of an orbital lobe, and a broad and distinct basicaudal bar with forward 


extensions. 
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The Dehong Dai and Jingpo Autonomous Prefecture 
Is situated at the westernmost corner of China's Yunnan 
Province and shares a border with Myanmar's Kachin 
State. This area 1s quasi-solely drained by the Irrawaddy 
River, also known as Ayeyarwady River, except for its 
southeastern corner which lies in the watershed of the 
Salween River. Yingjiang County as the largest of five 


counties stretches over the northern half of the prefecture. 


It has various climatic zones from temperate to tropical 
with altitudes ranging from 210 to 3 404 m (Wikipedia, 
2014). The region along its border with Myanmar 1s 
inhabited by people of the ethnic minorities Jingpo and 
Lisu. 

Examination of a jar of fishes in the collection room 
of the Kunming Institute of Zoology (KIZ), Chinese 
Academy of Science, from rivers around Tongbiguan 
town in Yingjiang County prompted an ichthyofaunistic 
survey to western Dehong. The survey yielded inter alia 
a lot of four specimens of a Schistura species that was 
compared with all congeners reported from the 
Irrawaddy. Nemacheilid loaches reported from the 
northeastern Irrawaddy are Schistura sikmaiensis (Hora, 
1921), S. vinciguerrae (Hora, 1935), S. polytaenia (Zhu, 
1982), S. yingjiangensis (Zhu, 1982), and S. malaisei 


Science Press 


Kottelat, 
Rendahl, (1948) is considered a junior synonym of S. 
Sikmaiensis (see Kottelat, 1990:208; 2012:118). In 
addition, the concerned lot was compared with affine 


1990. The taxon Nemacheilus putaoensis 


species from neighboring watersheds. The concerned 
specimens turned out to be a distinct species which is, 
scientifically described herein. 


MATERIALS AND METHODS 


Meristics, morphometrics and related terminology 
follow explanations given in Kottelat (1990). Termino- 
logy of cephalic sensory system differs in otic canal for 
‘infraorbital canal Al’. Lateral head length excludes 
opercular membrane. Eye diameter is measured horizon- 
tally. Ultimate and penultimate rays of the dorsal and 
anal fin that share the same pterygiophore are counted as 
one. Principal caudal rays are counted upper lobe plus 
lower lobe including unbranched rays. The widely used 
term ‘suborbital flap’ for the protuberance of the lateral 
ethmoid in some species of nemacheilid loach is 
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misleading. This feature is neither located exactly 
suborbital - anteroinferior would be a better fit - nor it is 
proved movable as the term Нар” anticipates. It is 
therefore herein replaced by the term ‘orbital lobe’. 

Measurements are taken point to point with a 
caliper and recorded to nearest 0.1 mm. Regional squa- 
mation densities and scale types were examined using a 
monocular Fenglin XSP-06 (50-200 times magnify- 
cation). 

Lateral and abdominal radiographs (x-rays) of the 
specimen were taken by a Kubtec Xpert 80 and used to 
count vertebrae and fin rays. Vertebral counts and 
associated terminology follow Clothier (1950) and Arratia 
et al (2010); the terminal compound centrum supporting 
the hypural series is counted as one vertebra. In addition, 
radiographs were used to noninvasively measure the width 
of the Weberian apparatus. It is measured radially at the 
center of the gas bladders over the entire complex. 

Location coordinates including altitudes were dete- 
rmined using a Global Positioning System (GPS) Garmin 
handheld device. 

Type series material was preserved in the field using 
a 10% formalin solution and after eight days transferred 
into 90% industrial ethanol. It is placed in the Kunming 
Institute of Zoology (KIZ), Chinese Academy of Scien- 
ces, Kunming, China. 

Morphometrics and meristics of comparative spe- 
cies are from Kottelat (1990) and Zhu (1989). 


Schistura megalodon species nova (Figures 1—6) 


| 
aaa 


Figure 1 Schistura megalodon sp. nov. 
Holotype, KIZ 2014001456, 40.7 mm SL, lateral, dorsal, ventral view. 
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Figure 2 Schistura megalodon sp. nov. 
Holotype, KIZ 2014001456, life coloration. 





Figure 3 Schistura megalodon sp. nov. 
Buccal enclosure, mouth closed, schematic depiction: bc, buccal cavity; ep, 
esophagus; lj, lower jaw; ll, lower lip; pc, palatal curtain; pd, processus 


dentiformes; rt, rostrum; ul, upper lip. 





Figure 4 Schistura megalodon sp. nov. 
Holotype, KIZ 2014001456, head, ventral view, scale bar: 1 mm. 


Type series material 

Holotype: KIZ 2014001456, 40.7 mm SL; Ergan- 
yahe River, Jieyanghe River subbasin, Irrawaddy basin, 
close to Erganya village, East of Tongbiguan Nature 
National Reserve, West of Tongbiguan town, Yingjiang 
County, Dehong Autonomous Prefecture, Yunnan, China; 
N24°36', E97°36', 1323 m above sea level; collected by 
М. Endruweit & ZY. Sha, 2014-01-28. 

Paratypes: KIZ 2014001457, 3 ex., 33-41.7 mm SL; 
KIZ 2014001458 with an opened visceral cavity; 
paratypes collected with the holotype. 


Diagnosis 
Key diagnostic characters of Schistura megalodon 
are a large processus dentiformes in the upper jaw, a 
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Figure 5 Schistura megalodon sp. nov. 
Paratype, KIZ 2014001458, gastrointestinal tract, ventral view, scale bar: 1 


mm. 





Figure 6 Schistura megalodon sp. nov. 
Holotype, KIZ 2014001456, gas bladder complex, ventral view, scale bar: 


1 mm. 


short pre-anus length of 65.4%-66.3% of SL, long paired 
fins (pectoral: 20.8%-24.2% of SL; pelvic: 17.9%- 
20.6% of SL), a wide body of 9.7%—11.3% of SL at anal 
fin origin, an incomplete lateral line, the absence of an 
orbital lobe, and a broad and distinct basicaudal bar with 
forward extensions. 


Description 

Body compact, cylindriform, posteriorly compres- 
sed with a short pre-pelvic length (51.3%-52.1% of SL), 
a short pre-anus length (65.4%—66.3% of SL), and a 
short pre-anal length (75.5%-78.2% of SL); body height 
nearly constant from nape to caudal fin base, slightly 


Kunming Institute of Zoology (CAS), China Zoological Society 


deeper in front of dorsal fin (14.7%—16.7% of SL); body 
wide (at dorsal fin origin 13.2%-15.3% of SL; at A 
origin 9.7%-11.3%). Caudal peduncle deep (1395-1496 
of SL) and long (16.4%—17.5% of SL), its length-depth 
ratio being 1.23—1.33; crests absent. 

Head moderately wide (max. 16.1%-17.5% of SL, 
at nares 10.3%-11.5%: 68.8%-71.6% of lateral HL) and 
slightly depressed (height at eye 11.1%—11.8% of SL; at 
nape 12.6%-14.4% of SL, 91.7%-107% of lateral HL). 
Eye located dorso-laterally, in large specimens (KIZ 
2014001456, 40.7 mm SL; KIZ 2014001457, 41.7 mm 
SL) conspicuously more dorsally; eye small (18.1%- 
18.8% of dorsal HL); interorbital width moderate 
(31%-40% of dorsal HL); ratio interorbital width in eye 
diameter large (1.69—2.12); eyes not visible when viewed 
ventrally. Orbital lobe absent. Nares situated antero- 
dorsally to eye, immediately next to each other; anterior 
nostril is located at the outlet of a valve-like tube. Snout 
blunt; its dorsal profile convex. Two pairs of rostral and 
one pair of maxillary barbels, singly rooted. When 
stretched inner rostral barbel does not reach anterior rim 
of eye, outer rostral barbel reaches midline through eye, 
maxillary barbel exceeds posterior rim of eye. Mouth 
inferior and moderately arched (ratio width to length: 
1.71-2); upper jaw with large, pointed processus 
dentiformes at symphysis; lower jaw gutter shaped and 
notched at symphysis when viewed ventrally. Processus 
dentiformes reaches over lower jaw when mouth is 
closed. A thin, vertical palatal curtain (Figure 3) shortly 
posterior to the tip of the lower jaw arranges a buccal 
enclosure, its distal edge 1s gutter shaped following the 
shape of lower jaw at the given position. Lips moderately 
thick, furrowed; upper lip axially vertically enlarged to 
muffle processus dentiformes, axially slightly notched; 
lower lip notched centrally, not forming lateral triangular 
pads, each side medially with about four lobes separated 
by furrows (Figure 4). Rostral, labial groove continuous 
and very deep; labial groove centrally discontinuous 
around lower lip. 

Body densely scaled, between pectoral fin bases 
sparsely scaled with naked patches. Scales minute, 
cycloid, imbricate. Head not scaled. Lateral line 
incomplete, to below D base; approximately 35—40 
pored scales. Cephalic sensory system well developed; 
supra-temporal and otic canals with three pores each, 
preo-perculo-mandibular canal with 10 pores, 
supraorbital canal with 8—9 pores, and infraorbital canal 
with 10—11 pores. An additional, well visible, short canal 
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of seven pores commences from the upper lip and runs 
nearly parallel to the infraorbital canal before they merge 
at the anterior rim of the eye. Infraorbital canal from 
rostrum to posterior rim of eye and anteriormost three 
pores of preoperculo-mandibular canal situated on 
(fixation-induced) elevated ridges. 

A short esophagus leads into a comparatively large 
stomach; stomach simple, sac-like, u-shaped, multiply 
ridged, nearly horizontally arranged, its posterior wall 
does not reach tip of pectoral fin, ventrally fully muffled 
by a thin layer of spleen, outlet located anteriorly on the 
right side. Intestine straight with two loops in the central 
section; upper loop does not reach wall of stomach 
(Figure 5). Peritoneum hyaline with plenty of silver 
speckles accumulated at visceral arch. 

Vertebral column with 4+(13-14) abdominal and 19 
caudal vertebrae. Weberian apparatus spreads over from 
the atlas to the fourth vertebra; gas bladder complex 
(Figure 6) manubrium shaped with a slight median 
commissure, massive, transversely arranged, located 
superior to the second and third vertebra, wide (ratio to 
interorbital width is 1.2—1.46); gas bladders paired, 
isodiametric, ovoid, ossified; secondary gas bladder 
absent. Urostylic centrum with three lower hypural plates 
directly detached and three upper hypural plates detached 
via a long pleurostylus. All principal caudal fin rays are 
supported via this structure. Dorsally there are 6—7 
rudimentary rays along the posterior part of the caudal 
peduncle associated to neural spines of preural centra 2 
to 4 and one epural; ventrally 4—5 rudimentary rays are 
associated to hemal spines of preural centra 2 and 3. 

Fin formula: D iii, 7; P i, 9; Vi, 67; A iii, 5; C i, 
9+8, i. Paired fins comparatively long (P: 20.8%-24.2% 
of SL; V: 17.9%—20.6% of SL). Pectoral fin horizontally 
inserted, falcate; tip reaches anteriormost tip of ligaments 
of pelvic girdle, but does not reach origin of pelvic fin; in 
KIZ 2014001457 first ray of pectoral fin thickened and 
ventrally flattened. Pelvic fin rays surpass anus by far, 
but do not reach origin of anal fin; origin of pelvic fin at 
vertebra 4-11; 
posteriorly free. Origin of anal fin at vertebra 4421. 


suprapelvic lobe well developed, 
Distal margin of dorsal fin slightly convex; insertion of 
dorsal fin at vertebra 4-11; insertion in larger specimens 
(KIZ 2014001456-7, 40.7—41.7 mm SL) slightly behind 
pelvic fin origin, in smaller specimens (KIZ 20140- 
01458-9, 33-34.1 mm SL) slightly in front of pelvic fin 
origin. Caudal fin slightly emarginated, almost truncate, 
lobes equally long, tips widely rounded. 


Zoological Research 


An overview of selected comparative morphome- 
trics is given in Table 1. 


Coloration of preserved specimens 

Body with 10—12 dark brown bars on light brown 
ground coloration; usually four predorsal, 2-3 subdorsal, 
and 4—5 postdorsal bars; predorsal bars 2—3 times wider, 
more irregular and inconspicuous than subdorsal and 
postdorsal bars. Humeral spot inconspicuous, superim- 
posed to first and second bar. Predorsal bars more 
conspicuous in smaller specimens (KIZ 2014001458-9). 
АП bars interconnected over dorsum, not connected at 
ventral extremity, except last bar between anal fin and 
basicaudal bar. Basicaudal bar conspicuously black, wide, 
straight, with forward directed extensions at both 
extremities; in some specimens (KIZ 2014001457) addit- 
ionally with backward directed extensions. Interspace 
between last body bar and basicaudal bar conspicuously 
light gray. Abdomen in front of anus light beige; 
posterior to anus light gray. Top of head with dark gray 
markings on light gray ground coloration; head ventrally 
light beige. Tip of anterior nostril tube black. Lower lip 
with a peculiar, dark gray mark on each side in larger 
specimens (see Figure 3). 

Dorsal fin with large black anterobasal spot 
reaching onto membrane between first and second 
branched fin ray followed by a dark brown basal band 
over the entire length of the fin base; a median band 
along rays; membranes hyaline. Caudal fin basally with 
black bar as described ranging far onto the fin; a 
transverse gray band at around halfway and another 
slightly lighter band submarginally. Pectoral, pelvic and 
anal fins hyaline with a whitish-yellowish hue. 


Life coloration 

Basically identical to the coloration of preserved 
specimens with differences аз follows: ground 
coloration light brown to purple; the first four pectoral 
fin rays yellow; unbranched and first two unbranched 
dorsal fin rays red distally; unbranched and uppermost 
two unbranched caudal fin rays of the upper lobe 


distally red. 


Etymology 

The specific epithet megalodon is a compound word 
borrowed from ancient Greek: péyac for 'big' and "ободс 
for 'tooth'; an allusion to the large processus dentiformes 
of the upper jaw; a noun in nominative, indeclinable. 
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Table 1 Comparative morphometrics of Schistura megalodon sp. nov. and of related species (values from Kottelat, 1990) 














Schistura megalodon sp. nov. S. malaisei S. vinciguerrae S. sikmaiensis 
types series; n=4 n=10 
Holotype Mean Min Max SD Holotype Mean Min Max SD Neotype 

SL [mm] 40.7 33.0 417 45.9 40.3 65.2 58.7 
In % of SL 

Dorsal HL 22.1 21.6 211 221 041 20 19.5 18.0 20.3 0.76 21.1 
Lateral HL 24.8 23.9 224 248 1.13 24.6 20.3 18.3 21.3 0.87 24.5 
Pre-dorsal length 53.6 532 521 542 0.90 52.3 51.4 48.9 52.6 1.09 53.5 
Pre-pelvic length 514 51.7 513 521 0.43 57.1 51.5 49.7 52.6 0.95 58.8 
Pre-anal length 75.7 76.6 75.5 782 123 79.5 76.1 74.1 77.7 1.00 80.8 
Pre-anus length 66.3 65.9 65.4 663 0.39 74.7 68.8 66.8 713 1.30 77.8 
Head depth (at eye) 11.3 11.4 11.1 11.8 026 11.1 9.9 8.9 10.7 0.50 12.8 
Head depth (at nape) 14.3 13.5 12.6 144 0.96 2.6 11.9 10.7 13.0 0.61 14.7 
Body depth 16.7 15.6 14.7 167 0.92 6.1 18.4 16.1 213 1.69 17.6 
Depth caudal peduncle 14.0 13.4 13.0 140 043 13.5 12.2 11.6 12.6 0.39 12.3 
Length caudal peduncle 17.2 17.0 164 175 0.48 3.3 15.6 14.6 16.2 0.51 12.6 
Snout length 9.3 8.9 8.2 9.4 0.57 8.7 8.9 8.2 9.7 0.45 9.9 
Head width (at nares) 11.5 10.8 103 11.5 0.66 0.7 9.4 8.2 10.1 0.50 11.1 
Maximum head width 17.4 16.8 161 17.5 0.80 15.0 14.1 13.1 14.9 0.48 15.0 
Body width (dorsal origin) 14.7 14.2 13.2 15.3 0.99 2.9 13.7 12.3 15.3 0.99 14.0 
Body width (anal origin) 10.8 10.4 9.7 11.3 0.73 8.9 8.7 7.6 9.7 0.64 9.0 
Eye diameter 42 4.0 3.8 42 0.16 4.6 4.0 3:2 4.5 0.38 3.9 
Interorbital width 8.8 7.8 6.7 8.8 1.10 6.5 7.0 6.3 7.9 0.39 6.1 
Length pelvic fin 17.9 19.2 17.9 206 1.20 16.6 18.7 16.9 20.7 1.02 17.2 
Length pectoral fin 21.4 21.9 208 242 1.59 18.1 22.3 19.6 24.2 1.32 19.3 
Ecology 


The type series was obtained from the Erganyahe 


River (Figure 7). A montane, 


subtropical biome 


predominates this region. At the time of visit (late 
January, 2014, dry season) the stream was approximately 
8-10 m wide and 0.3—0.5 т deep over riffles. Water 
parameter: 12 °C, 20 us/cm, pH 5.9. The water was clear 
and running moderately to fast over boulders, rocks, 
pebbles, sand, and twigs. Aquatic macrophytes were 
absent. Aside Schistura megalodon, there were Schi- 
zothorax malacanthus, Devario kakhienensis, Exostoma 
labiatum, Misgurnus anguillicaudatus, and Channa 
gachua. 

A full day of fishing just yielded the four specimens 
that compose the type series. Schistura megalodon is 
regarded as rare. This is in compliance with the 
statements of locals and Shuwei Liu (KIZ), who spend a 
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Figure 7 Erganyahe River, type locality of Schistura megalo- 
don sp. nov., view northwards, downstream. 


full spend a full day fishing at the Erganyahe River close 
to Erganya village in June 2013. He collected just one 
specimen that was identified conspecific basing on 
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photographs. The specimen concerned was not preserved 
(SW. Liu, pers. comm.). 

The large processus dentiformes may be a trophic 
adaptation. It is hypothesized that it may be used as a 
tool to detach eggs and other food items from the 
substratum or it may be used for efficient digging. It is 
equally well developed in the type specimens and 
therefore most likely not a sexual dimorphic character. 

The dorsal fin migrates backwards during ontogeny 
in Schistura megalodon. Its position relatively to the 
pelvic fin, thus, proves unsuitable as a diagnostic 
character. Smaller specimens of S. megalodon have the 
dorsal fin inserted slightly in front of the pelvic fin, while 
larger specimens have it inserted slightly behind the 
pelvic fin origin. This is manifested by the predorsal 
length: the two smaller paratypes (KIZ 2014001458—9, 
33-34.1 mm SL) have a predorsal length of 52.1%- 
52.8 mm SL, the holotype (KIZ 2014001456, 40.7 mm 
SL) and the largest paratype (KIZ 2014001457, 41.7 mm 
SL) have а predorsal length of 53.6%-54.2% of SL. 


Distribution 

Schistura megalodon is only known from the type 
locality (Figure 8). A few kilometers downstream of the 
type locality the Erganyahe River marks the borderline 
between China and Myanmar flowing northwards prior 
to its influx in the Jieyanghe River. The new species is, 
thus, expected to occur on Burmese territory as well. 


DISCUSSION 


Certainly Schistura megalodon has been hitherto 
overlooked by taxonomists working on this group of 
fishes due to its diminutive size, its rather limited 
distributional range, its occurrence close in the politically 
sensitive border region of China and Myanmar and the 
therefore restricted accessibility, and its absolute scarcity 
in the habitat. 

The palatal curtain in Schistura megalodon sepa- 
rates the pharynx into an anteriorly located prebuccal 
cavity and a posteriorly located buccal chamber. The 
buccal chamber is quasi sealed by the palatal curtain 
when the mouth is closed. This morphology reminds of 
the condition found in member of the genus Neono- 
emacheilus Zhu & Guo, 1985. However, in Neonoema- 
cheilus the prebuccal cavity is formed by hypertrophied 
lips which is not the case in S. megalodon. Furthermore, 
S. megalodon does not possess characters diagnostic of 
Neonoemacheilus such as an orbital lobe in males and a 
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Figure 8 Schistura megalodon sp. nov., type locality 


complete lateral line. It perfectly fits into the generic 
description of Schistura as provided by Kottelat (1990: 
24) and is therefore herein treated as a member thereof. 

Phenotypically close to Schistura megalodon 1s S. 
cincticauda (Blyth, 1860) from the Salween basin in 
Thailand. Schistura megalodon differs from S. cincti- 
cauda in having a smaller size (41.7 mm SL vs. 50 mm 
SL); a shorter pre-dorsal length (52.1%-54.2% of SL vs. 
54.7%-58.8%); a longer caudal peduncle (16.4%—17.5% 
of SL vs. 12.2%-16.4%); a wider body at anal fin origin 
(9.7%-11.3% of SL vs. 7.7%—9.4%; 45.1%—52.2% of 
dorsal HL vs. 34%-45%); a smaller eye diameter 
(3.8%-4.2% of SL vs. 4.6%-5.7%; 18.1%—18.9% of 
dorsal HL vs. 21%2—7%); suprapelvic lobe present vs. 
absent; and further differences in the color pattern. 

Schistura megalodon is distinguished from 5. 
malaisei by a shorter pre-pelvic length (51.3%-52.1% of 
SL vs. 57.1% of SL in S. malaisei); a shorter pre-anus 
length (65.4%-66.3% of SL vs. 71.1%-74.7%); a longer 
caudal peduncle (16.4%-17.5% of SL vs. 11.8%—13.3%), 
a wider head (16.1%—17.5% of SL vs. 15%-15.1%); а 
wider body (at dorsal fin origin: 13.2%—15.3% of SL vs. 
12.9%; at anal fin origin: 9.7%-11.3% of SL vs. 
8.0%-9.6%); a longer pelvic fin (17.9%-20.6% of SL vs. 
14.7%-16.6%); a longer pectoral fin (20.8%-24.2% of 
SL vs. 17.6%—18.1%); and a different color pattern. 
Schistura malaisei was described from Putao in northern 
Myanmar. 

Another similar nemacheilid loach described from 
the Irrawaddy basin in the Manipur State in India is S. 
sikmaiensis. Schistura megalodon 1s readily distinguish- 
able from S. sikmaiensis by having iii, 7 dorsal fin rays 
vs. iv, 8; 1, 9 pectoral fin rays vs. 11; lateral incomplete 
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with 35-40 pored scales vs. complete with 100—110 
pored scales; a shorter pre-pelvic length (51.3%-52.1% 
vs. 58.8%); a shorter pre-anus length (65.4%-66.3% of 
SL vs. 77.8%); a shorter pre-anal fin (75.5%—78.2% vs. 
80.8%); a less deep head (at eye: 11.1%-11.8% vs. 
12.896; at паре: 12.6%—14.4% vs. 14.7%); a less deep 
body (14.7%-16.7% vs. 17.6%); a deeper and longer 
caudal peduncle (depth: 1396-1496 vs. 12.3%; length: 
16.4%-17.5% vs. 12.6%); a shorter snout (8.2%—9.2% 
vs. 9.9%); a wider head (16.1%-17.5% vs. 15%); a 
wider body at anal fin origin (9.7%—11.3% vs. 9%); a 
wider interorbital distance (6.7%—8.8% vs. 6.1%); longer 
paired fins (pelvic fin: 17.9%—20.6% vs. 17.2%; pectoral 
fin: 20.8%—24.2% vs. 19.3%); processus dentiformes 
present vs. absent; and significant differences in color 
pattern. Kottelat (1990:210) designates the holotype of 
Nemacheilus putaoensis as neotype of S. sikmaiensis 
which makes N. putaoensis an objective junior synonym 
to S. sikmaiensis. Zhu (1989:62, Figure 42) lists one 
specimen of N. putaoensis from the Dayingjiang River in 
Yingjiang County and gives the species' distributional 
range as Jieyanghe River, the same subbasin from which 
S. megalodon is described. This river is an independent 
tier one tributary to the Irrawaddy itself and does not 
belong to the Dayingjiang subbasin. Zhu places 
putaoensis in Nemacheilus which actually implies the 
presence of an orbital lobe in males according to his 
generic diagnosis (1989:60). Both, the description and 
the accompanying figure, however, are lacking an orbital 
lobe. Zhu does not mention a processus dentiformes 
which fits well to the species' diagnosis. 

Schistura vinciguerrae (Figure 9a) was described 
from Meekalan within the Tenasserim mountain range in 
southern Myanmar and is reported from the Irrawaddy 
basin Yingjiang County in Dehong by Kottelat (1990:219) 
and Zhu (1989:52). Schistura megalodon is distinguished 
from S. vinciguerrae by having iii, 7 dorsal fin rays vs. iv, 
8-9; 1, 9 pectoral fin rays vs. 11-12 in total; lateral 
incomplete with 35—40 pored scales vs. complete with 
84-90 pored scales; a longer head length (dorsally: 
21.1%-22.1% of SL vs. 18%-20.3%; laterally: 
22.4%-24,8% vs. 18.3%-21.3%); a deeper and wider 
head (depth at eye: 11.1%-11.8% vs. 8.9%—10.7%; max. 
width: 16.1%-17.5% vs. 13.1%—14.9%); a shorter pre- 
anus length (65.4%—66.3% of SL vs. 66.8%—71.3%); a 
deeper and longer caudal peduncle (depth: 13%—14% vs. 
11.6%-12.6%; length: 16.4%—17.5% vs. 14.6%-16.2%); 
a wider body at anal fin (9.796—11.39 vs. 7.6%-9.7%), a 


Kunming Institute of Zoology (CAS), China Zoological Society 


smaller maximum size of 41.7 mm SL vs. 71 mm SL; 
and a far different color pattern. 





Figure 9 a. Schistura vinciguerrae, KIZ2012004064, 42.6 mm 
SL, Tanjiazhai, Tengchong County, Longchuanjiang River 
subbasin; b. Schistura polytaenia, К172012004077, 49 mm SL, 
male, Mahei village, Tengchong County, Dayingjiang River 
subbasin. 


Zhu (1982:107) describes Schistura yingjiangensis 
(as Nemacheilus yingjiangensis) from the Dayingjiang 
River in Yingjiang County and gives its distributional 
range as Dayingjiang апа Longchuanjiang rivers 
(1989:65). Schistura megalodon is distinguished from 5. 
yingjiangensis by 1, 9 pectoral fin rays vs. i, 10; lateral 
line incomplete vs. complete; orbital lobe absent vs. 
present; pelvic clearly surpasses anus vs. does not 
surpass; smaller maximum size (41.7 mm SL vs. 61 mm 
SL); and a longer pre-dorsal length (52.1%-54.2% of SL 
vs. 44%-50%). In S. megalodon the lateral HL is 
4.03—4.46 times in SL vs. 4.5—5 times; the eye diameter 
1.69-2.12 times in interorbital width vs. 1.4—1.5 times; 
the eye diameter 5.69—6.15 times in lateral HL vs. 4.2-5 
times; and interorbital width 2.81—3.36 times 1n lateral 
HL vs. 3.7-4 times. Besides, there are significant 
differences in the color pattern. 

Schistura megalodon is easily distinguished from 5. 
polytaenia (Figure 9b) described from the Long- 
chuanjiang River in having constant 17 (9+8) branched 
caudal fin rays vs. 15—16; vertebrae 4+(32—33) vs. 
4+(29—31); lateral line incomplete vs. complete; orbital 
lobe in males absent vs. present; pelvic fin clearly 
surpasses anus vs. does not reach; a broad and deeply 
black, conspicuous with forward 
extensions ай its 


basicaudal bar 
extremities vs. short, brown, 
inconspicuous, does not reach extremities; and a body 
pattern of 10—12 bars interconnected over the dorsum vs. 


18—30 bars and 9—13 saddles along the dorsum. In fact, 
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the color pattern along the flanks of S. polytaenia is 
highly irregular. It may be more appropriate to define it 
as mottled rather than barred. This is particularly true in 
specimens from the Mengjiahe River іп Sudian, 
Yingjiang County. 

Another nemacheilid species originally described 
from upper reaches of the Chaungmagyi River, a left 
bank, first tier tributary to the Irrawaddy, is the enigmatic 
Nemacheilus (described as Nemachilus) acuticephalus. 
Apparently, this species has not been collected since its 
original description (Hora, 1929). All subsequent authors 
1990:97; Talwar & Jhingran, 
1992:468) refer to this description. Schistura megalodon 


(inter alia: Kottelat, 


is distinguished from N. acuticephalus in having a broad, 
straight and distinct basicaudal bar with forward 
extensions and its extremities vs. a crescent-shaped and 
pale basicaudal bar not reaching extremities; 10 pectoral 
fin rays vs. 9; 9+8 branched caudal fin rays vs. 8+8 
(taken from Hora, plate 14; see Kottelat); pelvic fin 
surpasses anus vs. not reaching anus; crests on caudal 
peduncle absent vs. prominent crests supported by 
rudimentary rays on the dorsal and ventral extremity of 
the caudal peduncle; caudal fin slightly emarginate vs. 
forked. Besides, there are differences in the color pattern. 
Kottelat placed N. acuticephalus in Schistura. This 
generic allocation seems to be incorrect. The prominent 
caudal peduncle crests, which are in N. acuticephalus at 
least dorsally widely supported by ргосштепі 
rudimentary rays, are diagnostic for members of 
Homatula. Its mottled color pattern is also typical of 
members of Homatula. However, its small size of just 40 
mm TL and the hitherto absence of this genus in the 
Irrawaddy basin act contradictorily. In order to finally 
access the generic allocation of N. acuticephalus an in 
depth examination of either the type series or a large batch 
of topotypic fresh material is recommend. Until this 
question is satisfactorily resolved the taxon is tentatively 
treated as S. acuticephala following Kottelat (2012:104). 
Key to the species of the genus Schistura reported 
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Influence of a large dam on the longitudinal patterns of fish assem- 
blages in Qingyi Stream 
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Abstract: Using seasonally collected data (2009—2010) from 15 sampling sites that represent first- to fifth-order streams within 
the Qingyi watershed, we examined the spatio-temporal patterns of fish assemblages along two longitudinal gradients to explore 
the effects of a large dam on fish assemblages at the watershed scale. No significant variation was observed in either species 
richness or assemblage structure across seasons. Species richness significantly varied according to stream order and gradient. 
Dam construction appeared to decrease species richness upstream substantially, while a significant decrease between gradients 
only occurred within fourth-order streams. Along the gradient without the large dam, fish assemblage structures presented 
distinct separation between two neighboring stream orders, with the exception of fourth-order versus fifth-order streams. 
However, the gradient disrupted by a large dam displayed the opposite pattern in the spatial variation of fish assemblages 
related with stream orders. Significant between-gradient differences in fish assemblage structures were only observed within 
fourth-order streams. Species distributions were determined by local habitat environmental factors, including elevation, substrate, 
water depth, current discharge, wetted width, and conductivity. Our results suggested that dam construction might alter the 
longitudinal pattern in fish species richness and assemblage structure in Qingyi Stream, despite the localized nature of the 


ecological effect of dams. 


Keywords: Fish assemblage; Species richness; Spatiotemporal pattern; Longitudinal gradient; Dam building 


Stream fish assemblages are structured by abiotic 
factors, biotic interactions, and historical processes 
(Hoeinghaus et al, 2007). The physicochemical stream 
environment exhibits spatial heterogeneity and temporal 
variability, and may cause spatiotemporal variations in 
fish assemblages (Meador & Matthews, 1992). Typically, 
among the multiple spatial factors affecting fish 
assemblages, longitudinal-gradient variations in environ- 
mental conditions substantially influence the distribution 
and abundance of stream fishes (Araújo et al, 2009; 
Inoue & Nunokawa, 2002; Taniguchi et al, 1998; Vann- 
ote et al, 1980). From headwaters to downstream, fish 
species richness generally increases due to increases in 
stream size (Matthews, 1986), habitat diversity (Gorman 
& Karr, 1978), and shelter availability (da Silva Abes & 
Agostinho, 2001), and the different rates of fish immig- 
ration and extinction (Power et al, 1988). Concurrently, 
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fish species composition may vary longitudinally by 
species addition and/or species replacement (Boys & 
Thoms, 2006; Gorman & Karr, 1978; Roberts & Hitt, 
2010). The River Continuum Concept (RCC) describes 
the changes in community structure and species richness 
of organisms from headwaters to mouth waters and 
relates these changes to flow regime, water temperature, 
food availability, and substrate conditions (Vannote et al, 
1980). However, due to various human activities (e.g., 
damming, pollution, erosion), few riverine ecosystems 
remain free-flowing over their entire course. Overlaying 
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this pattern displayed by the RCC, the Serial Discon- 
tinuity Concept (SDC) claims that regulation by dams 
disrupts the underlying continuum and produces a series 
of lentic and lotic reaches, thereafter causing longitudinal 
shifts in abiotic and biotic parameters and processes 
(Ward & Stanford, 1983). 

Dams impact stream fishes and invertebrates in div- 
erse ways, such as blocking migratory pathways (March 
et al, 2003), fragmentizing habitats (Travnichek et al, 
1993), altering natural flow regimes (Bonner & Wilde, 
2000) and food webs (Power et al, 1996), decreasing 
water temperature downstream (Clarkson & Childs, 2000) 
and current velocity upstream (Bennett et al, 2002), 
disrupting riparian plant communities (Nilsson et al, 
1997), and shifting water chemistry (Humborg et al, 
1997). The extent to which stream fishes and invert- 
ebrates are affected by dams may be associated with the 
characteristics of dams (e.g., location, purpose, and 
management) and fauna (Cumming, 2004; March et al, 
2003). For example, large dams without spillways are 
impermeable barriers for migratory organisms and may 
extirpate all native migratory fish from upstream habitat 
(Holmquist et al, 1998). In contrast, large dams with 
spillways provide a possible passage for some native 
fishes. To date, investigations have primarily identified 
the effects of dams on lotic reaches directly below dams, 
mainstream reservoirs directly above dams, and lotic 
reaches upstream of impoundments (Clarkson & Child, 
2000; Cumming, 2004; Santucci et al, 2005; Travnichek 
et al, 1993). However, little attention has been paid to the 
effects of dams at the watershed scale. 

Stream ecosystems in China are experiencing mas- 
sive ecological perturbation due to diverse human 
activities, such as extensive agricultural and industrial 
production, urbanization, development, and hydraulic 
engineering construction, all of which threaten freshw- 
ater fish (Chen, 2005). Based on data collected from 
first- to fifth-order streams, representing two longitu- 
dinal gradients (one disrupted by a large dam and the 
other not) within a tributary (Qingyi Stream) of the 
lower reaches of the Yangtze River, the spatial and tem- 
poral patterns of fish assemblages were examined. Tota- 
lly, our three key goals were to: (1) to determine how 
fish assemblage patterns vary longitudinally from first- to 
fifth-order streams; (2) to identify the correlations 
between local habitat and fish assemblages; and (3) to 
assess how a large dam affects the longitudinal patterns 
of fish assemblages at the watershed scale. 
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MATERIALS AND METHODS 


Study area 

Qingyi Stream originates from the northern portion 
of Huangshan Mountain, and flows northeast toward its 
confluence with the lower reaches of the Yangtze River. 
This watershed is approximately 309 km long and covers 
an area of 7195 km’. Due to the subtropical monsoon 
climate, temperature and precipitation in this area are 
quite asymmetric across seasons. Annual air temperature 
ranges from —2.1 °C (January) to 27.5 °C (July), with a 
mean of 17.8 °C. Annual precipitation is high (approxi- 
mating 2000 mm/year), with most rainfall (79%) occur- 
ring in spring and summer (April to September) and only 
minor rainfall (less than 5%) occurring in the cold, dry 
winter (December to February) (Yan et al, 2011). 

Chencun Hydropower Station was constructed at the 
mid-reaches of the Qingyi mainstream in the 1950s. It is 
the largest station along the lower reaches of the Yangtze 
River, with 119 m deep impoundment and covering an 
area of 88.6 km’. Derived from Anhui Province 
topographic maps (1:300 000 scale), the Qingyi Stream 
is categorized as fifth-order (Strahler, 1957) at its 
mainstream. A total of 75 first-order, 34 second-order, 12 
third-order, 4 fourth-order, and 1 fifth-order stream are 
included in this watershed. Among the fourth-order 
streams, the Qingxi, Shuxi, and Machuan streams flow 
into an artificial reservoir (1.e., Taiping Lake) formed by 
the impoundment of Chencun Hydropower Station, while 
Huishui Stream directly flows into the mainstream with 
its confluence downstream of the station (Figure 1). 


Fish sampling 

A total of 15 segments representing first- to fifth- 
order streams were surveyed seasonally during April, 
August, and October 2009, and January 2010. The 
sampling segments were derived from two longitudinal 
gradients. One gradient (A) was set from the headwaters 
of Qingxi and Shuxi streams to the lower reaches of the 
mainstream, representing the gradient disrupted by 
Chencun Hydropower Station, and the other gradient (B) 
was set from the headwaters of Huishui Stream to the 
mainstream, representing the gradient not disrupted by a 
station. Each surveyed segment was set on the 
topographic map based on their stream orders. Four 
segments ranging from first- to fourth-order streams 
were selected within Qingxi, Shuxi, and Huishui streams, 
respectively. Three other segments representing fifth-order 
streams were set at the mainstream (Figure 1). Sampling 
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Yangtze 






Anhui Province 


Figure 1 Map of the Qingyi Stream watershed in Anhui Province, China 


Solid circles mark the 15 sampling sites. Sites 1-3, 4-6, 7—9, 10-12, and 13—15 represent first-, second-, third, fourth-, and fifth-order streams, respectively. 


Sites 1, 2, and 3 are located at Qingxi (QS), Shuxi (SS), and Huishui Streams (HS), respectively. у-< shows Chencun Hydropower Station. 


sites were selected in the field based on habitat represe- 
ntativeness and accessibility. Each shallow-water (<1.0 m 
depth) site of first-, second-, and third-order streams was 
sampled using backpack electro-fishing by wading in two 
passes. Each site encompassed at least two mesohabitat 
units (pool and riffle). Each deep-water (71.0 m depth) site 
of fourth- and fifth-order streams was sampled twice using 
boat electro-fishing along both riversides. Sampling sites 
were far (721.0 km distance) from conspicuous human 
disturbance, such as dam, farmland, and urban land, and 
each sampling was conducted with comparable effort (1.е., 
100 m long and approximately 40 min). Fish were 
identified to species level (except for Ctenogobius due to 
the deficiency in identifying tools), counted, and returned 
to sampling sites if alive. Voucher specimens were placed 
in 8% formaldehyde solution for further identification. 


Environmental survey 
Each sampling site was characterized using 10 


variables to describe local habitat: elevation (m), wetted 
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width (m), water depth (m), water temperature (°С), pH, 
conductivity (mS/cm), dissolved oxygen (mg/L), current 
velocity (m/s), discharge (m/s), and substrate. Elevation 
was determined by a portable GPS receiver. Wetted 
width was measured along five transects, regularly 
spaced across the stream channel. Water depth, temperat- 
ure, dissolved oxygen, pH, and conductivity were 
surveyed at four equal-interval points along each transect. 
Water depths «1.0 m were measured using a graduated 
wading rod, and depths 71.0 m were measured using a 
supersonic echo sounder. Current velocity was taken at 
60% of water depth at each point of each transect. 
Discharge of each channel was determined at the transect 
that yielded the most accurate measurement (smooth 
bottom and laminar flow). Along each transect, the 
proportion of substrate categories (particle size 1= 0-1 mm; 
2-1-5 mm; 325-25 mm; 4-25-50 mm; 5-50-100 mm; 
6-100-500 mm; 7-500-1 000 mm; 8->1 000 mm) 
was visually estimated, and an index of substrate coars- 
eness ranging from 1 to 8 was derived for each site 
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following Bain et al (1985). 


Data analysis 

Each fish species collected in this study was 
examined for frequency of occurrence (F) and relative 
abundance (RA), estimated from F;-100(S/S)?o and 
P=100(n/n)%, respectively, where 5; and S are the 
abundance of the samplings of which species i were 
collected and of total samplings, respectively, and n; and 
n are the individual numbers of species i and total fish 
species, respectively. The above F and P values were 
determined independently for three groups of sampling 
segments, i.e., first- to fourth-order streams along 
gradient A and B and fifth-order streams, respectively. 
Species richness and fish abundance were analyzed using 
a separate three-way ANOVA model with stream order 
(first- to fifth-order), season (January, April, August, 
October), and dam (gradient A and B) as factors. One- 
way ANOVA was used to test differences in species 
richness across stream orders along gradients A and B, 
respectively. The Newman-Keuls test was used for post- 
hoc comparisons after ANOVA. Independent sample t- 
tests were used to compare species richness within the 
same ordered streams between the two gradients. All data 
were log-transformed to meet the assumptions of 
normality and homogeneity of variances. 

Discrete spatial patterns in fish assemblages were 
identified using PRIMER 5 (Primer-E Ltd., 2001). 
Following a Bray-Curtis similarity matrix calculation, an 
analysis of similarity (ANOSIM) was used to test 
variations in fish assemblages across stream orders, 
seasons, and gradients. Firstly, variations associated with 
stream orders and seasons were analyzed separately for 
gradients A (two-way crossed ANOSIM) and B (one-way 
ANOSIM due to no replication), respectively. Secondly, 
variations associated with gradients and seasons were 
analyzed separately for each stream order using two-way 
crossed ANOSIM. The relationships among assemblages 
from each site and season were graphically presented 
using non-metric multi-dimensional scaling (NMS) 
analysis. The contribution of each species to differences 
among assemblage groups was identified using similarity 
percentages analysis (SIMPER) (Clarke & Warwick, 
2001). The variability in fish assemblages in relation to 
local habitat environment was evaluated by canonical 
correspondence analysis (CCA) using CANOCA 4.5. All 
variables entered the CCA after a forward-selection 
procedure, showing their importance in explaining total 
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variability in species composition. The significance 
(P<0.05) of the CCA gradients was assessed by Monte 
Carlo permutation tests, and their importance measured 
by the eigenvalues of the first two axes (ter Braak & 
Verdonschot, 1995). All variables of fish assemblages 
(1.е., species richness and abundance) and habitat 
log(X+1) transformed 
assumptions of multivariate normality and to moderate the 


environment were to meet 
influence of extreme data. Species that occurred in less 
than two sites were excluded from the above analysis to 


avoid negligible weighting (Gauch, 1982). 


RESULTS 


Overview of species diversity 

A total of 13 647 fish were captured throughout this 
study, representing 57 species, 15 families, and 5 orders. 
Species of family Cyprinidae comprised on average 
59.6% of total species richness Species richness per 








sampling site amounted to 11.9+6.6 (mean+SD) species, 





and abundance per site was 227.6+205.5 specimens. A 
total of 16, 24, 24, 47, and 44 species were collected in 
first- to fifth-order streams, respectively. Cyprinus carpio, 


Mylopharyngodon piceus, | Hemiculter  leucisculus, 
Sarcocheilichthys | sinensis, Abbottina obtusirostris, 
Macropodus chinensis, Mystus macropterus апа 


Hyporhamphus intermedius were only collected within 
fifth-order streams. Within the first- to fourth-order 
streams, 18 species were only collected from streams 
along gradient B, including Distoechodon tumirostris, 
Elopichthys bambusa, and Spinibarbus hollandi. Sauro- 
gobio dabryi was a unique species occurring in gradient 
A but not gradient B. Zacco platypus, Acrossocheilus 
fasciatus, Pseudogobio vaillanti, Vanmanenia steno- 
soma and Ctenogobius sp. were common (>40% of F) 
and relatively abundant (>1% of P) within both gradi- 
ents. In addition, Acheilognathus barbatulus, A. chan- 
kaensis, Rhodeus ocellatus, Squalidus argentatus and 
Misgurnus anguillicaudatus were more frequent and 
abundant within gradient A, while Opsariichthys bidens, 
Hemiculter bleekeri, Sarcocheilichthys parvus, Parab- 
otia fasciata, Siniperca chuatsi, Odontobutis obs- 
curus, Silurus asotus and Pseudobagrus albomargin- 
atus were more frequent and abundant in gradient B 
(Table 1). 


Local species richness 
Fish species richness significantly varied according 
to stream order and gradient, but not season (Table 2). 
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Table 1 Frequency of occurrence (%; left of dash) and relative abundance (%; right of dash) of fish collected within Qingyi 








Stream 
First to fourth-order Fifth-order 
Species Code 
Gradient A Gradient B 
Cypriniformes 
Abbottina rivularis ABR 50.0/1.7 31.3/2.3 83.3/2.9 
Abbottina tafangensis ABT* 3.1/<0.1 6.3/0.1 
Abbottina obtusirostris ABO* 16.7/0.2 
Acheilognathus barbatulus ACB 53.1/3.9 12.5/0.3 58.3/0.8 
Acheilognathus chankaensis ACC 18.8/0.9 6.3/0.3 41.7/0.5 
Acrossocheilus fasciatus ACF 68.8/7.7 87.5/10.1 25.0/0.3 
Carassius auratus CAA 28.1/1.8 31.3/2.1 91.7/5.8 
Cobitis rarus COR 53.1/5.6 37.5/2.2 
Cobitis sinensis COS 28.1/0.6 18.8/1.3 8.3/«0.1 
Cyprinus carpio CYC* 16.7/0.1 
Culter erythropterus CUE* 6.3/«0.1 
Distoechodon tumirostris DIT* 6.3/0.1 
Elopichthys bambusa ELB* 6.3/«0.1 
Erythroculter ilishaeformis ERI 6.3/0.2 41.7/2.3 
Gnathopogon taeniellus GNT 18.8/0.6 25.0/0.4 
Gobiobotia tungi GOT* 12.5/0.1 
Hemiculter leucisculus HEL 75.0/7.9 
Hemiculter bleekeri HEB 6.3/0.2 31.3/4.3 58.3/1.2 
Hemibarbus labeo HEL 18.8/0.4 
Hemibarbus maculatus HEM 31.3/0.7 16.7/0.1 
Misgurnus anguillicaudatus MIA 37.5/1.1 18.8/0.4 8.3/«0.1 
Mylopharyngodon piceus MYP 41.7/0.8 
Opsarüchthys bidens OPB 21.9/1.0 56.3/1.6 66.7/2.1 
Phoxinus oxycephalus PHO 12.5/0.2 18.8/3.8 
Pseudobrama simoni PSS 18.8/1.6 58.3/3.8 
Rhodeus ocellatus RHO 65.6/13.9 31.3/2.3 66.7/3.8 
Spinibarbus hollandi SPH* 12.5/0.1 
Zacco platypus ZAP 100/34.7 93.8/26.0 58.3/5.5 
Parabramis pekinensis PAP 6.3/0.1 16.7/0.1 
Parabotia fasciata PAF 18.8/0.6 31.3/2.0 33.3/5.0 
Pseudorasbora parva PSP 15.6/0.6 25.0/2.7 50.0/2.2 
Pseudogobio vaillanti PSV 43.8/2.5 43.8/4.8 33.3/0.5 
Sarcocheilichthys sinensis SAS* 8.3/«0.1 
Sarcocheilichthys parvus SAP 3.1/0.1 25.0/2.9 66.7/2.8 
Sarcocheilichthys nigripinnis SAN 25.0/4.0 75.0/2.3 
Squalidus argentatus SQA 344/1.2 12.5/0.7 33.3/1.0 
Saurogobio dabryi SAD 6.3/0.1 8.3/«0.1 
Varicorhinus barbatulus VAB* 12.5/0.2 
Vanmanenia stenosoma VAS 65.6/5.2 50.0/3.9 25.0/1.1 
Perciformes 
Channa argus CHA 18.8/1.4 91.7/6.1 
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(Continued) 
First to fourth-order Fifth-order 
Species Code 
Gradient A Gradient B 
Channa asiatica CHA* 6.3/0.1 8.3/0.1 
Siniperea roulei COR 12.5/0.3 33.3/0.3 
Ctenogobius sp. CTS 93.8/12.6 68.8/4.4 25.0/2.0 
Hypseleotris swinhonis HYS 12.5/0.2 9.3/0.1 
Mastacembelus aculeatus MAA 21.9/0.6 6.3/0.2 
Macropodus chinensis MAC* 16.7/«0.1 
Odontobutis obscurus ODO 12.5/0.9 37.5/2.7 91.7/7.9 
Siniperca chuatsi SIC 3.1/«0.1 25.0/1.2 75.0/2.4 
Siniperca obscura SIO 6.3/0.2 50.0/0.6 
Siluriformes 
Glyptothorax fukiensis GLF* 6.3/0.9 
Liobagrus styani LIS 31.3/1.0 31.3/0.6 
Mystus macropterus MYM* 8.3/0.1 
Pelteobagrus fulvidraco PEF 25.0/3.9 100.0/21.1 
Pseudobagrus albomarginatus PSA 3.1/«0.1 37.5/0.9 66.7/4.1 
Silurus asotus SIA 3.1/0.1 25.0/0.5 41.7/0.7 
Synbranchiformes 
Monopterus albus MOA 31.3/0.5 18.8/0.4 33.3/0.9 
Beloniformes 
Hyporhamphus intermedius HEI* 8.3/0.9 


*: Rare species occurring in less than two sites, not included in statistical analysis. 


Table2 Three-way ANOVA results of the spatial and temporal variations in fish species richness within Qingyi Stream 





Factors df SS MS F P Student-Newman-Keuls 

Stream orders 4 2.3 0.6 38.9 ЖЖ Ist<2nd=3rd<4th<Sth 
Seasons 5 0.1 0.02 1.9 ns 

Gradients 1 0.1 0.1 4.0 ® 

Ordersxseasons 12 0.2 0.01 0.9 ns 

Seasonsx gradients 3 0.02 0.01 0.5 ns A<B 
Orders*gradients 4 0.6 0.1 9.5 Er 

Огаегѕхѕеаѕопѕхргайіепіѕ 12 0.1 0.01 0.8 ns 


Factors are stream orders (first-fifth), seasons (January, April, July, October) and gradient (A, B; with and without the effect of a large dam, respectively). 


ж; Р<0.05, **: P<0.01, ns P>0.05. Student-Newman-Keuls test: P<0.05. 


Species richness significantly increased with stream 
orders, with the exception of no significant difference 
between second- and third-order streams. Gradient B 
had more species than that of Gradient A. A 
significant interaction effect was observed for stream 
order by gradient (Table 2). Within Gradient A, the 
lowest and highest species richness occurred in first- 
order and fifth-order streams, respectively, and species 
richness did not vary significantly among second-, 
third- and fourth-order streams. Within Gradient B, 
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however, the highest species richness was observed in 
fourth-order streams, while first-, second- and third- 
order streams were not significantly different in 
species richness (Figure 2). When comparing species 
richness within the same-order streams between the 
two gradients, significant differences were only 
observed in fourth-order streams (t-test, F=0.23, 
P«0.01), but not in first- (F=0.4, P=0.54), second- 
(Е-1.5, P-0.21) or third-order streams (F=7.2, 
P=0.19) streams. 
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Figure 2 Fish species richness along two longitudinal gradients 
by stream orders within Qingyi Stream (upper: gradient A; 
lower: gradient B) 

Different lowercases represent significant differences in species richness 

(One-way ANOVA, P<0.05) . 


Fish assemblages 

The two-way crossed ANOSIM results suggested 
that fish assemblages significantly differed across stream 
orders but not seasons, and were uniform for both 
longitudinal gradients of A and B (Table 3). Within 
gradient B, significant differences in inter-order fish 
assemblages were almost exclusively observed, with the 
exception of fourth- vs. fifth-order. Whereas, fish assem- 
blages in fourth-order streams were similar with those in 


Table 3 R-values and their significance levels (P) for 
comparisons of fish assemblages among stream orders and 








seasons using ANOSIM 
Gradient A Gradient B 
R P R P 

Order 0.50 0.001 0.66 0.001 
Ist vs. 2nd" 0.06 0.440 0.60 0.029 
Ist vs. 3rd 0.38 0.012 0.42 0.036 
Ist vs. 4th 0.29 0.049 1.00 0.029 
Ist vs. 5th 1.00 0.001 0.99 0.001 
2nd vs. 3rd 0.06 0.420 0.80 0.029 
2nd vs. 4th” 0.06 0.370 1.00 0.029 
2nd vs. 5th 0.98 0.001 0.96 0.002 
3rd vs. 4th" 0.19 0.296 1.00 0.029 
3rd vs. 5th 0.90 0.001 0.93 0.002 
4th vs. Sth” 0.85 0.002 0.06 0.310 
Season -0.25 0.988 -0.10 0.971 


One-way ANOSIM was conducted for gradient (B) without the effect of 
Chencun Hydropower Station, two-way crossed ANOSIM was conducted 
for gradient (A) disrupted by the station and both gradients (Both). 
Significance level was accepted at P<0.05 (in bold). *: opposite significance 


level occurring between the two gradients. 
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second- and third-order but not fifth-order streams within 
gradient A (Table 3). For fish assemblages within the 
same-order streams, significant between-gradient differe- 
nces were observed in fourth-order streams, but not in 
first-, second, and third-order streams (Table 4). 


Table 4 R-values and their significance levels for compari- 
sons of fish assemblages among seasons and between two 
gradients disrupted by a large dam or not using ANOSIM 








R P 
Season Gradient Season Gradient 
Ist-order —0.33 -0.25 0.876 0.815 
2nd-order —0.04 0.00 0.562 0.617 
3rd-order -0.31 -0.73 0.800 0.988 
4th-order —0.38 0.75 0.790 0.020 
5th-order -0.19 0.872 


Two-way crossed ANOSIM was conducted for first- to fourth-order streams, 
one-way ANOSIM was conducted for fifth-order streams. Significance level 
was accepted at P<0.05 (in bold). 


Based on NMS analysis, fish assemblages within 
gradient B varied substantially from headwaters to main- 
stream, despite the overlap between fourth- and fifth- 
order streams (Figure 3). In contrast, within gradient A, 
substantial overlaps were observed between neighboring 
orders from first- to fourth-order streams, while fish 
assemblages in fourth-order streams were distinct from 
those in fifth-order streams. 
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Figure 3 Non-metric multi-dimensional scaling (NMS) by stre- 
am orders and dams for fish assemblage data in Qingyi Stream 
Black symbols, gradient disrupted by Chencun Hydropower Station; Open 
symbols, gradient without damming; Grey symbols, mutual part of both 
gradients. Normal triangle, reverse triangle, square, diamond, and rotundity, 


fish assemblages from first- to fifth-order streams, respectively. 


SIMPER analysis revealed that Z. platypus was 
more abundant in first-, second- and third-order streams 
and less in fourth- and fifth-order streams within gradient 
B. The occurrence and abundance of C. rarus, P. fasciata, 
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R. ocellatus, and S. nigripinnis determined how fish 
assemblages varied substantially from first- to fourth- 
order streams. Ctenogobius sp. and C. argus contri- 


buted most to the divergence in fish assemblages betw- 
een fourth- and fifth-order streams within Gradient A 
(Table 5). 


Table 5 Diagnostic species by SIMPER analysis for first- to fifth-order streams along gradients A and B within Qingyi Streams 











Gradient A Gradient B 
AvAbu. AvSim. Contri (%) AvAbu. AvSim. Contri (Yo) 

1st-order AvSim., 46.3% AvSim., 44.1% 

Zacco platypus 60.9 18.5 40.0 167.3 9.4 30.2 

Acrossocheilus fasciatus 39.6 10.2 22.1 17.3 1.5 17.9 

Vanmanenia stenosoma 16.3 5.1 11.0 

Cobitis rarus 23.5 24 19.3 
2nd-order AvSim., 51.6% AvSim., 52.9% 

Zacco platypus 63.1 16.8 32:5. 37.0 741 33.5 

Acrossocheilus fasciatus 15.3 9.7 18.9 15.0 3.1 24.7 

Ctenogobius sp. 14.4 7.4 14.4 

Parabotia fasciata 16.0 7.6 14.4 
3rd-order AvSim., 48.6% AvSim., 57.8% 

Zacco platypus 44.8 11.9 24.6 95.8 15.9 27.6 

Rhodeus ocellatus 22.1 9.9 20.5 69.7 11.4 19.7 

Ctenogobius sp. 28.5 6.6 13.6 28.8 8.4 14.6 
4th-order AvSim., 51.0% AvSim., 64.3% 

Ctenogobius sp. 24.6 11.9 23.4 

Zacco platypus 46.1 8.5 16.8 

Acrossocheilus fasciatus 12.6 6.1 11.9 51.8 22 8.1 

Sarcocheilichthys nigripinnis 36.5 5.7 8.8 

Abbottina rivularis 20.3 4.8 7.5 
5th-order AvSim., 52.9% AvSim., 52.9% 

Pseudobagrus albomarginatus 58.0 7.0 15.1 58.0 7.0 15.1 

Channa argus 16.8 4.5 92 16.8 4.5 9.2 

Carassius auratus 15.8 3.4 7.3 15.8 3.4 7.3 


The first three species contributing most to the average similarity within each order are shown. AvSim., average similarity (%); AvAbu., average abundance; 


Contri., contribution. 


Relationships between species and environment 

The first and second axes of the CCA ordination 
accounted for 49.7% and 12.496 of total variance in 
species richness and abundance among sampling sites, 
and for 54.7% and 13.6% of variance in species data and 
environmental relationship, respectively. From CCA, 
gradient A separated fish assemblages based on elevation, 
substrate, and flow velocity on the right, and wet width, 
water depth, discharge, and conductivity on the left. The 
species related to this gradient were Е. ilishaeformis, P. 
pekinensis, S. nigripinnis, H. labeo, C. rarus, C. argus, S. 
chuatsi, P. parva, C. auratus, A. chankaensis, R. 
ocellatus, C. sinensis, Ctenogobius sp., A. fasciatus, and 
Z. platypus (left to right). Gradient B was caused by 
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water temperature and dissolved oxygen on the top and 
pH on the bottom, along which P oxycephalus, P. 
fasciata, A. chankaensis, P. vaillanti, M. albus, G 
taeniellus, and S. argentatus were distributed from top to 
bottom (Figure 4). 


DISCUSSION 


Our results showed that along a longitudinal 
gradient without the direct effect of a large dam, fish 
species richness was highest in fourth-order streams and 
showed no significant difference among first- to third- 
order streams. Following the increase in stream size and 
habitat diversity from headwaters to downstream, fish 
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Figure 4 Canonical correspondence analysis (CCA) diagrams 
for environmental variables and fish assemblages in Qingyi 
Stream 
Black arrows indicate the most important factors and grey arrows show the 


other factors. Species codes as in Table 1. 


species richness generally increases downstream (Gorm- 
an & Karr, 1978; Matthews, 1986). Based on fish speci- 
mens collected from 89 shallow (first- to third-order) 
tributaries within the Qingyi watershed, Yan et al. (2011) 
discovered that local species richness was determined by 
local habitat variables (e.g., wetted width and water 
temperature) but not tributary spatial variables (e.g., 
stream order and link magnitude). Yan et al (2010) also 
revealed that stream order was not the optimal frame- 
work explaining spatial variation in fish assemblages 
within Puxi Stream, a tributary above Chencun Reservoir 
in the Qingyi watershed. The results in our study that fish 
species richness showed no substantial variation across 
first- to third-order streams appears to support fish 
species responding to local habitat features but not to 
stream classification schemes of humans (Matthews, 
1998). In addition, our findings that the highest species 
richness occurred in fourth-order, not fifth-order, streams 
is consistent with the river continuum concept (RCC) 
claiming that maximum species diversity often occurs in 
mid-sized, not large, streams. Oberdorff et al (1993) also 
found a decline in fish species richness in the lower 
reaches of the Seine River, which was explained by 
anthropogenic disturbances decreasing habitat diversity. 
In addition to species richness, fish also exhibit an 
alteration in zonation along the upstream-downstream 
gradient (Lasne et al, 2007; Miranda & Raborn, 2000). 
One of the most recognized theoretical classifications 
dividing running waters based on fish is longitudinal fish 
zonation proposed by Huet (1959); the brown trout, 
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grayling, barbel and bream zones exhibit alteration, in 
turn, from headwaters to lower reaches. However, asso- 
ciated with continuous, not abrupt, variation in the 
natural environments of many streams, fish assemblages 
vary by “addition” rather than by “replacement” (Matthe- 
ws, 1998), which suggests a deficiency in the application 
of this zonation concept within these streams. In our 
results, despite some species only being collected in one 
ordered stream (e.g., P. oxycephalus and С. molitorella in 
first-order streams, C. carpio, M. piceus, and H. 
leucisculus in fifth-order streams), most fish species 
were distributed in at least two different ordered streams. 
Based on our SIMPER analysis, along the longitudinal 
gradient without the effects of a large dam, changes in 
the relative abundances of Z. platypus, C. rarus, P. 
fasciata, R. ocellatus, and S. nigripinnis contributed most 
to longitudinal species replacement. Our results suggest 
that “addition”, not “replacement”, was the main underl- 
ying mechanism explaining species distributions from 
upstream to downstream in Qingyi Stream. This spatial 
variation in species distributions could be explained by 
the relationship between species and the local habitat 
environment based on CCA. The spatial distributions of 
fish species in this study were affected by a series of 
local habitat environmental factors, such as elevation, 
substrate, water depth, current discharge, wetted width, 
and conductivity. At the watershed scale, Yan et al (2011) 
reported that fish assemblage structure in this study area 
was related to the combined effects of local habitat (i.e., 
elevation, substrate, and water depth) and spatial stream 
position (i.e. magnitude link and confluence link), 
suggesting that fish species distribution was associated 
with their ecological requirements, such as habitat prefe- 
rence and trophic ecology. 

In Qingyi Stream, fish assemblages presented diffe- 
rent patterns in their spatial variations between the two 
longitudinal gradients, which could be viewed as the 
effect of Chencun Hydropower Station. Fish showed 
significantly lower species richness and different species 
composition within the gradient disrupted by the dam 
compared with the gradient without the large dam. 
However, the substantial variations in fish assemblages 
between gradients were only observed in fourth-order 
streams. Numerous studies have revealed that dams may 
impact stream fish by blocking migratory pathways 
(March et al, 2003), decreasing current velocity upstream 
(Bennett et al, 2002), altering natural flow regime 
(Bonner & Wilde, 2000), and decreasing water temp- 
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erature downstream (Clarkson & Childs, 2000). In this 
study, when considering first- to fourth-order streams 
within the two gradients independently, a total of 18 fish 
species (e.g., D. tumirostris, E. bambusa, S. hollandi) 
were collected in gradient A only, and eight species (e.g., 
O. bidens, H. bleekeri, S. parvus) were more frequent 
and abundant in gradient A than in gradient B. This may 
be explained by two causes resulting from the effects of 
large dams on fish assemblages. Firstly, most of the 18 
species only occurring in gradient A, including Н. 
bleekeri, E. ilishaeformis, H. maculatus, and S. roulei, 
inhabited the lower reaches of the streams. The Chencun 
Hydropower Station may play a role in constraining their 
upstream movements, resulting in the local extinction of 
these fish in areas above the station. Secondly, dams may 
habitat | conditions 


modify local upstream апа 
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Population sizes and group characteristics of Siberian Crane (Leuco- 
geranus leucogeranus) and Hooded Crane (Grus monacha) in Poyang 
Lake Wetland 


Ming-Qin SHAO', Hong GUO, Jian-Hong JIANG 


College of Life Science, Jiangxi Normal University, Jiangxi Provincial Key Lab of Protection and Utilization of Subtropical Plant Resources, Nanchang Jiangxi 
330022, China 


Abstract: Both the Siberian Crane (Leucogeranus leucogeranus) and Hooded Crane (Grus monacha) have limited population sizes 
and are considered endangered by domestic Chinese and international agencies. To document the current size of their respective 
populations and characterize their groups, between October 2012 and April 2013 we undertook fieldwork at four nature reserve areas 
within the Poyang Lake wetlands. We divided Poyanghu National Nature Reserve (PYH) into the Wucheng (PWC) and Hengfeng 
areas (PHF), because each are each located in different counties. Our fieldwork showed that the Siberian Crane occurred mainly in 
PYH (364 in the PHF, 158 in the PWC) and the Nanjishan Wetland National Nature Reserve (NJS, with 200 individuals). The 
Hooded Crane was mainly distributed in PYH (302 in the PHF and 154 in the PWC). Family groups accounted for more than 50% of 
the total number of groups among both species, with Hooded Cranes forming more family groups than Siberian Cranes. Typically, 
these groups were formed of two adults with one offspring (Siberian Crane), and two adults with two offspring (Hooded Crane), with 








the mean family group size of the Siberian Crane and Hooded Crane being respectively 2.65+0.53 (n=43) апа 3.09+0.86 (n=47) 





individuals per group. The mean collective group size of the Siberian Crane and Hooded Crane included 28.09+24.94 (n=23) and 





28.94+27.97 (n=16) individuals per group, respectively, with the proportion of juveniles among Hooded Cranes being more than 


double that seen among the Siberian Cranes. 


Keywords: Siberian Crane; Hooded Crane; Group characteristics; Poyang Lake wetland 


The Siberian Crane (Leucogeranus leucogeranus) 
and Hooded Crane (Grus monacha) are frequent winter 
visitors of the Poyang Lake wetland area, but both 
species have suffered dramatic losses. Currently, the 
Siberian Crane is listed in the first category of the 
nationally protected wildlife species in China and 
regarded as a critically endangered species by IUCN, 
with a global population size ranging from 3,800 to 
4,000 individuals (BirdLife International, 2008; Li et al, 
2012). Though the Siberian Crane breeds in northwest 
Siberia, it frequently winters around Poyang Lake in 
Jiangxi Province and Dongting Lake in Hunan Province, 
with more than 90% of the known global population 
wintering around Poyang Lake (Zheng & Wang, 1998; 
Shan et al, 2012). The Hooded Crane has a somewhat 
larger global population—around 11 600 individuals— 
only ~1000 winter in China (BirdLife International 2008; 
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Zhou et al, 2010). The Hooded Crane is also listed in the 
first category of the nationally protected wildlife species 
in China, but is regarded as a vulnerable species by 
IUCN. Typically, the Hooded Crane breeds in eastern 
Siberia and Heilongjiang Province in China, but winters 
in an array of areas in eastern Asia, including Korea, 
Japan and China, predominately along lakes beside the 
middle and lower reaches of Yangtze River, including 
Poyang Lake in Jiangxi Province, Dongting Lake in 
Hunan Province, Shengjin Lake and Longgan Lake in 
Anhui Province, and Chongming Island in Shanghai 
(Jing et al, 2002b). 
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To date, research on the wintering ecology of 
Siberian Crane and Hooded Crane has mainly focused on 
numbers and distribution (Liu, 2001; Li et al, 2012; Shan 
et al, 2012; Yan & Ding, 1988; Zhou & Ding, 1982;), 
genetic structure (Zhang et al, 2012), behavior (Jing et al, 
2002b; Li et al, 2007; Zhou et al, 2010), diet (Zhao et al, 
2002) and habitat selection (Jing et al, 2002a; Sun & 
Huang, 2010; Zhang et al, 2010). However, despite the 
popularity and importance of the Poyang Lake are in the 
north of Jiangxi Province (and the largest freshwater lake 
in China), no comprehensive and quantitative assessment 
on population dynamic or age structure of these two 
species has been conducted. In the present study, we 
report the findings of several field observations of this 
area undertaken between October 2012 and April 2013, 
where we group 
characteristics and age structure. To facilitate more targ- 


documented population size, 


eted conservation and management of these two crane 
species, we focused on both documenting population size 
dynamics of the two crane species at different lakes in 
the area, and comparing the group characteristics and age 
structure of the two species. 


STUDY AREAS AND METHODS 


Poyang Lake (E115?47'-116?45', М28922-29945” 
is a key wintering or stopover sites for migratory 
waterbirds in the East Asian-Australasian Flyway. 
Poyang Lake had a humid subtropical climate with the 
annual temperature averaging 16.7-17.7 °С, and average 
annual precipitation at 1400—1900 mm (Huang & Guo, 
2007). Between October and April, the Poyang Lake 
undergoes a drought period, during which many 
habitats—e.g., grassland and mudflats—form and attract 
-500000-600000 waterbirds annually. During the entire 
year, the wetland vegetation at the Poyang Lake is 
dominated by Carex spp. Phragmites australis and 
Triarrhena lutarioriparia (Liu et al, 2006). For the 
present study, we selected four nature reserves in Poyang 
Lake to survey the population sizes and group 
characteristics of the Siberian and Hooded Cranes. To 
facilitate more accurate measurement, the Poyanghu 
National Nature Reserve (PYH) was divided into two 
areas, the Wucheng (PWC) and Hengfeng areas (PHF), 
since they are located in two counties and were 
investigated at different times. In the PWC, four lakes 
(Dahuchi, Baziqiang, Changhuchi, Zhushihu) were selec- 
ted as observation sites, with additional surveys at two 
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lakes (Shahu, Banghu) in the PHF. We also surveyed 11 
lakes at the Nanjishan Wetland National Nature Reserve 
(NJS; Zhanbeihu, Sanhu, Changhu, Fengweihu, Sanni- 
wan, Baishahu, Linghu, Shangbeijia, Xiabeijia, Beishen- 
hu, Shentanghu). A further 10 lakes in the Duchang 
Migratory Bird Nature Reserve (DC; Jishanhu, Zhongba, 
Xinmiaohu, Mayinghu, Huangjinzui, Binhu, Shuhu, 
Huamiaohu, Henggang, Xihu) and eight lakes from the 
Poyang Baishazhou Nature Reserve (BSZ; Zhuhu, Che- 
men, Rongqi, Sishilijie, Stwanghu, Xiaominghu, Dam- 
inghu, Biaoen) were also surveyed (Figure 1). 
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Figure 1 Surveyed areas around Poyang Lake 

Abbreviations are as follows: Dahuchi: DHC, Baziqiang: BZQ, Changhuchi: 
CHC, Zhushihu: ZSH, Shahu: ShaH, Banghu: BH, Zhanbeihu: ZBH, Sanhu: 
SH, Changhu: CH, Fengweihu: FWH, Sanniwan: SNW, Baishahu: BaiSH, 
Linghu: LH, Shangbeijia: SBJ, Xiabeijia: KBJ, Beishenhu: BeiSH, 
Shentanghu: STH, Jishanhu: JSH, Zhongba: ZB, Xinmiaohu: XinMH, 
Mayinghu: MYH, Huangjinzui: HJZ, Binhu: BH, Shuhu: SH, Huamiaohu: 
HMH, Henggang: HG, Xihu: XH, Zhuhu: ZH, Chemen: CM, Ronggi: RQ, 
Sishilijie: SSLJ, Siwanghu: SWH, Xiaominghu: XMH, Daminghu: DMH, 
Biaoen: BE. 


Between October 2012 and April 2013, six surveys, 
each lasting one or two days, were conducted from at 
each reserve at intervals of 20—30 days. For comparison, 
we divided the wintering season into three stages based 
on the population dynamics: early stage (October to 
December); middle stage (January to February); and late 
stage (March to April), in line with a previous study 
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(Zhou et al, 2010). Totally, between one and five fixed 
observation points were selected at each lake to count the 
number of adult or juvenile cranes, group size, group 
types of two crane species with binoculars (8x) and 
spotting scope (20—60x). The fixed points varied some- 
what according to water level. We noted that the group 
types comprised of family, collective and special group 
(Liu et al, 2008), with the family group being sub- 
divided into four types: 1) two adults; 2) two adults with 
one offspring; 3) two adults with two offspring; and 
4) one adult with one offspring. The collective group 
referred to groups including five or more individuals, 
while the special group referred to a single adult or 
juvenile crane. Groups were defined as cranes that were 
usually within less than 50 m of each other and appeared 
to coordinate their activities, 1.e., they stayed near one 
another and moved in the same direction (Lcheureux et 
al, 1995; Liu et al, 2009; Liu et al, 2010). 

Our observation data were first checked with the 
Kolmogorov-Smirnov Test and were found to fit a 
normal distribution. Next, one-way ANOVA was used to 
test for differences in the average size of the family or 
collective groups at different stages, and the size of the 
family or collective groups between Siberian Crane and 
Hooded Crane. All data were analyzed with SPSS 19.0 
(SPSS inc., Chicago, IL, USA). 


RESULTS 


Population size and distribution 


The Siberian Crane was recorded at 18 of the 35 
surveyed lakes located in the four nature reserve areas, 
but were mainly distributed at PYH and NJS. A large and 
steady population of Siberian Cranes was observed in 


Shahu, Banghu, Dahuchi, Bazigiang within the PYH, and 
30 O Family Ш Specia © Collective 


25 


ы 
e 


Percentage (%) 
о 

















10 
5 
0 MEE iue: П 
Early Middle 
Time 





Percentage (%) 


375 


we also observed a large population in Sanhu, Baishahu, 
Shangbeijia, Beishenghu, Shentanghu at NJS. The 
Siberian Crane was not observed regularly at any of the 
other surveyed sites in the other reserves (Table 1). 

The Hooded Crane was recorded at 10 lakes located 
within three of the nature reserve areas, but not in the DC 
(Table 1). The Hooded Crane was mainly distributed in 
Shahu, Banghu, Dahuchi, Changhuchi and Zhushihu of 
PYH, with a large and steady population of the Hooded 
Crane also observed in Banghu and Zhushihu. Though 
we also observed the Hooded Crane at the BSZ, the 
population was comparatively small. 

The number of the Siberian Crane peaked in late 
December and early January when they were mainly 
distributed in Shahu (333 individuals), Baziqiang (151 
individuals) and Beishenghu (182 individuals) but then 
decreased in mid-March (Table 1). Conversely, the 
number of Hooded Cranes was high in late January and 
peaked in later February and early March, when the 
cranes were mainly distributed in Banghu (226 individ- 
uals) and Zhushihu (135 individuals). 


Group characteristics 

During the course of our field observations, a total 
of 56 groups of Siberian Crane were observed, with 
family and collective groups being the most frequently 
encountered, accounting respectively for 55.54% and 
41.07% of the total Siberian Crane groups observed 
during this investigation. The family groups of the 
Siberian Crane were the most frequently observed during 
all three stages. For the Hooded Cranes, 65 groups were 
observed during this study, with family and collective 
groups respectively accounting for 72.31% and 24.62% 
of the total groups. Similarly, the family groups of the 
Hooded Crane were also the most often observed groups 
during all three stages (Figure 2). 
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Figure 2 Temporal variation in group type among Siberian Cranes (left) and Hooded Cranes (right) around Poyang Lake 
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Table 1 Spatial and temporal distribution of Siberian Crane and Hooded Crane around Poyang Lake 




















Study area Site 10/15* 11/19 1/2 1/27 3/7 4/13 
Shahu 0 21 333 15 0 0 
шыла Banghu 0 23 31 10 0 51 
rane 
total 0 44 364 25 0 51 
PHF 
Shahu 0 0 2 0 76 0 
Hooded Crane Banghu 0 28 0 113 226 0 
total 0 25 2 113 302 0 
10/14 11/12 12/24 1/16 2/25 3/8 
Dahuchi 0 3 7 11 22 25 
Siberian С; 
Стапе Baziqiang 0 15 151 0 3 9 
total 0 18 158 11 25 34 
Dahuchi 0 2 9 33 15 4 
PWC 
Baziqiang 0 0 4 0 4 0 
Hooded Crane — Chahuchi 0 57 0 0 0 4 
Zhushihu 0 15 15 27 135 0 
total 0 74 28 60 154 8 
10/19 11/11 12/23 1/20 2/22 3/9 
Sanhu 0 0 3 30 0 0 
Fengweihu 0 0 0 0 3 0 
Baishahu 0 0 0 35 0 0 
Siberian Linghu 0 0 0 2 0 0 
Crane Shangbeijia 0 0 15 48 39 20 
Beishenghu 0 0 182 14 0 4 
NJS 
Shentanghu = = = 39 0 0 
total 0 0 200 168 42 24 
Changhu 0 0 0 9 4 0 
Shangbeijia 0 0 0 4 0 0 
Hooded Crane 
Xiabeijia 0 0 0 0 0 4 
total 0 0 0 13 4 4 
11/4 12/2 12/29 1/23 3/14 4/10 
Binhu 0 0 8 0 0 0 
Siberian Xiaoyang 0 0 5 0 0 0 
DC 
Crane Chian 0 0 2 0 0 0 
total 0 0 15 0 0 0 
10/28 11/23 1/9 1/25 3/2 4/2 
Chemen 0 5 3 0 0 0 
Rongqi 0 0 3 3 6 0 
pd Xiaominghu 0 0 6 0 0 0 
BSZ rane 
Daminghu 0 0 3 7 8 0 
total 0 5 15 10 14 0 
Hooded Crane Віаоеп 0 0 0 4 0 
Siberian Crane 0 67 752 214 81 109 
Total 
Hooded Crane 0 99 30 186 464 12 


* date, e.g., 10/15 = 15 October; — not surveyed; РНЕ: Poyang National Nature Reserve (Hengfeng); PWC: Poyang National Nature Reserve (Wucheng) ; NJS: 
Nanjishan Wetland National Nature Reserve; DC: Duchang Migratory bird Nature Reserve; BSZ: Poyang Baishazhou Nature Reserve 
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While family groups were the most common 
grouping for both the Siberian and Hooded Cranes, there 
composition differed between species. For the Siberian 
Crane, of four types of family groups, those with two 
adults and one offspring (51.61%) were most frequently 
observed, followed by family groups of two adults 
(38.71%). For the Hooded Cranes, family groups were 
typically two adults with two offspring (40.43%), though 
groups with two adults and one offspring (27.66%) and 
two adults (31.91%) were also regularly observed. No 
family group with one adult and one offspring was 
observed for the Hooded Crane at any of our survey sites 
(Figure 3). 

The mean family group size of Siberian Crane and 
Hooded Crane were respectively 2.65+0.53 (n=43) and 
3.09--0.86 (n=47) individuals per group (Figure 4). The 
mean collective group size of the Siberian Crane and 
Hooded Crane were respectively 28.09+24.94 (n=23) 
and 28.94--27.97 (n=16) individuals per group. The mean 
family group size of the Hooded Crane was significant 
larger than that of the Siberian Crane (F=8.191, а/-88, 
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P<0.05). The collective group size between the two 
species of crane showed no significant difference 
(F=0.01, df=37, P>0.05). 

During our survey, we observed no significant 
difference in the mean family group size for the Siberian 
Crane (F=2.023, df-40, P>0.05) throughout the three 
stages, but the mean family group size of the Hooded 
Crane was significantly lower in the early as compared to 
the late stage (P«0.001). Surveys of both species showed 
no significant differences in their mean collective group 
size between the three stages (Siberian Crane: F=0.396, 
df=20, P>0.05; Hooded Crane: F=0.834, df=13, P>0.05). 


Age structure 

Totally, 737 Siberian Cranes and 243 Hooded 
Cranes could be categorized as either being adults or 
juveniles. The percentage of juveniles among Siberian 
Cranes was 12.48%, while among Hooded Cranes was 
29.22%. The numbers of juveniles displayed the same 
pattern among both species, being low in early or mid- 
winter, and high in late winter (Table 2). 
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Figure 3 Statistics of four family group types of Siberian Cranes (left) and Hooded Cranes (right) at three different stages 
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Figure 4 Temporal variation of family group size (left) and collective group size (right) of two Crane species in Poyang Lake 
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Table 2 Temporal variations in age structure of two species cranes in Poyang Lake 








Species Early stage Middle stage Late stage Total 
"T Proportion of juveniles 9.89% 6.92% 42.86% 12.48% 
Siberian Crane e 
No. of observed individuals 364 289 84 737 
Proportion of juveniles 23.73% 27.61% 40.00% 29.22% 
Hooded Crane Који 
No. of Observed individuals 59 134 50 243 
breed, and may have delayed their migration after the 
> y y g 
DISCUSSION had acquired enough energy at Poyang Lake, which may 


Population sizes and distribution 

In 1980, Zhou & Ding (1982) recorded the first 
sighting of more than 100 Siberian Cranes at Poyang 
Lake. Over the next six years, their numbers have incre- 
ased gradually, peaking at 1 600+ individuals in 1987 
(Yan & Ding, 1988). Between 1998 and 2010, the num- 
ber of the Siberian Crane at Poyang Lake fluctuated from 
1627 to 3954 individuals (Shan et al, 2012). Recent rese- 
arch indicates that the global population size somewhere 
between 3 800-4000 individuals (Li et al, 2012). This 
increase suggests that Poyang Lake is an ideal wintering 
habitat for the Siberian Crane. During our study, we 
noted that these cranes are primarily distributed across 
two national nature reserves, with maximum numbers in 
PHF, PWC and NJS being 364, 158 and 200 individuals, 
accounting for 9.10%, 3.95% and 5.00% of the global 
population size, respectively. This shows that these two 
reserves are important habitats for the Siberian Crane. Com- 
pared with previous surveys of the Siberian Cranes, our 
results showed a smaller and more scattered distribution. 
For example, previous studies indicated that the Siberian 
Crane was mainly distributed in Banghu and seldom 
observed in the NJS (Shan et al, 2012), but we observed a 
large number of these cranes in the NJS, and detected 
frequent local movements of the species among different 
lakes in NJS. 

In recent years, continuous low water levels during 
spring and early summer in Poyang Lake have resulted in 
long-term drought and a decreased availability of food 
for many of the water birds that migrate there. The 
scattered distribution pattern of the Siberian Crane we 
observed during our surveys may accordingly be due to 
an insufficient food supply capable of support all the 
cranes in Banghu. Similarly, the departure time from 
Poyang Lake recorded for the Siberian Crane during this 


study was later than the time recorded in previous reports. 


One potential explanation for this difference may due to 
a delay in the wintering period due to climate change. 
Additionally, the juveniles we observed did not need to 
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subsequently have decreased their mortality. 

Our observations also showed that the Hooded 
Crane was mainly distributed in the PYH, with maxi- 
mum numbers in PHF and PWC respectively at 302 and 
154 individuals, indicating a larger number than previ- 
ously observed in the entire Poyang Lake (Li et al, 
2012).While the Siberian Crane was mainly distributed 
in Shahu, Banghu, Dahuchi, Beishenghu, the Hooded Crane 
was only found in Banghu, Dahuchi and Zhushihu. 
Similarly, the Siberian Crane mainly selected shallow water 
areas or mudflats for forging habitats, while the Hooded 
Crane foraged mainly in the grasslands near the lakes. The 
migration of the Siberian Crane also took place somewhat 
earlier than that of the Hooded Crane. These ecological 
differences demonstrate the temporal and spatial separation 
indicate a positive coexistence of the two species. 


Group characteristics 

On the whole, both Siberian and Hooded Crane 
groupings were dominated by family and collective groups. 
The proportion of family groups for the Hooded Crane 
recorded in this study was higher than that observed in 
Chongming Island (Jing et al, 2002b). Family groups for 
the Siberian Crane were dominated by two adults with one 
offspring and two adults with two offspring for the 
Hooded Crane. Compared with the Black-necked Crane, 
the Siberian Crane had a higher proportion of two adults 
with one offspring and a lower proportion of two adults 
with two offspring. The clutch size of the Hooded Crane 
was two (n=4) with a hatching rate of 100% (Guo et al, 
2005). A previous study in 1981 noted that family groups 
of the Hooded Crane were dominated by two adults with 
one offspring in Poyang Lake, but by two adults with two 
offspring in Japan in 1967 (Wang, 1988). The higher 
proportion of two adults with two offspring for the 
Hooded Crane that we noted during our observations 
indicates a high survival rate of juveniles. This assessment 
is also bolstered by the observed higher proportion of two 
adults with two offspring for the Hooded Crane at Poyang 
Lake than previously observed in Shengjinhu, or among 
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the Black-necked Crane (Liu, 2001; Liu et al, 2008). 


Age structure 

Our characterization of the age structure of both 
crane species showed some interesting differences 
between Hooded Cranes and Siberian Cranes. The 
proportion of juveniles of the Hooded Crane were more 
than double those seen among Siberian Cranes, 
suggesting that the proportion of juveniles of the 
Siberian Crane recorded in this study was similar to the 
numbers recorded in 1981, and also in line with 
observations of Black-necked Cranes (Zhou & Ding, 
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Abstract: The Asian swamp eel (Monopterus albus) is one of the most economically important freshwater fish in East Asia, but data 
on the immune genes of M. albus are scarce compared to other commercially important fish. A better understanding of the eel’s 
immune responses may help in developing strategies for disease management, potentially improving yields and mitigating losses. In 
mammals, interferon regulatory factors (IRFs) play a vital role in both the innate and adaptive immune system; though among 
teleosts /RF4 and IRF10 have seldom been studied. In this study, we characterized JRF4 and IRF10 from M. albus (maIRF4 and 
malRF10) and found that maIRF4 cDNA consists of 1716 nucleotides encoding a 451 amino acid (aa) protein, while maIRF10 
consists of 1 744 nucleotides including an open reading frame (ORF) of 1236 nt encoding 411 aa. The maIRFI0 gene was 
constitutively expressed at high levels in a variety of tissues, while maIRF4 showed a very limited expression pattern. Expression of 
malRF4 and та!КЕ10 in head kidney, and spleen tissues was significantly up-regulated from 12 h to 48 h post-stimulation with 
polyinosinic: polycytidylic acid (poly I:C), lipopolysaccharide (LPS) and a common pathogenic bacteria Aeromonas hydrophila. 
These results suggest that IRF4 and IRF10 play roles in immune responses to both viral and bacterial infections in M. albus. 


Keywords: Monopterus albus; IRF4; IRF10; Poly I:C; Up-regulation 


Interferon regulatory factors (IRFs) are ancient 
molecules conserved throughout the evolution of metaz- 
oans and play a vital role in the innate and adaptive 
immune system (Tamura et al, 2008). To date, 11 IRF 
family members (IRF1-11) have been described in 
vertebrates and invertebrates (Huang et al, 2010), with 
IRF1-10 being present in most vertebrate species and 
IRF11 being found in non-vertebrate deuterostomes 
(Huang et al, 2010). Phylogenetic analysis of these 11 
IRF proteins demonstrated that they can be subdivided 
into four groups that reflect their evolutionary history 
(Nehyba et al, 2009; Xu et al, 2010). Previous studies 
showed that in humans (Homo sapiens) and mice (Mus 
musculus), IRF4 was expressed in most types of immune 
cells, and has critical functions in B cell differentiation 
and immunoglobulin production (De et al, 2012). 
Likewise, IRF4 seems to also play important roles in the 
development and function of T helper cells, regulatory T 
(Treg) cells, dendritic cells (Xu et al, 2012) and Ср4 T 
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cell differentiation (Suzuki et al, 2004). Other studies on 
chickens (Gallus gallus) found that IRF4 was mainly 
expressed in the bursa, bursal lymphocytes (Nehyba et al, 
2002), and thymus (Dougherty et al, 2009) and it is 
capable of repressing the expression of ovalbumin gene 
(Dougherty et al, 2009). ConA can induce the expression 
of IRF4 in splenic cells, while IFNs cannot induce the 
expression of IRF4 (Nehyba et al, 2002). Furthermore, 
among teleost, pathogen-associated molecular patterns 
(PAMPs) were found to stimulate the IRF4 expression in 
rainbow trout (Oncorhynchus mykiss) (Holland et al, 
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2010) and rock bream (Oplegnathus fasciatus) (Bathige 
et al, 2012). 

Belonging to the same IRF4 subfamily, IRF10 was 
previously found to have been eliminated or rendered 
non-functional in both mice and humans (Nehyba et al, 
2009). Curiously though, IRF10 can elevate the expre- 
ssion of major histocompatibility complex (MHC) class I 
molecules and guanylate-binding protein (GBP) and 
interfere with the induction of the type I IFN target genes 
in chickens (Nehyba et al, 2002). Both MHC class I and 
GBP play important roles in viral infections (Nehyba et 
al, 2002; Hu et al, 2011). Another feature of IRF10 is that 
similar to IRF4, it can also repress the expression of 
ovalbumin gene (Dougherty et al, 2009) and ConA can 
induce the expression of IRF10 in splenic cells (Nehyba 
et al, 2002). Furthermore, type I IFN and IFN-y seem 
capable of inducing the expression of IRF10 in primary 
fibroblasts, though the expression of IRF10 is induced 
relatively late and needs other protein synthesis (Nehyba 
et al, 2002). In the Japanese flounder (Paralichthys 
olivaceus), bacteria or viral hemorrhagic septicemia virus 
(VHSV) can increase the expression of IRF10 in kidney 
tissue (Suzuki et al, 2011), though in zebrafish (Danio 
rerio) IRF10 was found to be a negative regulator of 
interferon transcription (Li et al, 2013). Clearly the 
diversity of roles and effects of IRF10 could use some 
clarification, especially regarding the role IRF10 plays 
in immune responses following viral or bacterial 
infections. 

In recent years, the Asian swamp eel (Monopterus 
albus) has become one of the most economically impo- 
rtant freshwater fish in East Asia. Unfortunately, data on 
immune genes of М. albus are scarce as compared to 
other commercially important fish. Given the eels econ- 
omic importance, a clearer understanding of its immune 
responses may help to develop strategies for disease 
management, which may potentially aid in aquaculture, 
increase yields, or decreases loses of this species. In this 
study, we cloned full-length cDNAs of М. albus IRF4 
(maIRF4) and maIRF10 and then investigated the tissue 
distribution of these two genes’ expression. Paired with 
this analysis, we also opted to investigate a common 
pathogen afflicting the Asian swamp eel, Aeromonas 
hydrophila, the causative agent of Septicemia in this 
species (Yang et al, 2008; He et al, 2010). Pathogens of 
some disease such as stigmatosis in Asian swamp eel had 
not been investigated clearly. For this analysis, we 
measured the change in maIRF4 and maIRF10 gene 
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expression in the head kidney (HK), spleen (SP), and 
gills (GI) following stimulation with PAMPs, e.g. 
polyinosinic: polycytidylic acid (poly I:C) and lipopol- 
ysaccharide (LPS), and a common pathogenic bacteria, 
A. hydrophila. 


MATERIAL AND METHODS 
Fish 

Eels of 150—200 р body weight were obtained 
from the aquaculture base of Yangtze University, China. 
Eels were maintained in fiberglass tanks supplied with a 





continuous flow of recirculating freshwater at 24-1 °C 
Fish were fed once daily on water earthworms, and 
were acclimated for at least two weeks prior to experi- 
mentation. 


Gene cloning and sequence analysis 

cDNA samples were prepared from head kidney 
with the first strand cDNA synthesis kit (Fermentas, 
Canada). First, degenerate primers IRF4-F/IRF4-R and 
IRF10-F/IRF10-R were designed against conserved regi- 
on and used for PCR amplification of the internal region 
of maIRF4 and maIRF10 that was obtained from cDNA 
samples. PCR products were respectively isolated using 
a Gel Extraction Kit (Tiangen, China), cloned into a 
pMD18-T vector (TaKaRa, Japan) and transformed into 
Escherichia coli strain DH5a competent cells. Putative 
clones were then screened via PCR using the aforem- 
entioned primers, and the selected clones were seque- 
nced. 5'-RACE and 3'-RACE were performed with gene- 
specific primers and adaptor primers to obtain the full- 
length cDNA sequence of maIRF4 and maIRF10. Unive- 
rsal primers mix (UPM) was obtained from mixing UPM 
Long and UPM Short at a ratio of 1:100. For 3'-RACE, 
PCR was initially performed with primers UPM/3-F1 
followed by a nested PCR with primers UPM/3-F2. For 
5'-RACE, the adaptor primers AAP and AUAP were used. 
RNA from M. albus spleen was reverse-transcribed at 42 
°С using the gene-specific primer 5-К1. Following 
synthesis of the first strand cDNA, the resulting purified 
cDNA was used in the TdT-tailing reaction, and the 
talled cDNA was then amplified by primers 5-R2 and 
AAP. A dilution of the original PCR (0.196) was re- 
amplified using AUAP and a nested 5-R3 primer (all 
primers listed in Table 1). 

The nucleotide sequences generated were 
assembled and analyzed with AlignIR (LI-COR, Inc.). 
Sequence identities were calculated using the DNAStar 
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Table 1 Oligonucleotide primers used in the study 





Name Sequence (5'—3”) Usage 
UPM Long CTAATACGACTCACTATAGGGCAAGCAGTGGTATCAACGCAGAGT 3'-Race PCR universal primers 
UPM Short CTAATACGACTCACTATAGGGC P 
AAP GGCCACGCGTCGACTAGTACGGGIIGGGIIGGGIIG 5'-Race PCR universal primers 
AUAP GGCCACGCGTCGACTAGTAC P 
IRF4-F GA(G)AGC(T)CAGCTGGACATCTC Я { , 
IRF4-R ACAGGAGCTGCCTGGCA(G)AAC Cloning for maIRF4 internal fragment 
IRF4-5-R1 CTCCCGCATCACAG 5’ RACE 1“ round PCR 
IRF4-5-R2 GAGAAAGGTGAAGGCTGGCTGT 5' RACE 2" round PCR 
IRF4-5-R3 GGGATGATGCGGTAGACTTTGT 5' RACE 3" round PCR 
IRF4-3-F 1 ACGGAGACAAGCCCAACAAGC 3’ RACE 1% round PCR 
IRF4-3-F2 GTCGGCTGCTGCCTCGTTTCC 3’ RACE 2" round PCR 
IRF10-F GAGCGC(G)AA(G)CCAGCTGGACATC : 2 
IRFIO-R CTCC(G)CGCTCCAGC(T)TTGTTGGG Cloning for maIRF10 internal fragment 
IRF10-5-R1 CTGGGCTCTTGGTG 5’ RACE 1“ round PCR 
IRF 10-5-R2 GCTCTTGGTGGTCACTTTCATT 5' RACE 2" round PCR 
IRF10-5-R3 AAGCGGGCTGAAGAAGGTGATA 5' RACE 3" round PCR 
IRF10-3-F1 AAGTGACCACCAAGAGCCCAGAT 3' RACE 1“ round PCR 
IRF 10-3-F2 CAATGGCTCGCCTTCTTTGTCA 3' RACE 2" round PCR 
B-actin F CAGTCCTCCTAAGGCGATAA а er 
B-actin R GCATCATCTCCAGCAAAGC Real-time quantitative PCR control 
IRF4-RT- F ACGGAGACAAGCCCAACAAG Real-time quantitative PCR 
IRF4-RT -R ACAAGTGGCTGCCCGTCTG 4 
IRF10-RT- F ACAATGGCTCGCCTTCTTT Real-time quantitative PCR 
IRF10-RT- R TGGGACCACTCCAATACAC 4 


5.0 (http://www.mega software. net/index.html/). Multi- 
ple sequence alignments were generated using Clutal W 
2.02. Phylogenetic trees were created by the neighbor- 
joining method using MEGA 5.01 and were bootstrapped 
1,000 times. 


Tissue distribution 

Intestines (D, HK, trunk kidney (TK), GI, brain 
(BR), liver (L), SP, muscle (M), gonad (GO), whole 
blood (BL) and skin (SK) from four fish were collected 
and used for RNA isolation using TRIzol (Invitrogen, 
USA) in order to analyze the expression of ma/RF4 and 
malRF10 in healthy eels. Real-time PCR was conducted 
to detect the expression of maIRF4 and maIRF10 in 
these tissues. The expression of maIRF4 and maIRF10 in 
different tissues was calculated relative to the expression 
level of B-actin. 


Real-time PCR analysis 

The maIRF4, maIRF10 and house-keeping gene В- 
actin cDNA fragments were generated via RT-PCR. A 
common reference with purified PCR products of the 
three genes was separately used for quantification. Prim- 
ers used for detection of gene expression are detailed in 
Table 1. PCR reactions were performed using Chromo 
4™ Continuous Fluorescence Detector (MJ Research). 
Amplifications were carried out at a final volume of 20 uL 
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containing 1 uL DNA sample, 10 uL 2xSYBR green 
Real time PCR Master Mix (Toyobo, Japan), 2 uL of 
each primer and 5 uL H,O. PCR amplification consisted 
of 5 min at 95 ?C, followed by 40 cycles consisting of 10 s 
at 94 °С, 20 s at 58 °С, 20 s at 72 °С and read plate at 82 °С. 
Melting curve analysis of amplification products was 
performed at the end of each PCR reaction in order to 
confirm that a single РСК product was detected. Each 
sample was run in duplicate to ensure accuracy. 


Challenge trial using different stimulus 

To study the effect of immunostimulants on the exp- 
ression of maIRF4 and malIRF10, four groups of fish 
(three fish each group) were respectively injected intra- 
peritoneally (i.p.) with either 500 uL 2 mg poly I:C/mL 
(Sigma), 500 uL 2 mg LPS/mL (E. coli O127:B8, Sigma), 
1x10* CFU A. Hydrophila/mL, or 500 uL PBS to serve as 
a control. Consecutively, at 0 h, 12 h, 24 h, 48 h after 
injection, three fish were killed, and tissues from the 
head kidney, gill and spleen were collected for total RNA 
extraction and gene expression detection via real-time 
PCR. 


Statistical analysis 
Fold change was calculated as (Ts/Tn)/(Cs/Cn) wh- 


ere Ts equals the treated sample assayed for the specific 
gene and Tn equals the treated sample assayed for the 
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normalizer gene (B-actin), and Cs and Cn respectively 
equals the calibrator group with the specific and normal- 
izing gene (Purcell et al, 2004). One way-analysis of 
variance (ANOVA) and an LSD post hoc test was used to 
statistically analyze the expression data, with P«0.05 
being considered statistically significant. 


RESULTS 


Cloning and characterization of maIRF4 and maIRF10 

Analysis showed that та!КЕ4 cDNA (GenBank 
accession no. JX463267) is approximately 1.7 kb in 
length, with an open reading frame (ORF) of 1 356 
nucleotides, encoding a protein of 451 amino acid (aa) 
residues with a putative molecular weight of 51.7x10? 
and an isoelectric point (pI) of 9.61. Furthermore, the 
malRF4 cDNA contains 22 microsatellite dinucleotide 
(T-G) repeats in the 3” UTR. The full-length cDNA of 
maIRF10 meanwhile comprises 1 744 nucleotides 
(GenBank accession no. JX463268). The ORF encodes a 
411 aa protein, with a calculated molecular weight of 
58.97x 10" and a pI of 8.85. The 3'-UTR has two ATTTA 
motifs that mediate mRNA degradation (Iwai et al, 1991) 
and the polyadenylation signal (AATAAA) lies 27 bp 
upstream of the polyA tail. 

Phylogenetic analysis showed that the IRF4 mole- 
cules diverged from the other IRF4 subfamily members 
(ІКЕ8, 9, 10), with eel IRF4 clustering with the other 
teleost IRF4s (Figure 1). Furthermore, teleost IRF4 dive- 
rged into two clades, with maIRF4 having the highest 
similarity with medaka IRF4-2 (Figure 1). The putative 
malRF4 protein was 45.6% to 84.0 % identical to IRF4 
proteins from mammals, chickens, frogs and other fish 
(Figure 2), while maIRF10 is 43.4% to 76.6% identical 
to its homologous proteins in mammals, birds, amphib- 
ians and other fish (Figure 3). Additionally, flounder 
IRF10 was the closest to maIRF10. MaIRF4 and malIR- 
F10 had the DNA-binding domain and IRF association 
domain (Figure 2, Figure 3), both of which were 
important domains in the IRF family (Lohoff & Mak, 
2005). 


Tissue expression 

Real-time PCR analyses showed that the expression 
of maIRF4 in all tested tissues was quite low, with 
relative abundant in intestine (0.21x10 7-fold), head 
kidney (0.16x10-fold) and trunk kidney (0.14x10 ?-fold) 
(Figure 4). However, constitutive expression of maIRF10 
transcripts was detected in all tissues studied. Compared 
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with p-actin| expression was predominantly in head 
kidney (7.2х10 3-Ғо14) and whole blood (6.9x 10 ?-fold). 
And a relatively low level of maIRF10 expression was 
detected in brain (1.1x10°-fold) and spleen (0.4х107- 
fold) (Figure 4). 


Induced expression 

In order to further investigate the differences pattern, 
malRF4 and maIRF10 expression induced by poly Г.С, 
LPS and A. hydrophila were examined. Stimulation by 
poly EC up-regulated maIRF4 expression at 12 h post- 
stimulation (hps) in head kidney (23.2-fold), spleen 
(15.3-fold) and gill (3.8-fold) tissue (P«0.05). Subseq- 
uently, expression increased, reaching a peak in head 
kidney (45.8-fold), spleen (27.5-fold) and gill (5.4-fold) 
at 24 hps (P«0.05) (Figure 5). Following LPS stimul- 
ation, та1КЕ4 mRNA levels were initially increased in 
head kidney and spleen at 12 hps, and remained elevated 
at 24 hps before decreasing to basal levels in spleen. In 
gill only at 24 hps was a significant increase seen. After 
A. hydrophila stimulation, maIRF4 transcript levels were 
enhanced in all three tissues studied at all time-points 
and peaked at 24 hps, with a fold change of 17.7 in head 
kidney, 13.6 in spleen and 6.51 in gill (Figure 5). 

The level of maIRF10 induction was always lower 
than that of maIRF4. MAaIRFI0 expression was 
significantly up-regulated at all times-points by poly СС 
in head kidney and in spleen but not in gills. The 
malRF 10 expression level peaked at 24 h in head kidney 
(5.3-fold) and spleen (5.8-fold). Following LPS 
stimulation, maIRF10 mRNA expression increased from 
12 hps to 48 hps (up to 3.6-fold) in head kidney, but 
significant up-regulation occurred at 24 h (3.6-fold) in 
the spleen, with no significant change in gills. Following 
A. hydrophila stimulation, maIRF10 transcript levels 
were also enhanced in head kidney and spleen at 12 and 
24 hps, with a peak fold change of 6.3 and 4.9 at 24 hps, 
respectively. Notably, levels remained elevated to 48 hps 
in head kidney. No significant changes were detected in 
malRF10 expression in the gills at any time-point 
following stimulation by А. hydrophila (Figure 5). 


DISCUSSION 


Phylogenetic analysis showed that teleost IRF4 
genes aggregated into two branches, suggesting that the- 
se two fish IRF groups arose through teleost-wide whole 
genome duplication (Figurel). However, only one IRF4 
gene was cloned in the Asian swamp eel. This may be 


Volume 35 Issue 5 


384 XU, et al. 


explainable by a twice genome duplication in the species 
while a third genome duplication occurred in other 
teleost (Zhou et al, 2002). The maIRF4 gene is highly 
conserved compared with flounder IRF4 and medaka 
IRF4-2, respectively exhibiting 84.2% and 83.7% iden- 
tity over the entire protein. The serine-rich domain is a 
target for virus induced phosphorylation that facilitates 
interaction with other IRF members and subsequent acti- 
vation of virus clearance signaling pathways (Sun et al, 
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2007). This domain is present in zebrafish IRF4-1, fugu 
IRF4-1, medaka IRF4 (Figure 2) and rock bream (О. 
fasciatus) (Bathige et al, 2012) but absent from IRF4 in 
many other species including humans, chickens, xenopus 
(Xenopus tropicalis) (Figure 2), rainbow trout (Holland 
et al, 2010) and eel, indicating significant evolutionary 
divergence of IRF4. All IRF10 genes from birds, 
amphibians and fish are closely clustered and closely 
related to IRF4 (Figure 1). Compared with these 
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Figure 1 Phylogenetic tree analysis of IRF4, IRF8, IRF9 an IRF10 among vertebrates 


A neighbor-joining phylogenetic tree of IRF proteins based on protein sequences analyzed with Clustal W and MEGA 5.0. Data were analyzed using Poisson 


correction, and gaps were removed by pairwise deletion. The degree of confidence for each branch point was determined by bootstrap analysis (1 000 times). 


The sequences of IRFs used for the analysis are derived from the GenBank and Ensembl Databases, with accession numbers added after the genes. 
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DNA Binding Domain (DBD) 
Human IRF4 MNLEGGGRGGEFGMSAVSCGNGKLRQOMLIDOIDSGKYPGLVWENEEKSIFRIPWKHAGKODYNREEDAALFKAWALFKGKFREGIDKPDPPTWKTRLRCALNKSNDFEELVERSQLDISD 120 
Chicken IRF4 MNLEPG----ECGMNSVSCGNGKLRQWLIDQIDSGKYPGLMENDEKSIFRIHQKHACKODYNREEDAALFKAWALFKGKFREGIDKPDPPTWKTRLRCALNKSNDFEELVERSQLDISD 116 
Chicken АЕА ----- ----MTPLSCGNGKLROWLIDQIDSGKYPGLVWENEERSIFRI GKQDYNREEDAALFKAWALFKGKYREG1DKPDPPTWKTRLRCALNKSNDFEELVERSQLDISD 107 
Fugu IRF4-1 ----LD- -YGGSGSSGNGKLRQWLIEQVDCGKYPGLWWENDEKSIFRI GKODYNRDEDAALFKAWALFKGKFREGIDKPDPPTWKTRLRCALNKSNDFEELVERSQLDISD 110 
Zebrafish IRFA-1 | mnLpGp------CIMSVSCGNGKLRWLIEQIDSGEYSGLYWMENDEKTIFRI GKODYNRDEDAALFKAWALFKGKYREGLDKPDPPTWKTRLRCALNKSNDFDELVERSQLDISD 114 
Zebrafish ІКЕ4-2 мсврвв------ RGMTGSSGNGKLRQWLIEQVDTGKY PGLYWENDEKS I FRI HWKHAGKODYNRDEDAALFKAWALFKGKFREGVDKPDPPTWKTRLRCALNKSNDFEELVERSQLDISD 114 
Fugu IRF4-2 - -EDGGMSVSCGNGKLRQWL INQIDSRKY PGL IWENADKS LFRI 4iWKHAGKQDYNREEDAALFKAWAMFKGKFKEGVDKPDPPIWKTRLRCALNKSNDEDELVDRSQLDISE 110 
Medaka IRF4-1 - ------GCGNGKLROWLIDOIDSRRYPGLVWENEEKTIFRIPWKHAGKODYNREEDAALFKAWAMFKGKYKEGIDKPDPPTWKTRLRCALNKSNDFDEIVERSQLDISD 103 
Flounder IRF4 GSGSGGNGKLROWLIEQVDCGKY PGLVWENDEKS I FRI PWKHAGKQDYNRDEDAALFKAWALFKGKFREGIDKPDPPTWKTRLRCALNKSNDFVELVERSQLDISD 114 
Medaka IRF4-2 -YGGSGSSGNGKLROWLIEQVDCGKY PGLVWENDEKS 1 FRI PWKHAGKODYNRDEDAALFKAWALFKGKFREGIDKPDPPTWKTRLRCALNKSNDFEELVERSQLDISD 114 
Eel IRF4 ---YGGSGSSGNGKLRQWLI DOVDCGKY PGLVWMENDERT I LRI PWKHAGKQDYNRDEDAALFKAWALFKGKFREGI DKPDPPTWATRLRCALNKSNDFEELVERSQLDISD 114 
ШЕКІТІТТІГІТІГІГЕЕГІГТІГІГЕТТЕГІГІГТІТІТІТІІТТІГІГІТТІРІТТІТІТТІРТІГІГІГІІТІТІТІТІТІТІТІТГІІТІТІТІГІГІ 
Натап ІКЕ4 PYKVYRIVPEGAKKGAK -QLTLEDPOMSMSHPYTMTTPYPSLPAQOVHNYMMPPLDRSWRDYVPDQ -PHPEIPYQCP-MTFGPRGHHWOGPACENGCQ 215 
Chicken IRF4 PYKVYRIVPEGAKKGAK NSMEE-QPLMNHPFPITSPYTSLPS-QVPNYMVP-| 'AARGHHWOGPGCENGCQ 209 
Chicken IRF4 PYKVYKIIPEGSKKGK---------------- -TSLIS-FFLLNVPLSYVTTGRSGLRQNINSLFLA----TMGDVRPP 173 
Fugu IRF4-1 PYKVYRIVPEGAKKTSAECORA! IILSAGSROEDSPLSPVS--YPVHPSYPALQTQ-VPOFMPIQEC-DWRDYCQE----AASMPEMSFTOCPCPP--RSLPWONPSMENGYO 220 
Zebrafish IRFA-1 PYKVYRIVPEGAKRG: IKAISMEENTTHVTPLS--YPMHSAYPALOPO-MSGFMLPQERRDWREFGSDPPETQTPHADLPYGQCPYPPS-RSLPWHTAPCDNGYO 214 
Zebrafish IRF4-2 PYKVYRIIPEGAKKG: LDSHKQDLRALSLFFFFLYLPCKSPSDENSTAYFQTFTEMKMSCIEEMPQHANAIGFISPYTSLRN-QFLKREGIGGITTLLOS-S 204 
Fugu IRF4-2 PYKVYRIVPEGSKKG -SRSIEDSQSNSGSPNYPMHPTYAPAPSQ-VCNYISPAER-GWREYP-- -TLSDISYSQSPYTS--R---WDP-----i GYQ 194 
Medaka IRF4-1 PYKVYRIIPEEAKKA( KMSIMKDAQSAPLGYIPPYTSLHTQVGGNPPSYVAPAFCLHPTISFKLES. )VHIKPETASLAGPGLYAKSGETGVEGLLATGT 213 
Flounder IRF4 PYKVYRITPEGAKK: -RPRQEDSPVSPMS--FQVHP-YPALQTQ-MPQYMTTPDG-SWRDFCPE: QAPLPELPYSQCPCPP--RSLSWOGPSMENGYO 202 
Medaka IRF4-2  PYKVYRIIPEGAKK: IRPROEDSPTSPVG--YOMHPSYPTMPTOEMTPYISSPEC-SWRDYCPE OTSLTELSFAQCPCPP--RSLPWOTPSIENGYO 209 
Eel IRF4 PYKVYRIIPEGAKK ---RPROEDSPISPMS--YOVPSSYPALOTQ-MPQYMAPVEC-GWRDYCQE----QASLPELPFTOCPCPT--RSLPWOSASIENGYO 203 
ТТТ 
Serine-rich domain 
Human IRF4 VTGTFYACAPPESQAPGVPTEPSIRSAEALAFSD--- NLDOVLFPYPEDNGORKNIEKLLS--- 309 
Chicken IRF4 VTGTFYACAPPESQTPGIPIEPSIRSGEALALSD-- -QGOSYEVS-----SLEQVIFPYPEDNSQORKNIEKLLS--- 303 
Chicken IRF4 LQEQLMLSNSINFRYSFVQFLSSINLTDVT--SD-- -HGHCYDIS-----RLEQVVFPYPEDSSLRKNIEKVLS--- 265 
Fugu IRF4-1 LRASIYSYSPADSQSSPFALDASIRSAETESD-- -FRLHVTVFYRDSPVREVTVSSPEGCHLT----PCSPEDKLYQQP----GNPEVVPLPVDSLSPPR-TDDdSSSPPS3 319 
Zebrafish IRFA-1 ISGSFYTYSP--SESHPVAMDPSMRSAEAMAISD-- ---CRLHVSLFYRESLVKELTTSSPEGCRISSSASPGSPSSPSSPSEERLYGGAEPVLFPFPYPOSQRRGAEKLPN--- 319 
Zebrafish ІВЕ4-2 KLCPVHIYMGHTQNCSTASVPTHHQSAELGPGHTQKMVRWSFSYFSLHVALYYRESLVKEVTVTNPEGCRIT[SS-STCSPSSPSSCPEDKFHSGAEVVFFPFPYPESQROGAEMLPN--- 320 
Fugu IRF4-2 FSGSFYSCNASDPQPSPFTLDTSMRSAEAMALSD---------- YRLHVMVFYRDALVREVTVSSPEGCQLG----PSREG-QAYASP----GAPELVELPH 
Medaka IRF4-1 AASPSCTSPRLTCRPTIRSVP---LPAAVOQPTLVHIACRK-----CLHVSLYCRETLVKEVTTTSPEGCRIK- 
Flounder IRFA — LRASIYSYGPADSQAPPFTLDAGIRSAEALSD FRLHVSVYLRDNLVREVTTSSPKGCHIT----PCSPEEKLSLLP----GGPDVVPLPVDHLSVORRAEECSPNPPS 
Medaka IRF4-2 LRASIYSYGP-DSQPSPFSLDASIRSAEALSATY LRLHVSVYLRETLVREATVTSPEGCHLT----PCPPE-KHYLTA----GGPEVVPLPVDSLSALRRTDECPPSPPS 309 
Eel IRF4 LRASIYSYGPADSQPSPFSLDASIRSAEALPD --FRLHVTIYFRDTVMREVTTSSPEGCHIT----PCSPEEKHYLTP----GGPEVVSLPVDSLSAQRRADECPPSPPS 303 
Жк: ој: ка f. * z "EET 
IRF Association Donk) = `7 Бис ји 
Human IRF4 HLERGVVLWMAPDGLYAKRLCQSRIYWDGPLALCNDRPNKLERDQTCKLFDTQQFLSELQAFAHHGRSLPRFQVTLCFGEEF PDPQRORKLITAHVEPLLARQLYYFAQQNSGHFLRGYD 429 
Chicken IRF4 HLERGVILWMAPDGLYAKRLCOSRIYWDGPLALCSDRPNKLERDQTCKLFDTOOFLAELQAFAHHGRPLPRYOVALCFGEEFPDPORORKLITAHVEPMFAROLYYFAQONSGHLLRGYD 423 
Chicken IRF4 HLDRGVILWLAPDGMYARRMCQSRIYWEGPLSMYSDRPNKLERDIACKLFDTQOFLSELQAYALHGRPVPRFQIVLCFGEEFPDAQRGRKLITAHVEPVFARQLV--------------- 370 
Fugu IRF4-1 TLESGVLLWMGLDGLYARRLCOSRVYWORGSAQYADKPNKLEREVNYKLLHTQDCLAEIQSFGLHGRPPPRFQILLSFGDECLDPORHRRTLAVQVEPLFSROLLYYAQPVGGHYYRNYE 439 
Zebrafish IRF4-1 | VLERGVLLWLSPDGLYAKRLCOGRVYWEGPLAPYADKPNKLEKEQTCKLMDTOQFLTELOGFIHHGRPMPRSQVILCFGDEFPDPORQSKMITAQVEPMFAROLLYFASQTNGHYLRSYE 439 
Zebrafish IRF4-2 LEKGVLLWMTADGMYAKRLCQGRVYWEGPLAQFMDKPNKLEKEQPCKLFDTHOFLIDLODFAHNGRHLPRHOVVLCFGDEYPDPQRSRKMITAQVFKIFHKNMDSFTHLLNSSSIHSS- 439 
Fugu IRF4-2 -LERGVILWMAPDGLYARRCCPCRVYWTGAHAPPTDKPNKLEREQNCKLLDTHLFITELOSYTLHARPAPCSQVLLFFEDESTEGQRPRRTYTVQVEPLFAROLLILTHPGSMNYIRSHE 401 
Medaka IRF4-1 ^ VLETGVLLWVTPDGLYAKRLCRGPVYWEGPTAPYSDKPNKLEKAQPCKLLDTQOFLVELODFAOHGRHVPRHOVALCFGDEHLGPORPRKMFTAQV------------------------ 402 
Flounder IRF4 TLERGVLLWMGADGLYACRLCQSRVYWQGGPSPYGDKLNKLERDVTCKLLHSQDYLTELQSFGLEGRPLPRLQVLLSFGDECLDPORQRRTLSVQVEPLFARQLLYYAQQTGGHYYRSYD 422 
Medaka IRF4-2 ^ VLERGVLLWMSADGLYARRLCOSRVYWQGGLSPYGDKPNKLEREVTCKLLETQDFITEIQSFGFEGRPPSRFQVLLSFGDECLDPQRQRRTLTVQVEPLFARQLLYYAQOTSSHYYRSYE 429 
Eel IRF4 TLERGVLLWMGADGLYARRLCQSRVYWQGGLSPYGDKPNKLEREVTCKLLHTQGYLTEIQSYGLHGRLLPRFQVLLGFGDECLDPQRORRTLTVQVEPLFARQLQYYAQOTGSHLYRSYE 423 
ET PESEN, EERE OR CIR a ж: KER, SAGs ae SE es e Se ao CRIT 
Human IRF4 
Chicken IRF4 
Chicken IRF4 
Fugu IRF4-1 HAGVT--EHVSPSEDYQRTITHHHSS: SLOE 467 
Zebrafish IRF4-1 -SLPVEDYORSLOHLTE----------- 460 
Zebrafish IRF4-2 ------------------------- 
Fugu IRF4-2 LOHLPPEHSLSPTODYHRVITHHHNS------GPQN 431 
MedakaIRF4-] ------------------------------------ 
Flounder IRF4 LPGVT--DHFNASEDFQRVVTHHHHHSSSSSSSLOE 456 
Medaka IRF4-2 ^ HPGVT--DHINASEDYORAITHHHAAA----- SLOE 458 
Eel IRF4 HPVVT--EHINSSEDYORAITHHHSS------ SLOE 451 


Figure 2 Multiple alignment of maIRF4 amino acid sequence with other IRF4 proteins 


Symbol (*) represents identical residues, (:) conservative substitution and (.) similar residues. Missing amino acids are denoted by hyphens. The DNA-binding 


domain (DBD) and IRF association domain (IAD) are highlighted in grey. The conserved tryptophan (W) residues that comprise a “tryptophan cluster" are 


boxed. A serine-rich domain is also boxed. 


IRF10 family members, the maIRF10 protein sequence 
showed the greatest homology with Japanese flounder, 
medaka (Oryzias latipes) and stickleback (Gasterosteus 
aculeatus), ranging from 59.6% to 76.6% identity across 
the whole protein. 

Generally, /RF10 is expressed in all tissues while 
IRF4 has a much more limited expression pattern. For 
example, among chickens /RF10 is detected in all tissues, 
while only faint signals of /RF4 mRNA were detected in 
all the previously tested tissues in chicken or mice 
(Dougherty et al, 2009; Takaoka et al, 2008). Consistent 
with these results, maIRFI0 mRNA was found to be 
abundant in all the tissues we analyzed, with predomi- 
nant expression in whole blood, skin and head kidney, 
while maIRF4 expression was significantly lower than 
that of maIRF10 in all the same tissues (Figure 4). The 
highest expression of maIRF4 was detected in the intes- 
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tine, head kidney and gill, major immune tissues/organs, 
trunk kidney and brain. A previous study reported that 
among the Japanese flounder, 7RF/0 mRNA was stro- 
ngly expressed in gill, head kidney, heart, peripheral 
blood lymphocytes, spleen and trunk kidney (Suzuki et al, 
2011), while another study on rainbow trout found the 
highest expression levels of /RF4 in lymphomyeloid-rich 
fish tissues (e.g., spleen, head kidney and gills; Holland 
et al, 2010). MaIRF 10 also shows relatively high expre- 
ssion in non-immune tissues such as muscle, gonad and 
brain, indicating that IRF10, as well as IRF4 and IRFS, is 
present in at least some non-immune tissues (Dadoune et 
al, 2005; Eguchi et al, 2008). 

Previously, it was reported that IFN1, IFN-y and 
ConA can up-regulate IRF10 expression (Nehyba et al, 


2002). However, induction expression of IRF10 occurs 
relatively late, and is dependent on protein synthesis. 
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DNA Binding Domain (DBD) 
MNLEGGGRGGEFGMSAVSCGNGKLROWILIDOIDSGKYPGLMWENEEKSIFRIHWKHAGKODYNREEDAALF| 
MNLETGSRGSEFGMSAVSCGNGKLROWILI DOIDSGKY PGLYWENEEKSVFRIAWKHAGKODYNREEDAALF! 
- NAQIESGRYAGLRWEDAGKTLFRIRHWKHAAKQGYQVRQDAALF| 
IAQIDSGRYPGLHWENRERTLFRIHWKHAAKQDYRQOODAALF| 
LEQIDSGRYPGLHWENPEKTVFRIHWKHAAKODYRQOADAALFE! 
IAQIDSGKYAGLSWENEEKTMFRIPWKHAAKODYRONODAAL] 
IAQIDSAEYPGLSWENAEKSMFRIFHWKHAAKQDYRONODAALF| 
IAQIDSGKYDGLHWEDERKSMFRIPWKHAAKKDYKEEEDAALF: 
VAQIESGRYEGLQWEDEAKTMFRIHWKHAAKKDYKOTEDAALF| 
ISOIESGRYEGLSWEDEDRTMFRIPUWKHAAKKDYKQTEDAALF| 


ЕКСКЕКЕСІРКРРРР! ‘TRLRCALNKSNDFEELVERSQLDISD 120 
LFKGKFREGIDKPDPP' ‘TRLRCALNKSNDFEELVERSQLDISD 120 
IYKGKYLEGIDKEDPS' ‘TRLRCALNKSADFCEVPERSQLDISN 106 
\VYKGKYHEGTDKADPS' TRLRCALNKSTDFQEVPERSQLDISE 106 
[YKGKFREGSGRADPS TRLRCALNKSPDFOEVMENSOLDLTE 104 
KGKFOEGRDKADPS TRLRCALNKSTDFOEVPERSOLDISE 106 
IYKGKFQEGRDKADPS' ‘TRLRCALNKSTDFQEVSERSQLDISE 106 
\VYKGKFREGRDRADPT! ‘TRLRCALNKSTDFREVPERNQLDISE 103 
\VYKGKYKEGRDRADPT! TRLRCALNKSTDFOEVPERSOLDISE 106 
\VYKGKY IEGRDKADPT! ‘TRLRCALNKSTDFQEVPERNQLDITE 105 
IAQIESGKYAGLGWENEEKTMER I HWKHAAKKDY SOTADAALFI \VYKGKYREGRDKADPT! ‘TRLRCALNKSTDFEEVPERNOLDITE 99 

IGOIESGRYEGLSWENEEKTMFRI HAKHAAKKDYROTADAALF| WVYKGKYIEESDHOEPTI ‘TRLRCALNKSTDFQEVPERNQLDITE 105 


215%; HELE, X XE kE, ;;;Жжжжжжжж,ж; ,, KAKAK DU DIO Ж Cl DE. ЖЖЖЖ жж ж; ж, ЖЖЖ; 























































-DYVPDOPHPEIPYOCPMTFGPRGHHWQGPACENGCQVTGTFYA 222 
DYAPDOSHPEIPYQCPVTFGPRGHHWQGPSCENGCQVTGTFYA 


PYKVYRIVPEGAKKGAKOLTLEDPOMSMSHPYTMTT PY PSLPAQOVHNYMMPPLDRSWR- 
PYKVYRIVPEGAKKGAKQLTLDDTOMAMGHPYPMTAPYGSLPAQOVHNYMMPPHDRSWR: 
PYKVYRLLSDGAHDPVTRACAPGKEE--GILQ--SQOKPPGGVTEPACR-TEEQVLS: 
PYKVYQIVCDG: 
PYKLYRILTPSE 
PYKVYRILDDSGR- 
PYKVYRILEDSAR- -KALRETHLSMQED----ESDSGAGWTVNGTASVCSGTGAK 
PYKVYLIOQDGGS- -KTROESHQS-HORLSWSWSGSVHK-------- OYTMSHFSS 163 
PYKVYRVKPFPRSPSCPGKQTHFDENFDIEATIAYTCRRLDPCTSSFSLSELAAGSKLYGWITPRNSSLCCPKKRPPGPKGRVPNNQXXXXSLTLFLVASVSKGSVSDKQRRKTISCHNK 226 



























PYKVYRIOQDS--- SPNSPDILDEKRPFQNESFQANIEEEK-TWHVDLMSEHMYCDIK 179 
PYKVYRVSSSGEP- PENPVAPVKQESTSSGEELK-EINEMLMKEHMYCGSS 163 
PYKVYCIOQDNPPKT- KLEIKSCILSKESPQANDQVTVOG- SPEKPSFLDREVHYQKESFKAKEEEEK-PVSGDLMREHMYCELR 187 
закра: 

IRF Association Domain (IAD) 
CAPPESQAPGVPTEP -SIRSAEALAFSDCRLHICLYYREILVKELTTSSPEGCRISHGHT-- -YDASNLDQVLFPYP- -EDNGQRKNIEKLLSHLERGVVLWMAP 32] 
CAPPESQAPGIPIEP: -SIRSAEALALSDCRLHICLYYRDILVKELTTTSPEGCRISHGHT -YDVSNLDQVLFPYP- -DDNGORKNIEKLLSHLERGLVLWMAP 32] 















FLLPSPLADHSYQAONPVHPWSPLPSEDFSNPDCWLEVRLFYRAELVREATALAAEGCSLSPRAA -TAAAERLLGPPPRVAQVRFPEPPASAHVLRRLLPHLERGVLLWVAP 302 
HLLP-- КАР-. -SFLPAEDVNHSDCWLHIRLYYCDVLVKELTTRTAEGCRIASRTEG- PYGPSCMEQIEFPPPRALGGGG-WIVAVTEVLERLLPHLERGVLLWVAP 285 
ISIPPIQIQVRKROGSPEPYGSKSWKSDRLCPPKATTOHSYYRALDLSKVTTKTAEGCRIVPWSAAPQOCGPCPPLSLEVVSLPSPOELPGHL--TAEMGGLIQKLLWHLERGVLLWVAP 318 
--YPSHVPIS-DCRLEVRLFYHGNLVQSLTTASPDGCFILOGCAPVGNERIYGPCEAEKVFFPRP: -DTIRLPPGIAEAMSRLLPHLEKGVLVWVAP 280 
YPPOVHLSVDCRLELRVFYYGRVVESVLSCSPDGCFLLOGCAPVGSERIYGPCAATQLFFPPP: -NAAMLPTGICEAMTRLLPHLEKGVLLWVAP 275 












































DKKPQSQVPSPI- 





-TFFSP-LTI SDERMOVMLLYQGORVHKVITKS PDGCFILOGRVPLGNERIYGPCNAQOLSFPSP----ASLSLPSHVAETMARLLCHLERGVLLWVAP 292 


ЖЕГЕН гж ж 2 IX жжжржррржраж 


DGLYAKRLCOSRIYWDGPLALCNDRPNKLERDQTCKLFDTQOFLSELOAFAHHGRSLPRFQVTLCFGEEFPDP--QRORKLITAHVEPLLARQOLYYFAQONSGHFLRGYDLPEHISNPED 439 
DGLYAKRLCOSRIYWDGPLALCSDRPNKLERDOQTCKLFDTOOFLSELOVFAHHGRPAPRFOVTLCFGEEFPDP--QRORKLITAHVEPLLARQLYYFAQONTGHFLRGYELPEHVTTP-D 438 
EGVFAKRLCQGRVYWRGPLAPHRAQPNKLERERTCOLLDTRRFLAELRAHLODGCPEPEYOIRLCFGEEYPGPPDOPEERLIMAHVEPVFARELLLHSKIHGQSLEGASAQQAGFGLALM 422 
EGVFMKRQCQGRVYWNGPLAPHODWPNKLEREKTYKLLDTQQOFLEELRRYLSHGQPAPQYQIHLCFGEEYPTSTGRHLOKLIMAHVEPVFARELFHHAQRLGPALLRGCPQPCAP - 400 
EGVFIKROCNARVYWSGOLAPHSDRPNKLEREKTCKVLDTEQFLOELOMYVTOGTPEPQYOIQOLCFGEEYPDTSRVNPKILMMALV----------------------------- 
DGVFIKRFCOGRVYWDGPLAEHROKPNKLERERTCKLLDMTIFMQELOSHQOATGPEPRYTVDLCFGEEFPDPSQPKNKKLITAQVIPLFAVECLRRHNASNNVEMKOSPPHRKTND- 
DGLFIKRFCOGRVYWDGPLAEHRHKPNKLERERTCKLLDMKIFSQELLNYROGIGPEPQYIVELCFGEEFPDPTOPKNKKLIRAQVTPMFAVDALRKLKADNNVEMKPPHPLAQENQ: 
DGVFIKRFCQGRVYWSGPMARHADRPNKLEREKTFKLLDVPAFLNDLRGCLOGKGPAPSYEIELCFGEEYPDPNVPKSRKLIMAKVVPLFAVELLER----------------------- 
DGVFIKRFCOGRVYWSGPMAEHTNRPNKLEREKTFKLLDIPIFLHGLOSCLOGKGPPPPYKIELCFGEEYPDPGVPKSKKLIMAQVVPLFAEELLQRFN--------------------- 
DGVFIKRFCOGRVYWSGPLAOHTDAPNKLEREKTFKLLDIPRFVSELORSLWGKGPAPSYEIELCFGEEYPDPHVVKTRKLIMAQVVPLFAVELLOKFNPGASEEKRSNLSSNSVGEKL- 404 
EGLFIKRFCOGRVYWSGPLARNTDOPNKLEREKTFKLFDVNTFFTELOSYLRGKGSKPSIEVELCFGEEFPDPKLPKNRKLIVAQVIPSFALELLORFSLGGSADRSALLFVKLDENDFS 389 


EGVEIKRECOGRVYNSGEMALHTDRPNKLEREKTFKLLDIPTELSEFOTCIQGRGPSPSCHEIELCPGEEYPDPSVEKTRKLIIAQUVVPLFRVELLOKFNLVETEENQUTLSSNTOREKM- 411 
i*i: жижу Kp** * рж ******llk Dll г жжжжжжрж | жр ж ж 
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Figure 3 Multiple alignment of maIRF10 amino acid sequence with other IRF10 proteins 


Symbol (*) represents identical residues, (:) conservative substitution and (.) similar residues. Missing amino acids are denoted by hyphens. The DNA-binding 


domain (DBD) and IRF association domain (IAD) are highlighted in grey. The conserved tryptophan (W) residues that comprise a “tryptophan cluster" 


boxed. 
0.30 IRF4 " 
8 6 IRF10 
B 0.25 8 8 
ЗА E 7 
Zo 020 өт 
о = HS 6 
35 55 
8 6015 4-5 
9% ef 4 
og oo 
$4 010 283 
25 
8 0.05 RO 2 
ш : E: 1 
0 0 
I HK TK GI BR L SP М GO BL SK I HK TK GI BR L SP M GO BL SK 
Tissue Tissue 


Figure 4 Expression analysis of maIRF4 and maIRF10 in different tissues 


are 


Real-time quantitative PCR was performed with cDNA samples prepared from whole blood (BL), brain (BR), gill (GI), skin (SK), muscle (M), intestines (1), 


spleen (SP), liver (L), gonad (GO), head kidney (HK) and trunk kidney (TK). Bars are averages plus standard deviations of tissue data from 4 fish. 


Neither ІЕМІ or IFN-y induce IRF4 expression, though 


IRF4 is known to be 


CD3, phorbol-12-myristate-13-acetate (PMA) and IgM 
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cross-linking in chicken and mouse (Ma et al, 2006; 


induced by ConA, plant lectins, Nehyba et al, 2002). In humans, dengue virus and T cell 


receptor cross-linking has likewise been known to 
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Gene cloning and induced expression pattern of IRF4 and IRF10 in the Asian swamp eel (Monopterus albus) 387 
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Figure 5 Expression levels of maIRF4 and maIRF 10 induced by poly I:C, LPS and Aeromonas hydrophila in head kidney, spleen and gill 


Fish were injected with each stimulant and 12 h, 24 h and 48 h, and tissues were then isolated for RNA extraction and cDNA synthesis. The transcript levels 


were obtained using real-time PCR and expressed relative to the 0 h healthy (control) fish and presented as the fold change after being normalized to the 


expression of B-actin. Bars denote averages plus standard deviations of tissues from 3 fish. * over the bar indicates a significant difference after a treatment by 


one-way ANOVA (7: P«0.05). 


increase IRF4 mRNA levels (Chen et al, 2008; Matsu- 
yama et al, 1995). Among the Asian swamp eel, we 
found that poly Г.С, LPS and А. hydrophila induced 
malRF4 and | maIRF10 
intraperitoneal injection. of poly LC significantly 


expression, wherein 
induced maIRF4 and maIRF10 expression from 12 hps 
to 48 hps (Figure 5). Although maIRF4 and maIRF10 
expression induced by LPS or A. Aydrophila was lower 
than that induced by poly I:C, there was a significant 
increase in expression in both the head kidney and 
spleen following both treatments (Figure 5). maIRF4 
expression peaked at a 45.8-fold increase in head 
kidney and maIRF10 reached a maximal expression 


level (5.8-fold) at 24 hps following poly I:C stimulation. 


Similar results have been found in other fish: for 
example, in rainbow trout, IRF4 was significantly up- 
regulated by PMA (4-fold above control levels) in 
splenocytes in vitro (Holland et al, 2010) while in the 
Japanese flounder, IRF10 expression was up-regulated 
7.8-fold at 6 hps by LPS and peaked at a 23.6-fold up- 
regulation at 6 hps following poly EC stimulation of 
peripheral blood leukocytes (Suzuki et al, 2011). 
Further comparison of IRF10 and IRF4 inducible 
expression highlighted a greater increase in malRF4 
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expression as compared with that of maIRF10. These 
results are consistent with previous expression studies 
conducted on birds (Nehyba et al, 2002), which found 
that IRF4 constitutively expressed at a low levels in most 
tissues and organs among chickens (Dougherty et al, 
2009) and that inducible expression of IRF10 requires 
protein synthesis, with IRF10 playing a unique role in 
the later stages of antiviral defense (Nehyba et al, 2002). 
By contrast, in rainbow trout, LPS significantly down- 
regulated IRF4 expression in splenocytes (Holland et al, 
2010). Similarly, estrogen was also found to decrease the 
expression of IRF4 and IRF10 in chickens (Dougherty et 
al, 2009). 

In summary, in this study the complete та КЕ4 and 
malRF10 genes of the Asian swamp eel were cloned for 
the first time. Analysis of the overall amino acid 
sequence, identity and phylogenetics, showed that the 
cDNA encoding malRF4 and maIRF10 exhibited 
homology with other IRF4 and IRF10 proteins of tetr- 
apod species and fish sequences currently deposited in 
public databases. Furthermore, maIRF10 was constit- 
utively expressed at high levels in all the tissues invest- 
igated, while the transcript level of maIRF4 was found to 
be quite low in all organs. Stimulation by Poly Г.С, LPS 
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and A. hydrophila were also shown to induce malRF4 
and та1КЁ10 expression, indicating that these two fac- 
tors likely play a part in the immune reaction to both 
bacterial and viral infections. 
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Abstract: The B cells translocation gene 1 (BTG/) is a member of the BTG/TOB family of anti-proliferative genes, which have 
recently emerged as important regulators of cell growth and differentiation among vertebrates. Here, for the first time we cloned the 
full-length cDNA sequence of Hyriopsis schlegelii (Hs-BTG1), an economically important freshwater shellfish and potential 
indicator of environmental heavy metal pollution, for the first time. Using rapid amplification of cDNA ends (RACE) together with 
splicing the EST sequence from a haemocyte cDNA library, we found that Hs-BTG/ contains a 525 bp open reading frame (ОКЕ) 
encoding a 174 amino-acid polypeptide, a 306 bp 5' untranslated region (5' UTR), and a 571 bp 3' UTR with a Poly(A) tail as well as 
a transcription termination signal (AATAAA). Homologue searching against GenBank revealed that Hs-BTG1 was closest to 
Crassostrea gigas BTGI, sharing 50.57% of protein identities. Hs-BTGI also shares some typical features of the BTG/TOB family, 
possessing two well-conserved А and B boxes. Clustering analysis of Hs-BTG1 and other known BTGs showed that Hs-BTG1 was 
also closely related to BTG1 of C. gigas from the invertebrate BTGI clade. Function prediction via homology modeling showed that 
both Hs-BTGI and C. gigas BTG1 share a similar three-dimensional structure with Homo sapiens BTG1. Tissue-specific expression 
analysis of the Hs-BTG/ via real-time PCR showed that the transcripts were constitutively expressed, with the highest levels in the 
hepatopancreas and gills, and the lowest in both haemocyte and muscle tissue. Expression levels of Hs-BTG1 in hepatopancreas 
(2.03-fold), mantle (2.07-fold), kidney (2.2-fold) and haemocyte (2.5-fold) were enhanced by cadmium (Cd”") stress, suggesting that 
Нѕ-ВТСІ may have played a significant role in Н. schlegelii adaptation to adverse environmental conditions. 


Keywords: Hyriopsis schlegelii; BTG1; Gene cloning; mRNA expression; Cadmium stress 


Exposure to heavy metals—cadmium in partic- 
ular—is widely known to be toxic to humans and higher 
vertebrates, negatively affecting several key organs. 
Unfortunately, industrial pollution run-off often results in 
these metals making their way into water sources, where 
they are absorbed by the ecosystem's resident species. 
For example, oysters were previously found to 
accumulate environmental cadmium, suggesting that 
they may serve as potential indicator organisms for 
marine cadmium contamination (Lu et al, 1998). 
However, the toxicological activities and effects of such 
heavy metals among invertebrates such as freshwater 
shellfish are not well understood. This situation is 
especially acute in developing countries and regions such 
as China, where heavy levels of water pollution that has 
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accompanied industrial development. 

One key reason that shellfish species like oysters 
can serve as viable indicators of environmental heavy 
metals is that they have remarkable adaptability to 
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different environments and are resistant to the toxic 
effects of metals like cadmium. Among freshwater 
shellfish, both the Japanese Hyriopsis schlegelii (H. 
schlegelii) and Chinese H. cumingii species of mussels 
exhibit similar characteristics. In 1997, H. schlegelii 
was introduced to China by the Reservoir Development 
Company in Fuzhou City, Jiangxi Province. Over the 
last decade, this species has become a key economic 
feature of freshwater shellfish production, largely due 
to decades of selective breeding that have results in 
strong disease resistance and environmental adaptability, 
as well as its enlarged shell width and other physical 
characteristics contributing to senior pearls cultivation 
(He et al, 2013; Peng et al, 2012a, b; Xie et al, 2011). In 
order to develop a better ecological farming, exploring 
molecular indicator capable of monitoring enviro- 
nmental pollution is particularly important for this 
mussel. 

Testing for the presence and effects of cadmium 
among potential indicator species is not always straigh- 
tforward, but gene expression analysis under various 
environmental factors stress is a commonly employed 
experimental method to explore candidate gene’s corres- 
ponding function (Duan et al, 2013; Wang et al, 2011) 
and in doing so illustrating the effects of certain 
environmental pollutants on different organisms. In the 
present study, we sought to examine the Н. schlegelii В 
cell translocation gene 1 (BTG1), which has been cloned 
and studied in several species, including Rattus norveg- 
icus, Mus musculus, Gallus gallus, Xenopus laevis, and 
Danio rerio (Fu et al, 2012), but not among many lower 
animals. BTG/ is a member of the BTG/TOB family that 
was identified from B lymphoblastic leukemia with 
chromosomal translocation (Rouault et al, 1992). The 
BTG/TOB family is a class of anti-proliferative proteins 
involved in negative regulation of cell cycle (Winkler, 
2010) and tumor cell growth (Mauxion et al, 2009). 
BTG/TOB proteins have a conservative BTG region of 
100 — 120 amino acids at its amino-terminal (N terminal), 
containing two well-conserved A and B boxes (Mauxion 
et al, 2009; Winkler, 2010). Among mammals, the 
BTG/TOB family includes six members: BTG1, BTG2/ 
PC3/Tis21, BTG3/ANA, BTG4/PC3B, TOBI/TOB, and 
TOB2 (Winkler, 2010), with both BTGI and BTG2 
belonging to the same subfamily with an approximate 
carboxyl terminal (C terminal) length (Mauxion et al, 
2009). Previously, the BTG/TOB proteins were reported 
to be involved in several different physiological activities, 
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including regulation of embryonic development, cell 
differentiation and apoptosis, among others (Mauxion et 
al, 2009). Similarly, BTG/TOB proteins were also found 
to affect gene transcription and mRNA stability in cells 
via deadenylation (Yang et al, 2008; Mauxion et al, 
2009). In this study, we obtained the full length cDNA of 
Hs-BTGI by library screening and RACE cloning to 
investigate whether ВТСІ is involved in the response 
against heavy metal stress. 


MATERIALS AND METHODS 


Experimental materials and reagents 

We obtained healthy H. schlegelii from the national 
seed market of Hyriopsis at Reservoir Development 
Company in Fuzhou City, Jiangxi Province. The H. 





schlegelii, with a shell length of 107-6.5 mm, was 
cultured in a laboratory aquarium for one week (18-25 'C) 
prior to testing. 

RNAase inhibitors, DNA Polymerases and Markers 
were purchased from TaKaRa. Gel extraction kit and 
plasmid extraction kit were purchased from Axygen, 
while DNase I (RNase Free) and a first strand cDNA 
synthesis kit were purchased from Promega. TRIzol was 
from Invitrogen. E. coli DH5a cells were stored in our 
lab. Primers were synthesized by Shanghai Sangon 
Company. Cadmium chloride (СаСІ,:5Н;О, AR) was 
produced by the Shanghai Chemical Reagent and other 
conventional reagents of analytical grade were obtained 
from Sinopharm Group. 


cDNA 

Total RNA of different tissues was extracted using 
TRIzol prior to assessing its purity and concentration. 
Primary sample digestion by RNase-free DNase I was 
performed before cDNA synthesis. The organization- 
specific cDNA synthesis was conducted as described 
previously (He et al, 2013; Peng et al, 2012a, b). 


cDNA library screening 

We performed PCR amplification of monoclonal 
colony using random primers (M13) according to our 
established cDNA library from blood cells of Н. 
schlegelii (Xie et al, 2011), under the following 
conditions: 94 °С 5 min; 94 °C 30 s, 45 °C 30 s, 72 °C 
90 s, 33 cycles; 72 °С 7 min. Amplified fragments more 
than 500 bp in length were sequenced using M13 
sequencing primers (Sangon), and then we obtained 
expressed sequence tags (ESTs) sequence of Hs-BTG1. 
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Full-Length cDNA cloning and molecular character- 
rization 

The 3' and 5' cDNA ends were amplified using 
SMART TM RACE Amplification Kit and Advantage 2 
PCR Kit. Primers (BTG/-GSP1 and BTG/-GSP2; shown 
in Table 1) were designed using Primer Premier 5.0 
according to the ESTs sequence of Hs-BTG/. Touchd- 
own PCR amplification was performed using universal 
primers (UPM), 5' specific (BTGI-GSP1) and 3' specific 
(BTGI-GSP2) primers, with 5'-RACE-Ready and 3'- 
RACE-Ready cDNA as the template. PCR reaction 
system was as follows: 10x Advantage 2 PCR Buffer 
5 uL, dNTP Mix 1 uL, 5'-RACE-Ready cDNA or 3'- 
RACE-Ready cDNA 2.5 uL, ОРМ 5 uL, BTGI1-GSPI or 
BTGI-GSP2 1 uL, Advantage 2 Polymerase Mix 1 uL, 
PCR-Grade water 34.5 uL, total volume 50 uL. Reaction 


conditions were as follows: 94 °C 5 min; 94 °C 30 s, 72 °C 
3 min, 5 cycles; 94 °C 30 s, 70 °C 30 s, 72 °C 3 min, 5 
cycles; 94 °С 30 s, 68 °С 30 s, 72 °С 3 min, 28 cycles ; 
72°С 10 min. Amplified products were recovered and 
purified using a gel extraction kit before being sequenced 
by Shanghai Sangon Company. 

NCBI ОКЕ Finder and BlastP were used for ident- 
ifying the largest ORF and homologous amino acid 
sequences of Hs-BTG. SignalP 4.0 was used for 
predicting the signal peptide. ClustalW was used in 
multiple comparisons of the amino acid sequences from 
different species. Protein homology modeling was 
conducted by using SWISS-MODEL. Chimera 1.8 was 
used for the tertiary structure analysis. The phylogenetic 
tree was constructed using Neighbor-Joining (NJ) 
method with 1000 bootstraps in Mega 4.1. 


Table1 Primers used for cloning Hyriopsis schlegelii ВТСІ cDNA and expression quantification 


Primer name 


Sequence (5'-3') 





MI3F 
MI3R 
SMART ШМ A Oligonucleotide (12 umol/L) 


3' -RACE CDS Primer А (3'-CDS; 12 pmol/L) 


5' -RACE CDS Primer A (5' -CDS; 12 umol/L) 


10xUniversal Primer A Mix (UPM) 
Long Primer (0.4 umol/L) 


Short Primer(2 umol/L) 
Hs-BTG1 5GSP1(10 pmol/L) 
Hs-BTG1 3GSP2(10 pmol/L) 
Н8-ВТСІ qPCR-F 

Hs-BTG1 qPCR-R 

f-actin qPCR-F 


В-асіп qPCR-R 


Expression profile of Hs-BTG1 

Totally, we tested 10 tissues from three mussels, 
including intestines, gills, gonads, mantle, ax foot, 
adductor muscle, blood cells, hepatopancreas, kidneys, 
and hearts. Total RNA was extracted and treated by 
RNase-free DNase I (Promega). 1.5 ug total RNAs were 
reverse transcribed using MMLV RT Kit (Promega) 
under following reaction conditions: 37 °C, 60 min, 75 °C 
5 min, 4 °C, hold. Expression profiles of 8701 from 
different tissues were detected by quantitative PCR using 
SYBR Premix Ex Taq П kit (TaKaRa) with f-actin as an 
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TGTAAAACGACGGCCAGT 
CAGCAACAGCTATGAC 
AAGCAGTGGTATCAACGCAGAGTACGCGGG 
(Т),5У № (N=A, С, G, or T; V=A, G, or C) 
AAGCAGTGGTATCAACGCAGAGTAC(T)30V N 


(N=A, C, G, or T; V=A, С, or C) 


CTAATACGACTCACTATAGGGCAAGCAGTGGTA 
TCAACGCAGAGT 


CTAATACGACTCACTATAGGGC 
CGTGGTCCTCGTGCGAAATAGAATG 
TGCCACTTTGGTTTCCAGTTGAGG 
GCTGGCTTATCGTGTGGAAT 
CGGTACTGGCATCCTCTGAT 
AAGGTTACGCCCTTCCTCAT 


GCCATTTCCTGCTCAAAGTC 


internal reference. The amplification condition was as 
follows: 95 °C 5 min; 95 °C 30 s, 60 °C 30 s, 72 °C 30 s, 
40 cycles; 72 °С 10 min. 2“4 relative quantification 
method was used. Tissues that showed the lowest 
expression were set as controls. Each sample had three 
replicates, and all experiments were repeated in triplicate. 


Cadmium (Cd?) stress test 

We determined the test concentrations following the 
method reported by literature (Kim et al, 2012). The 
standard Са?” concentration (<0.005 mg/L) according to 
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"Water quality standard for fishery" GB7471 (0.005 mg/L) 
as well as 1000 times the standard concentration (5mg/L) 
were respectively set as the control and experimental 
groups. The exposure method was static contact with 
CdCl, · 5Н;О infected water. Totally, 20 healthy mussels 
in each group were cultured in an aquarium with 20 L 
test solution. The water was changed by 50% every day 
and the concentration was maintained during the trial 
period. We randomly selected 3—4 mussels at exposure 
points of 0 h, 6 h, 12 h, 24 h, 48 h, 96 h, and from these 
total RNA from the hepatopancreas, kidneys, blood cells, 
and mantle was extracted. This experiment was repeated 
three times. 


Target gene expression analysis under Cd?" stress 

We performed SYBR  green-based real-time 
quantitative PCR assay (He et al, 2013; Peng et al, 2012a, 
b) to detect the expression changes of target gene in the 5 
mg/L concentration of cadmium stress conditions. The 
final relative expression result is the product of two 
relative expression values: the 2747 value of the 
experimental group (5mg/L) at different time points 
relative to Oh, and the 27407 value of the experimental 


group relative to the control group at the same time point. 


The mRNA expression difference of Hs-BTGI was 
analyzed using SPSS 16.0 (SPSS, inc. Chicago, USA). 
P«0.05 and P«0.01 were set as significant and highly 
significant, respectively. 


RESULTS 


Sequence analysis of Hs-BTG1 

The full length cDNA of Hs-BTG/ was found to be 
1402 bp and involves a complete coding sequence of 525 
bp, including an ORF encoding 174 amino acids, 306 bp 
5'-ОТК, 571 bp 3-UTR, with a typical AATAAA and 
poly A tail (Figure 1; sequence submitted to NCBI, 
accession number KF015272). BlastP analysis (http: // 
blast.ncbi.nlm.nih.gov/Blast.cgi) showed that the Hs- 
ВТСІ has one typical conserved N terminal domain 
shared by BTG Superfamily (i.e; BTG/TOB domain) 
which is characterized by two well-conserved boxes 
named A and B (Figure 1, Figure 2) (Fu et al, 2012; 
Sakaguchi et al, 2001; Winkler, 2010); SignalP-4.1 
prediction identified no signal peptide sequence. 


Homology and cluster analysis 
Amino acid sequences of BTGI from H. schlegelii, 
C. gigas, D. rerio, X. laevis, G. gallus, M. musculus, and 
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H. sapiens were compared using Clustal W (Figure 2), 
showing that Hs-BTG1 has the highest identity (50.57%) 
with BTGI from C. gigas. Hs-BTGI has lower identity 
with ВТСІ in vertebrates, with an identity of 35.67% 
with human ВТСІ. We next performed a homology 
modeling of ВТСІ from Н. schlegelii and C. gigas 
according to reported 3D structure of human BTGI 
protein (Yang et al, 2008; Winkler, 2010), and found 
each presents a similar 3D structure (Figure 2) with five 
a-helices and two antiparallel B-sheets. 

Cluster analysis (Figure 3) using Mega4.1 (Neighbor- 
Joining, NJ) showed that BTGs of vertebrates and 
invertebrates clustered into two branches. In invertebrates, 
the Hs-BTG1 was clustered firstly together with the С. 
gigas ВТСІ, and then BTGI of L. gigantea, L. Singor- 
iensis, A. florae, and C. formosanus along with the 
aforementioned two BTGI were clustered into one large 
branch. In vertebrates, another large branch was also 
formed mainly by some ВТС15 from fish to humans, as 
well as BTG2 from G gallus, M. musculus, and H. sapiens. 


Expression profile of Hs-BTG1 

Expression profile of Hs-BTG/ in different tissues 
was analyzed via real-time quantitative PCR (RT-qPCR) 
(Figure 4). Expression of Hs-BTG/ mRNA was detected 
in all 10 tissues—intestines, gills, gonads, mantle, ax 
foot, adductor muscle, blood cells, hepatopancreas, 
kidney, and heart—but was most highly expressed in the 
hepatopancreas and then the gills, with while blood cells 
exhibited the lowest expression of Hs-BTG1. 


Expression change under Cadmium stress 

Expression changes of Hs-BTG1 in hepatopancreas, 
kidneys, blood cells, and mantle under Cadmium (Cd?) 
exposure for Oh, 6 h, 12 h, 24 h, 48 h, 96 h were shown 
in Figure 5. In the hepatopancreas, Hs-BTG/ expression 
significantly increased at 6 h, decreased at 12 h, then 
increased at 24 h, reaching a peak at 48 h (2.03-fold of 
the standard control group) and then falling back at 96h. 
The overall expression pattern in blood cells was quite 
similar, with the only difference being a continued 
decrease at 24 h, though the expression level also peaked 
at 48h (2.5-fold of the standard control group). In kidney 
tissue, the Hs-BTG/ transcription dropped until an 
increase at 24 h (2.2-fold of the standard control group) 
before slowly returning to normal levels between 48 h-96 
h. Similarly, Hs-BTG/ expression level in mantle also 
reached the peak at 24 h (2.07-fold of the standard control 
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1 GATTTACATAGCGTAGGCT GCCAGT GGAGGTTTCACAGCGAACA 
45 AAGAGTGAACAGTCTGGTGGGAAGTATTTGGGTGCTCTCGAGGAAAAGTCATTTGCGATCAGCTACATTTCGAAA 
120 ATTTTAAGCAAAATATAATGTTCAACATCACCGTCTATACAAAAAGTACTCTATAAAGCTGGAAAAACAACATTG 
195 AGTGGTTCGTGAAGTGTTATCTTTGTGCAGTTGGTAGCCATCTTGTTTGGGGCTTGTTGTCAGCGGAGAAAAAAC 





























270 TGTCGAAGAGATAATATATGATTAAAAACTGCTCGTTATGAAGAAAGAAGTTAAAAGTGCTGTGGATTTTTTAAC 
1 ME GU ЕЕ Ж ҮК Ж С ТАР ҮКЕ ЕТ 
345 GAACATTCTCAGATCCAGGAATGTAGATCGTGTACAAACACATTTCTTCAATACAAAGCTGCAAAATCTGCTTCA 
14 МЕТЕ КИ ЕРУ ову Ке ШҮ (а) бб НИЕТ ФЕ КТ С МЕ ЈЕ ЈБ ја 
Вох А 
420 TAGTCGQTACGAAGACCATTGGTTTCCTGACAAACCATGCAAAGGTAGTGGTTATAGATGTATCCGATTCAACCA 
39 S ROC e i јр И ЈУ је 9) Је Је E Қона GNO MA (rT О (6 M R ЕШ БЕП 
495 CTCGTTAGACCCTATAGTCCGGGAAGCT GGCTTATCGTGTGGAATTCAAGAGACTACGATGAAATCCTTITTCCC 
64 © Jb i le mm Wie ie A Jb S Је © wa) РЕТ ТИЕ Је S dale IP 
Box B 

570 CAGGGAGTTGACGATGTGGGTAGATCCCCACGAAGTTTCCTACAGGATAGGAGAQAACGGCAGTGTTGGTGTGTT 
89 cry М G S V GV L 
645 GTTTTCTAGTGAAAACGATGCAGATCAAGAAAATACCTCACAATCAGAGGATGCCAGTACCGGAGAAACATCGGA 
114 FS SE N DADQE N T35Q S ED AS TGE T S D 
720 CCTTAACAACTCCTACAACTTAACCTGTAAGGAACAATTTTTGAGTGCTTTGCCAGGGCTTTCAATGAGCGGGGT 
139 L NN SY NL T C K EQFL S AL P GL SM $86 V 
795 AAACCTTAATAAATTTGCCACTTTGGTTTCCAGTTGAGGAAATACGAGACTGTTTCTATCACCGTCTCATTCATT 
164 NL NK F A TL ИУ 5 5 ж 

870 ATATCATTCTATTTCGCACGAGGACCACGACACTGGCTTTACTTCCATGGACTGCGATTTTATGTATTTCAGTAA 
945 TTGTTCGACATTTCGTCTTGCATGGTGTCGTAGATAAAATGGTAATGATGTGCTTGTCATTGTGATAGAATAATT 
1020 GATCATGATGAATGTAGGATATTTTAATATCATTAATATCATGAATGATCAATTGATATTTTTCGTGAATTCACA 
1095 TTCATTCAATTCGAGATTTAATTTAATGGCATTATTTACAATGTCGGGAGCTCAGGATGAGAAGTATTTTATTGC 
1170 ATTTTTAATTTTCACACTGCCTTACAAGCAGCGATGAATGTCTTTCAGACTAAATTTTATGCACATCCTGATCTC 
1245 AAGTTATTTTATGTGTAAATTTGTGTACATTGCTTGTACAGATGACCAAATGTTGAAGCTGAAAGCTATTTTGTC 
1320 TCTGTACAAGCAGATATATAGTGTTTCTACAAAAAACAAAAAAATAATGGAATAAAAATTACCTACTTAAAAAAA 
1395 АААААААА 

1 25 50 75 100 125 150 174 
Query seq. Та ата а шышы cita аи e ene c Ld Ra a a c a 
Specific hits BTG 


Superfamilies BTG uperfamil; 


Figure 1 cDNA and deduced amino acid sequence of BTGI in Hyriopsis schlegelii 
Nucleotide and deduced amino acid sequence were numbered on the left. The BTG/TOB function domain was shaded and the two highly conserved domains 
(BTG Boxes A and B) were indicated with a rounded rectangle. The stop codon TGA and the polyadenylation signal AATAAA were tagged with an asterisk and 
underline, respectively. The lower image depicts the BlastP result for Hs-BTG1. 


group), but decreased at 48 h, and then showed a slight 
increase at 96 h. 


DISCUSSION 

To date, there are 6 identified members of the 
BTG/TOB family in vertebrates, named for its BTG tag 
sequence in the N terminal. The BTG region contains 
two highly conserved homologous short motifs, namely 
A box (YKHHWFPDRPNKGSGYRCIRIN) and B box 
(LPSELTLWVDPFEVSYRIGE). These two boxes are 
separated by 20-25 nonconservative amino acids 
(Guehenneux et al, 1997; Sakaguchi et al, 2001; Fu et al, 
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2012). The A box seems to have an anti-proliferation 
function, while the B box appears act as the binding sites 
of its target molecule (Yang et al, 2008). Among mam- 
mals, the BTG/TOB family proteins can be classified 
into 3 subfamilies based on different sequence length of 
the C terminal: ВТСІ and BTG2 as a subgroup, BTG3 
and BTG4 as another subfamily, while ТОВ1 and TOB2 
as the third subgroup (Matsuda et al, 2001; Mauxion et al, 
2009; Winkler, 2010), with both ВТСІ and BTG2 having 
a shorter C terminal and an extra 8—10 amino acids in the 
N terminal compared with other two subgroups (Winkler, 
2010). 
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50.57% 


35.67% 
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Figure 2 Multiple alignment of BTG1 between Hyriopsis schlegelii with other species 


(A) Alignment of the deduced amino acid sequence of Н. schlegelii BTG1 with the corresponding sequences from other species; (B) the tertiary (3D) structures 
of Homo sapiens ВТСІ (а), Crassostrea gigas ВТСІ (b), and Н. schlegelii BTG1 (c) modeled in SWISS-MODEL. The highly conserved A Box and B Box 
regions are marked іп (A) and a schematic representation of the a-helical (yellow) and B-sheet (purple) structural elements are present in (B). a. H. sapiens: 
МР 001722.1 BTGI, b. C. gigas: EKC27510.1 ВТСІ, с. H. schlegelii BTG1: AGT79958.1; C. gigas: EKC27510.1; Danio rerio: МР 956314.1; Xenopus laevis: 


NP 001080825; Gallus gallus: CAA45507.1; Mus musculus: NP. 031595.1; H. sapiens: NP. 001722.1. “...” 


default at the corresponding sites. 


Though ВТС1 is comparatively well characterized 
in mammals and other vertebrates, little has been done 
on freshwater shellfish. Here, for the first time, we 
identified a homologous gene (Hs-BTG/) of the 
BTG/TOB family in H. schlegelii. The motif structure 
of the deduced amino acid sequences also contains 
highly conserved A and B box. BlastP analysis showed 
that Hs-BTGI has the highest identity (50.57%) with C. 
gigas ВТСІ. Homology modeling further showed that 
both H. schlegelii and C. gigas BTG1 present similar 
3D structure with the reported human ВТСІ protein 
(Yang et al, 2008; Winkler, 2010). Interestingly, both H. 
schlegelii and C. gigas BTG1 have a truncated (for 
8-10 amino acids) N terminal as compared with 
vertebrate BTGI, but the length of their C terminal is 
roughly consistent with the corresponding part of 
vertebrate BTG1. We also found that Н5-ВТСІ and C. 
gigas ВТСІ were clustered together, both belonging to 
the invertebrate BTG1 large branch, while vertebrate 
ВТС1 and BTG2 fell into another large branch. These 
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indicates the same amino acid; “---” indicates the 


results concur with both BlastP and traditional 


classification, so consequently the sequence we 
obtained was classified into the first subgroup in the 
BTG/TOB family and termed as the homologue (Hs- 
ВТСІ) to other species BTGI. 

It is worth noting that there are a very limited 
number of BTG/TOB members identified in inverteb- 
rates, making it difficult to accurately position the 
relationship between Нѕ-ВТСІ and other BTG/TOB 
members from the constructed cluster tree (Fu et al, 
2012), which is supported by the low Bootstrap value in 
Additionally, ВТСІ from 
invertebrates (such as C. gigas and H. schlegelii) showed 


relatively low identity compared with vertebrate BTGI, 


some branches. while 


they have conserved motif and similar 3D structure 
shared by the entire BTG family, implying that BTG1 of 
lower vertebrate animals may retain some conservative 
function and also possess different functions as 
compared to their vertebrate homologs (Matsuda et al, 
2001). 
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Figure 3 Phylogenetic tree of BTG1 from Hyriopsis schlegelii and other animals 
Mus musculus ВТСІ: NP_031595.1; Homo sapiens BTG1: NP_001722.1; Gallus gallus BTG1: NP_990681.1; Danio rerio BTG1: NP_956314.1; Oryzias 
latipes BTG1: ХР 004083105.1; Takifugu rubripes BTG1: ХР 003972839; Xenopus laevis BTG1: NP_001080825; С gallus BTG2: ХР 418053; M. musculus 
BTG2: NP_031595.1; H. sapiens BTG2: NP_006754.1; Coptotermes formosanus BTG1: AGM32529; Apis florae ВТС1: ХР 003697190; Lycosa singoriensis 


ВТС1: ABX75488; Гоша gigantea BTG1: Е5088805; H. schlegelii BTG1: AGT79958.1; C. gigas BTG1: ЕКС27510.1. The number in the branch is the 
bootstrap value. 
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Figure 4 Expression level of Hs-BTG/ mRNA in different tissues 
The vertebrate BTG family was previously reported keeping cell cycle arrest in the GO phase. In vitro 
experiments further demonstrated that overexpression of 
ВТСІ can inhibit cell proliferation and promote 


to be involved in a variety of biological processes, e.g., 
cell proliferation and differentiation. For example, ВТС1 


expression is maximal in the G0/G1 phases of cell cycle 
and down-regulated throughout the G1 phase (Rouault et 
al, 1992; Matsuda et al, 2001), suggesting that the 
expression of BTGI may inhibit cell proliferation by 
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apoptosis (Rodier et al, 2001; Hata et al, 2007; Lee et al, 
2003; Matsuda et al, 1992; Rouault et al, 1992). Several 
studies have shown that ВТС1 is also involved in cell 
differentiation and early embryonic development, for 
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Figure 5 Fold inductions of Hs-BTG/ mRNA levels in hepatopancreas (a), kidney (b), haemocytes (c) and mantle (d) after exposure 
to waterborne Cd exposure (5 mg/L) 


Data represent group mean fold induction+SD (n=3 in each time). Asterisks indicate mean fold induction, which are significantly (P«0.05) different from 


control values (0.005 mg/L). 


example in regulating angiogenesis (Iwai et al, 2004), 
promoting muscle cell differentiation (Busson et al, 
2005), and the gastrulation movement in early Xenopus 
embryo (Wessely et al, 2005). 

Several lines of evidence suggest that the biological 
function of BTG is correlated to its expression, such as in 
adult pigs and sheep where constitutive expression 
of BTG1 could be detected in several different tissues 
(Feng et al, 2011; Zhang et al, 2009). Another study 
found that the highest expression of human 8702 was in 
the renal proximal tubule, alveolar epithelial cells and 
prostate basal cell layer (Melamed et al, 2002). During 
the early embryonic development of Xenopus, abnormal 
XBTGI expression may cause embryonic gastrulation 
failure (Saka et al, 2000), while overexpression of a 
homologous gene x-BTG-x may induce a double axis of 
the embryo (Wessely et al, 2005). Moreover, expression 
of BTG-b was detectable in embryonic induction region, 
as well as the forebrain, hindbrain and paraxial meso- 
derm in zebrafish (Sakaguchi et al, 2001). Meanwhile, 
the highest level of 8701 expression in adult grass carp 
was found in the liver (Fu et al, 2012), consistent with 
our present results, suggesting that a maintained high 
expression of Hs-BTG/ may play an important role for 
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implementation of normal physiological function of the 
hepatopancreas. 

As we noted earlier, gene expression analysis under 
various environmental stress factors is a common method 
of studying candidate gene's corresponding function(s) 
(Duan et al, 2013; Wang et al, 2011). In the present study, 
the highest Hs-BTG/ transcription level was in the 
hepatopancreas, with the expression being significantly 
up-regulated when placed under conditions of heavy 
metal cadmium stress, and expression in other tissues 
(e.g., blood cells, mantle, and kidney) also increased 
under stress, though to a lesser extent. Since ВТС1 was 
previously reported to have anti- proliferation and 
apoptosis effect in higher vertebrates (Corjay et al, 1998; 
Lee et al, 2003; Rouault et al, 1992), we speculate that 
Н5-ВТСІ may also be involved in a similar function 
under cadmium stress. One intriguing observation 
garnered from this study was that Hs-BTG/ in kidney 
tissues was notably up-regulated, but then fell back. A 
potential explanation may be that the kidney serves a 
temporary storage of accumulated environmental 
cadmium. If so, then inhibition of Hs-BTG/ expression 
may contribute to proliferation of the renal cells, and 
then the increased cells would be able to accumulate 
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greater amounts of cadmium, while conversely a rapid 
increase of Нѕ-ВТСІ caused by an exceeding threshold 
level for renal cells tolerance against cadmium may 
promote cell apoptosis. Though an intriguing possibili- 
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Cloning the sterol carrier protein 2 genes of Japanese toad (Bufo 
japonicus formosus) and Chinese toad (Bufo gargarizans) and its tissue 
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Abstract: In this study, to clarify the bioactive polypeptides included in the skins and secretions of Bufo, we screened the 
Japanese toad (Bufo japonicus formosus) skin cDNA library by colony polymerase chain reaction (PCR), and obtained a 
transcript of 1 075 bp consisting of 137 bp 5' untranslated region (UTR), 515 bp 3' UTR and a 423 bp open reading frame 
(ORF) encoding a polypeptide of 140 amino acid residues (GenBank accession number: KF359945). Homolog analysis 
showed a 70%-96% homology with sterol carrier protein-2 (SCP-2) present in other animals, which is implicated in lipid 
metabolism of other organisms. The gene SCP-2 of Chinese toad (B. gargarizans) was cloned from a first strand cDNA of 
Bufo skin (GenBank accession number: KF381341) via PCR, whose encoding polypeptide has only one amino acid difference 
from that of Japanese toad. Tissue distribution analysis showed that SCP-2 expressed in all organs tested, though in the liver 
and spleen it manifested lower expression than in other organs. These findings might indicate SCP-2 being one of the active 
ingredients in toad skin. These findings may in turn have implications for further drug development from traditional Chinese 


medicine sources. 


Keywords: Bufo gargarizans; Bufo japonicus formosus; SCP-2 cDNA cloning; Tissue expression 


Amphibian skin and their secretions have been 
shown to contain large amount of biologically active 
compounds, suggesting a new potential source for drug 
discovery (Clark, 1997; Lai et al, 2002a, 2004; Novković 
et al, 2012; Rash et al, 2011; Zhao et al, 2014). The 
potential for novel therapeutics derived from these 
tissues is not unexpected; skin from the Bufo toad 
(Chan'pi), its cortex (Chan'yi) and secretions (Chan'su) 
have long been important components included in many 
prescriptions of traditional Chinese medicine (TCM) 
used in clinical treatments of several diseases, especially 
tumor control (Efferth et al, 2009; Liu et al, 2009; Tong, 
2011; Xin et al, 2012). Previous reports showed that 
cinobufocini injection (water soluble extracts of toad 
skin) possessed excellent anti-tumor curative effects (Qi 
et al, 2010, 2011; Zhou et al, 2009), likely due to their 
unique polypeptides (Wu et al, 2012). Despite these 


Science Press 


promising findings, little has been done to advance the 
use of the active compounds in toad skin, especially in 
China. 

China itself is in rich of amphibian species and has 
a long history of developing traditional medicine from 
unorthodox sources. Similarly, China is up and coming 
player in the pharmaceuticals and drug development. The 
key challenge in leveraging these two advantages is first 
gaining a clearer understanding of the underlying genetic 
and molecular mechanisms in many traditionally used 
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treatments, and second, finding better ways to utilize the 
available resources more efficiently (Lai et al, 2002b). 
We previously sought to elucidate the active polypeptide 
components included in toad-skin and its related 
materials by first screening the skin plasmid cDNA libr- 
ary of Japanese toad (Bufo japonicus formosus) via col- 
ony polymerase chain reaction (PCR) (Yuan et al, 2013; 
Zhang et al, 2013; Zhuge et al, 2013). As part of our 
other research efforts, we have also begun cDNA cloning 
from Chinese toad (B. gargarizans) skin first strand 
cDNA (Hu et al, 2013). During these processes, we were 
able to isolate the cDNAs encoding sterol carrier protein- 
2 (SCP-2) from both Japanese toad and Chinese toad. 

SCP-2 was initially segregated and purified from 
mouse (Mus musculus) liver tissue (Osumi et al, 1980), 
which is distributed in peroxidase, mitochondria, 
endoplasmic reticulum and cytoplasm, and functions as 
intracellular transporters such as cholesterol, lecithin, 
fatty acid, ester acyl CoA etc (Kriska et al, 2010; 
Schroeder et al, 2000, 2007). Here we cloned SCP-2 
from two Bufo species and conducted RT-PCR analysis 
on the Chinese toad, and found that SCP-2 is likely one 
of the effective polypeptides included in toad-skin origin 
materials. These findings confirm our earlier observation 
that toad skin (and potentially other tissues) may be 
viable targets for future development of drug treatments 
and novel therapeutics. 


MATERIALS AND METHODS 


Experimental materials and reagents 

The proprietary Japanese toad skin plasmid cDNA 
library held by the Japan Advanced Industrial Science 
and Technology (AIST, Tsukuba, Japan) was authorized 
for use by Zhejiang Agricultural and Forestry University 
(ZAFU) for research as part of a Material Transfer 
Agreement. Concerning this library, pSD64TR (3 250 bp) 
has been used as a vector, and EcoR I and Xho I as 
cloning sites. The upstream primer of the vector is SP6 
(5'-ATTTAGGTGACACTATAGAA-3’) and the downs- 
tream one is S.D.A. (5'-TTATGTAGCTTAGAGACTC- 
3^, respectively. The cDNA length ranged from 
500—2000 base pairs (bp). For further testing, Chinese 
toad individuals were obtained from the East lake 
Campus of ZAFU, and then ice anaesthetized prior to 
dissection, wherein the organs were removed and cut into 
small pieces before being frozen in liquid nitrogen. The 
resulting samples were kept in a —70 ?C refrigerator prior 
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to total RNA extraction. 

The RNA extraction kit was purchased from 
Shanghai Bocai Biotechnology Company; Quantscript 
RT kit, pGM-T vector and Escherichia coli competent 
cells (DH5a) from Tiangen Biotechnology Limited 
Company; Тад PCR kit Кот TaKaRa Biotechnology 
Limited Company; and primer synthesis and DNA 
sequencing were commissioned by Shanghai Sang'ni 
Biotechnology Company. 


B. japonicus formosus SCP-2 screening 

Japanese toad cDNA screening was performed as 
described previously (Yuan et al, 2013). In brief, 
Japanese toad skin plasmid cDNA library was transfo- 
rmed into Е. coli (DH5o), and colony PCR was perfor- 
med using colony suspension as templates, and SP6 and 
XhoTT (5-АСАТСТСТССАСТТТТТТТТТТТТ-3, a 
self-designed primer complementary with the area 
compassing the connection point of cDNA polyA tail and 
the downstream cloning site of Xho I) as primers. 
Following this process, the recombinant plasmids were 
collected and double enzyme digested with EcoR I and 
Xho 1 to further confirm positivity, and then sent for 
sequencing with vector primers SP6 and S.D.A. 


B. gargarizans SCP-2 cloning 

For SCP-2 cloning from Chinese toad, total RNA 
was extracted from its dorsal skin, and a first strand 
cDNA synthesized based on the manufacturer's protocols. 
Meanwhile, based on the Japanese toad SCP-2 sequence, 
an upstream primer (SCP-2-S: 5-CGTGGTCGTTACG 
TTATACAAG-3’) and a downstream primer (SCP-2-R: 
5'-GAAATTAGTGGCTTTTATTAAGTG-3’) were desig- 
ned for use in RT-PCR. The PCR product was ligated 
into a pGM-T vector and then sequenced with vector 
upstream primer T7 and downstream primer SP6. 


Sequence Analysis 

DNAstar/EditSeq was used to find the open reading 
frame (ORF) and deduce their encoding protein amino 
acid sequence. Potential phosphorylation sites were 
predicted via Net Phos 2.0. A further 17 SCP-2 protein 
sequences from other animals were downloaded by 
NCBI blast program (http: //blast.ncbi.nlm.nih.gov/Blast. 
cgi) and aligned with DNAstar/MegAlign. Phylogenetic 
tree based on SCP-2 amino acid sequences was 
constructed using the neighbor-joining method by 
MEGAS.1 (bootstrap with 1 000 replications). 
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SCP-2 tissue expression analysis by RT-PCR 

Total RNA samples were extracted from different 
Chinese toad organs including brain, heart, lung, liver, 
spleen, kidney, stomach, intestines, fallopian tube and 
skin, and their first strand cDNA were synthesized as 
mentioned above. Samples of the different tissue were 
analyzed via PCR for SCP-2 expression using a refer- 
ence gene of ф-асап (upstream primer: 5-ТТСАСАС 
CTTCAACACC-3'; downstream primer: 5'-CTTGATGT 
CACGCACAA-3’). 


RESULTS 


SCP-2 screening and sequence analysis of B. japo- 
nicus formosus 

Restriction enzyme digestion with EcoR I and Xho 
I showed one recombinant plasmid screened from 
Japanese toad skin plasmid cDNA library had a cDNA 
insert of about 1000 bp, which was later confirmed by 
sequencing analysis (Figure 1). The transcript is 1 075 
bp consisting of 137 bp 5'UTR (untranslated region), 
515 bp 3' UTR and a 423 bp ORF encoding a 
polypeptide of 140 amino acid residues, which showed 
high homology with sterol carrier protein-2 (SCP-2) 
found in other animals. The clone we screened in the 
present study (named В. japonicus formosus SCP-2), 
has been deposited into GenBank (accession number: 
KF359945). 


SCP-2 cloning and sequence analysis of B. gargarizans 
From Chinese toad skin first strand cDNA, a 920 bp 


transcript was obtained consisting of 22 bp 5' UTR, 475 bp 
3' UTR and 423 bp ORF encoding a polypeptide consisting 
of 140 amino acid residues (Figure 2). The only difference 
between two Bufo SCP-2 proteins 1s that Thr130 in Chinese 
toad was substituted by Ser in the Japanese toad. This clone 
was designated as B. gargarizans SCP-2 and deposited into 
GenBank (accession number: KF381341). 


Phosphorylation site prediction of Bufo SCP-2 

From the analysis of phosphorylation site prediction, 
10 potential sites (Serl0, Зег11, Зег19, Ser77, Ser79, 
Ser82, Ser91, Ser93, Thr66 and Tyr44) were found in 
both Bufo species (Table 1), suggesting that SCP-2 
expression might be regulated by the upstream factors. 


Homology analysis of SCP-2 amino acids 

Phylogenetic analysis showed that two Bufo species 
had a homology as high as 96% with Xenopus (Silurana) 
tropicalis and 92% with Xenopus laevis, with lower 
homology ranging from 70% to 91% among 15 other 
animals (Figure 3). The phylogenetic tree we constructed 
showed 9 different mammals (Homo sapiens, Bos taurus, 
Rattus norvegicus, Mus musculus, Mesocricetus auratus, 
Capra hircus, Camelus ferus, Ictidomys tridecemlineatus 
and Sus scrofa) gathered in a branch, 5 different birds 
(Gallus gallus, Falco peregrinus, Falco cherrug, 
Pseudopodoces humilis and Melopsittacus undulatus) in 
a branch, 1 reptile (Anolis carolinensis) in another 
branch, 2 frogs (X. laevis and X. Silurana tropicalis) in a 
branch, 2 toads (B. gargarizans and B. japonicus 
formosus) in a branch, and 2 fish in a branch (Figure 4), 
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CCACAGGCCTGGC TCAGTGTGCAGAGC T CTGC TGGC AGCTCCGAGGTCAGGC TGGGAAAAGACAGGTTCCGGGA GCCAAGC TCGCTC TGC 
AGC ACAAC ATCGGACTAGGGGGGGCCGT GGTCG TTACG TTATAC AA GCT TTCCGGACGC AGCAGCC AGA TCTTCCCGTATCCAGC 
HGFPDAAARGSSRIQ 


TGAATCCGACCAGCGCCGAAGACGGAT T TAAAGCCCAG TTTGTA TTCAAGGAAA TTGAGAAGAAAT TGAAAGAGGAAGGGGACC AGT ATG 
L NP TS AED GEK AQFVEFKE I EKKLKEE GE Q Y 


TTAAGAAGATTGGAGGAGTCT TTGCCTT TAAAGTGAAGGAT GGACC TGG TGGAA AAGAGGCAAC TTGGG 16611 GATGT GAAGAACGGC A 
VKKIGGVEF AF KVKDGPGGKEAIK.wWVVDVKNG 


AMGGCTCTGTGTCCTTCGACTCCGATAAGAAMAGCGGAC TGTACGATCTC AATGTCCGAC TC TGACC TAT TGGCT CTGAT GACCGGCC AGA 
ко зу SFEDSDKKADCTIISMSDSDLLALHMTGO 


TCAATCCACAGACCGCTTTCTTCCAGGGCAAGC TGAAAGTCACT GGAAA TATGGGTC TGGCCAT GAAGT ТОСАСАСССТССАСС TGC AGC 
I NPOüTAFFQOGKLIKVvrIGNHGLAMKLQOSLOQILGQ 


CTGTGAAAGCCAAGCTCT GAAGAGATCGTGCGT GCGAT ACATGAGAAAT GCCAGCAGTACCGCAGATGCCAACCCT TCT GATGC CCAGAT 
PVKAKL . 


GICATCGA TGCTGGCGTC TAGAGGACCAATCAC AACTGCTGCTT TCCC T ACCCAAGAAT AAAGAATCAT GTCCCCTTTATTTCTAAGCAG 
GIGGCAGT GCATC TGCGC ACGTGTACATGTATATAGGTGTACGAATCCGGAA TC CTGGC CTACGATGCC AGTCTGTTTAGAAAT TAAGGC 
GTTTGTAGGGATGGGGC T TCTAACGGTGTGACT TAATCTGATGTCTTCTGTCTG TC TTC TCAAC AC TAA TGGAT TTAAATAATT TCGGAG 
GCGGATCT ATAGC GCAGC AGTGAT TTAATCCCTGTAATCTTGGGAT AGGAGGCAGGTCT CC TGC CCGGC GCC TGAACGC CACGGAT GC AG 
AATTCTCTAACCTGTTAATAGTC AGAAT AACAC TTAATAAAAGCCACTAATTTC GAAAA AAAAA AAAAA AAAAAAAAAAAAAAAA 








Figure 1 SCP-2 cDNA and its deduced amino acid sequence of Bufo japonicus formosus 
Start and stop codons were enclosed by line box. —: Polyadenylation signal; =: Poly(A) tail. 
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CGTGGTCG TTACGTTATACAAGA TGEGC TTTCCGGACGCAGCAGCCAGATCTTC CCGTATCCAGCTGAATCCAACCAGCGCCGA AGACGG 90 
I— MH GFPDAAARSSRIQLNPTSAED 6 


ATTTAAAGCCCAGTTTGTATTCAAGGAAATCGAGAAGAAAT TGAAAGAGGAAGGGGAGC AGTACGTTAAGAAGATTGGAGGAGTCTTTGC 180 
FK AQF VF K Е Т ЭЕ КОК LK EE G@EaQyYV¥K К IG б VF А 


CTT TAAAG TGAAGGATGG ACCTGGTGGAAAAGAGGCAA CT TGGG TGGTT GATGT GAAGA ACGGCAAAGGCTCGGTGTCCTTCGACTCCGA 270 
FKVKDGPGGKEAIWwWVVDVvVk NGKGS3VvVsFDSsD 


TAAGAAAGCGGAC TGTACGATCT CAATGTCCGACTCTGACCTAT TGGCT CTGATGACCG GCCAGATCAATCCACAGACCGCTTTCTTCCA 360 
ККАПСТІ 5 М5П5ПІІЛІ МН ТЕСІМНМРОТАҒРҒҒ 4 


GGG CAAACTGAAAGT CAC TGGAAAT ATGGGTCT GGCCA TGAAGT TGCAG ACT СТ CCAGC TGCAGCCTGT GAAAGCCAAG CTCTGAAGAG А 450 
а КІ КУ TGNMGL AMKL OT L OL OP V KAKI < 


TCGTGCGTGCGATACATG AGAAAAGCCAGCAGT AGAGGAGATGGGAAGGGT ICT GATGCCCAGATGTCATCGATGCTGGCATCTACAGGA 540 
CCAATCACAACTGCTGCT TTCCC TACCC AAGAA TAAAGAAT CAT GTCCCCTTTATTTCT AAGCAGGTGG CAGTGCATCTGCECACGTGTA 630 
CATGTATATAGGTGTACG AAT CCGGAGT CTACG ATGCCAGTCTG TTTAGAAATT AAGGCGTTTGTAGGG ATGGGGCTTCTAACCGTGTGA 720 
CTTAATCTGGTGTCTTCTGTCTGTCTTCTCAACACTAATGGATT AAAATGATTT CAGAGGCGGA TCTAT AGCGCAGCAGTGATTTAATCC 810 
CTGTAATCTTGGGATAGG AGGCAGGTCT CCTGCCCGGCGCC TGA ACGCCACGGA TGCAG AATTCTCTAACCTGT TAATAGTCAGAATAAC 900 
ІШ TAATAAAAGCCACTAAT 920 


Figure 2 SCP-2 cDNA and its deduced amino acid sequence of Bufo gargarizans 
Start and stop codons were enclosed by line box. 1: Upstream primer, SCP-2-S; 2: Downstream primer, SCP-2-R; —: polyadenylation signal. 
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Figure 3 Bufo SCP-2 amino acid sequence alignment with other animal species 


which is consistent with the traditional animal tax- 
onomy. 


SCP-2 expression in different organs of B. gargarizans 

RT-PCR detection of B. gargarizans SCP-2 in the 
brain, heart, lung, liver, spleen, kidney, stomach, 
intestines, fallopian tube and dorsal skin (Figure 5) 
showed that SCP-2 was expressed in all tested organs, 
though the expression was lower in both liver and spleen 
as compared with other organs. 
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DISCUSSION 


Previous studies found that SCP-2 is involved in 
adjusting concentrations of cholesterol inside and outside 
cell membranes via activation of cholesterol hydrase, and 
participating in regulating cholesterol transport in the cell 
culture system (Kriska et al, 2010; Schroeder et al, 2000, 
2007). Many human diseases—diabetes, arteriosclerosis, 
Zellweger, NPC disease and gallstones—have similarly 
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Figure 4 Phylogenetic tree of amino acid of SCP-2 between two Bufo species and 17 other species 
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Figure 5 RT-PCR detection of SCP-2 expression in different 
organs of Bufo gargarizans 
M: DNA ladder; 1: brain; 2: lung; 3: heart; 4: liver; 5: spleen; 6: kidney; 7: 


stomach; 8: intestines; 9: fallopian tube; 10: dorsal skin. 


Table 1 Bufo SCP-2 phosphorylation site prediction 





Position Context Score Prediction 
10 AAARSSRIQ 0.552 *S* 
11 AARSSRIQL 0.598 *S* 
19 LNPTSAEDG 0.995 *S* 
77 NGKGSVSFD 0.943 *S* 
79 KGSVSFDSD 0.902 *S* 
82 VSFDSDKKA 0.994 *S* 
91 DCTISMSDS 0.992 *S* 
93 TISMSDSDL 0.955 *S* 
66 GKEATWVVD 0.577 *т* 
44 EGEQYVKKI 0.926 жүж 


been found to have associations with abnormal expre- 
ssion of SCP-2 (Castelli, 1984; McLean et al, 1996). For 
example, NPC (Niemann-Picktype C) disease seems to 
be caused by NPC-peak C protein mutations, and a 
corresponding drop in liver SCP-2 expression was 
indicated (Schroeder et al, 2007). In diabetic mice 
models induced by streptozotocin, the level of SCP-2 in 
the liver was reduced 60%-90%, alongside a 60% 
reduction in ovarian SCP-2 (McLean et al, 1996). SCP-2 
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expression is also related to the formation of cholesterol 
calculus (Cui et al, 2011). 

In the present study, our successful cloning of sterol 
carrier protein-2 (SCP-2) genes from both B. gargarizans 
and B. japonicus formosus (Figures 1, 2) indicated that 
SCP-2 expressed in toad skin as well as in other organs 
(Figure 5) has a high homology with that in other 
animals (Figure 3, 4). Due to the function of SCP-2 in 
adjusting lipid metabolism in numerous animal species 
and many SCP-2 expression-related disorders among 
humans, it may be reasonable to predict that SCP-2 is 
one of potentially several important ingredients within 
toad skin. 

Our study extends the basic knowledge necessary to 
assess the potential for Bufo skin and other organs for 
potential drug development. Previously, numerous 
reports summarized several descriptions on the clinical 
efficacy of Chan’su, such as detoxification, analgesia, 
anti-inflammation, antidiarrheal, and antitumor, etc. (Liu 
et al, 2009; Xin et al, 2012). Curiously though, we found 
that SCP-2, as a potentially important element in toad 
skin, is largely involved in adjusting lipid metabolism. 
However, few reports have ever noted the use of Bufo 
skin in treating lipid metabolism related diseases, aside 
from an old description concerning Chan’su that pointed 
out that administration of the medicine could make the 
symptoms of stasis and 
“BenCaoHuiYan”. This 
intriguing, given that it appears in a Ming Dynasty era 


stagnation disappear in 


reference itself is rather 
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book written by Zhu-mo NI published in 1624. Presu- 
ming our finding survive replication and further targeted 
studies are conducted to answer some of the remaining 
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Abstract: The grass carp (Ctenopharyngodon idella) is one of the most important cultivated fish species in China. Mounting 


evidences suggests that microRNAs (miRNAs) may be key regulators of skeletal muscle among the grass carp, but the knowledge of 


the identity of myogenic miRNAs and role of miRNAs during skeletal muscle anabolic state remains limited. In the present study, we 


choose 8 miRNAs previously reported to act as muscle growth-related miRNAs for fasting-refeeding research. We investigated 


postprandial changes in the expression of 8 miRNAs following a single satiating meal in grass carp juveniles who had been fasting 


for one week and found that 7 miRNAs were sharply up-regulated within 1 or 3 h after refeeding, suggesting that they may be 


promising candidate miRNAs involved in a fast-response signaling system that regulates fish skeletal muscle growth. 


Keywords: MicroRNA; Grass carp; Fasting; Refeeding; Skeletal muscle 


The grass carp (Ctenopharyngodon idella) is one of 
the most important cultivated species in Chinese freshw- 
ater aquaculture. Generally, fast skeletal muscle compri- 
ses the largest tissue of the fish body (Zhang et al, 2009), 
and forms the main edible part of this species. Given the 
complicated development process as well as the impor- 
tance of these tissues, gaining a clearer understanding of 
the mechanism underlying muscle development may 
provide key data for both developmental biologists and 
researchers attempting to improve grass carp muscu- 
lature for aquaculture. 

Several lines of research have implicated growth 
factors, regulatory proteins, and transcription factors as 
key actors involved in the regulation and maintenance of 
skeletal muscle mass among fish (Nihei et al, 2006; 
Steinbacher et al, 2006; Chu et al, 2010). Recent studies 
found that alongside transcriptional factors involved in 
muscle proliferation and differentiation, a set of micro- 
RNAs (miRNAs) may also play important roles in skele- 
tal muscle development among vertebrate animals (Ge & 
Chen, 2011; Giller & Russell, 2010; Chu et al, 2013). 
MicroRNAs (miRNAs) are approximately 22-nt nonco- 
ding RNAs that act as negative regulators of gene expr- 
ession, either via inhibiting mRNA translation or prom- 
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oting mRNA degradation through base pairing to the 3' 
untranslated region (UTR) of target mRNAs (Xie et al, 
2005; Liu, 2008; Zhang & Wen, 2010). Furthermore, 
miRNAs regulate the expression of transcription factors 
and signaling mediators critical to both cardiac and 
skeletal muscle development and function (Callis & 
Wang, 2008; van Rooij et al, 2008). Studies using the 
mouse myogenic C2C12 cell line demonstrated miR-1 
and miR-133 are involved in myoblast proliferation and 
differentiation via regulation of the expression of 
HDAC4 and SRF, respectively (Chen et al, 2006). 
Meanwhile, miR-1, miR-133a and miR-206 (muscle- 
specific miRNAs), were also noted to be differentially 
expressed in Japanese flounder (Paralichthys olivaceus) 
during metamorphosis, as well as in Nile tilapia (Oreo- 
chromis niloticus) among several developmental stages. 





Received: 25 March 2014; Accepted: 17 July 2014 

Foundation items: This study was supported by the National Natural 
Science Foundation of China (31230076; 31340054), the Natural 
Science Foundation of Hunan province (14JJ2135) and the State Key 
Laboratory of Freshwater Ecology and Biotechnology (2012FB01) 
“Corresponding authors, E-mails: jzhang@ccsu.cn; chentao_114@ 
163.com 


Volume 35 Issue 3 


MicroRNA signature in response to nutrient restriction and re-feeding in fast skeletal muscle of grass carp (Ctenopharyngodon idella) 405 


Taken together, these findings imply that miRNAs likely 
play an important role in regulating muscle development 
(Fu et al, 2011; Yan et al, 2012a). Further evidences 
suggest that miRNAs act as key regulators of 
myogenesis, but unfortunately characterizing the identity 
of myogenic miRNAs and delineating the role of 
miRNAs during skeletal muscle anabolic state remains 
unclear. 

The maintenance of skeletal muscle mass is a 
complex and controlled process, largely influenced by 
both nutritional and physiological states of different 
animals (Fuentes et al, 2012). Fasting-refeeding protoc- 
ols have been commonly used as models to investigate 
the regulation of muscle growth in fish species that 


experience the transition from catabolic to anabolic states. 


In the present study, we analyzed the expression of 8 
select miRNAs during skeletal muscle anabolic state 
using a fasting-refeeding experiment. The goal of our 
study was to better parse out the potential role of these 
miRNAs in skeletal muscle proliferation and different- 
iation. Тһе 8 miRNAs (miR-1a, miR-133a-3p, miR- 
133b-3p, miR-146, miR-181a-5p, miR-206, miR-214 
and miR-26a) were previously reported to act as 
muscle growth- related miRNAs (McCarthy & Esser, 
2007; Flynt et al, 2007; Kuang et al, 2009; Yan et al, 
2012b). In theory, these miRNAs and their target 
genes may comprise a coordinated regulating net to 
favor resumption of myogenesis as an early response 
to refeeding. 


MATERIALS AND METHODS 


Fasting-refeeding experiments and sampling 

All grass carp individuals were reared under stan- 
dard conditions at the Che Tian Jiang Reservoir in Loudi, 
Hunan, China. Two homogeneous groups of grass carp 
juveniles (average body weight 150 g, 90 days post- 
hatching (dph)) were reared respectively in two net cages 
(5 mx5 mx2 m) with fifty fish per tank. АП juveniles 
were fed under standard conditions for 3 weeks, after 
which the juveniles underwent fasting for 1 week, and 
were then fed a single meal, which was distributed to all 
individuals until they appeared to be visually satiated. At 
each time point from 0 h (before the recovery meal), and 
at 1, 3, 6, 12, 24, 48 and 96 h (hours after the single 
meal), six fish were sampled, wherein fast muscles were 
dissected from the dorsal myotome of individuals. The 
resulting samples were then snap-frozen in liquid 
nitrogen and stored at -80 °С until further processing. 
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Quantitative real-time PCR for the miRNAs 

Tissue samples were ground in liquid nitrogen, and 
total RNAs were extracted using TRIzol (Invitrogen, 
USA), and then treated with RNAse-free DNAse I 
(Promega, USA) in the presence of RNAse inhibitor 
(Sigma, China Branch) followed by ethanol precipitation. 
The obtained RNAs were polyadenylated by poly (A) 
polymerase, and then reverse transcribed with one step 
PrimeScript miRNA cDNA synthesis Kit (TaKaRa, Dal- 
ian, China) and a Universal Adaptor Primer (a poly (T) 
primer ligated with an adapter) for miRNA quantitative 
assays. 

The miRNA expression levels were quantified using 
real-time PCR with grass carp B-actin gene (GenBank 
No. DQ211096.1) as an internal control. The cDNA 
samples were used as templates for quantitative RT-PCR 
assays with SYBR Premix Ex Taq II (TaKaRa, Japan) 
and its amplification reaction was carried out on a Bio- 
Rad CFX96 system (USA). Each 2 uL cDNA template 
was added to a total volume of 25 uL reaction mix 
containing 12.5 uL SYBR Green mix, 1 uL of each 
miRNA or gene specific forward primer (as shown in 
Table 1, 10 umol/L) and ІШ, of universal downstream 
primer (Uni-miR qPCR Primer, 10 pmol/L, TaKaRa) or 
gene specific reverse primer, 8.5 uL nuclease-free water. 
The protocols used are as follows: (i) pre-denaturation at 
95 °C for 60 s; (ii) amplification and quantification, 
repeated 40 cycles of at 95 °C for 5 s and at 60 °C for 25 
S; (11) melting curve program (65-95 °C with heating rate 
of 0.1 °C/S and fluorescence measurement) (Zhou et al, 
2010). The relative expression ratio (R) of target miRNA 
was calculated by R=2“*“' (Livak & Schmittgen, 2001; 
Bustin et al, 2009), where Ct is the cycle threshold. The 
basic equation employed was: 


ААС! = (Ct. 
(Сї 


target gene | Снан ng gene Лана = 


target gene | Chrousckespi ng gene Jeon rol 


The miRNA expression levels were then analyzed 
by one-way ANOVA procedures and regression analysis 
of SPSS 17.0 (SPSS inc., Chicago, USA). Duncan’s 
multiple range tests were used to compare the control (0 
h before the recovery meal) and experimental (# hours 
differences were 


after refeeding) groups. The 


considered statistically significant when P<0.05. Data 





are shown as means+SE (n=6). Correlation of gene 
expression was analyzed by the Spearman rank order 
correlation test. Hierarchical clustering was performed 
using Cluster3. 
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Table 1 Primers used for miRNA detection 





Name Primer Sequence (5'—3”) Name Primer Sequence (5'—3”) 
miR-la-F TGGAATGTAAAGAAGTATGTAT miR-146-F CGTGAGAACTGAATTCCATAGATGG 
miR-133a-3p-F CGCGTTTGGTCCCCTTCA miR-133b-3p-F TTGGTCCCCTTCAACCAGCTA 
miR-206-F CGTGGAATGTAAGGAAGTGTGTGG miR-214-F ACAGCAGGCACAGACAGGCAG 
miR-181a-5p-F CGAACATTCAACGCTGTCGGT miR-26a-F CGTTCAAGTAATCCAGGATAGGCT 
p-actin-F GCCGTGACCTGACTGACTACCT B-actin-R CGCAAGACTCCATACCCAAGAAG 


Forward primers used for detection were detailed in the Methods; the reverse primer used for detection was universal downstream primer (Uni-miR qPCR 


Primer, 10 pmol/L, Takara). 


RESULTS 


Effect of fasting and refeeding on the expression of 
the miRNAs 

A significant up-regulation of each of the 8 
miRNAs was observed between 1-6 hours after the 
single meal (P<0.05) (Figure 1). MiR-1a was sharply up- 
regulated within 1 h after refeeding and peaked at 6 h 
(Figure 1A). The expression of miR-133a-3p, miR-133b- 
Зр, miR-146, miR-181a-5p, miR-206 and miR-214 were 
significantly increased at 3h (P«0.05), reached the 
maximal levels at 3, 6 or 48 h postprandial (Figure 1B- 
G). While miR-26a responded slowly to refeeding, it 
significantly increased at 6 h (P«0.05) and peaked at 48 
h (Figure 1F). However, miR-1a, miR-133b-3p and miR- 
181a-5p significantly decreased at 12 h. All 8 miRNAs 
returned to the initial baseline values at 96 h. 


Heat map summary of hierarchical clustering of 
miRNAs in skeletal muscle during fasting-refeeding 
periods 

Hierarchical clustering analysis of miRNAs in 
muscle was done according to the similarity in their 
expression across different postprandial times (0-96 h). 
Hierarchical clustering of the miRNAs throughout the 
trial showed three clades (Figure 2). The first clade 
including 2 pairs of closely linked miRNAs (miR-133b- 
3p and miR-133a-3p, miR-206 and miR-181a-5p) that 
clustered together. The second clade clustered miR-26a 
with miR-146 expression. 


DISCUSSION 


MicroRNAs (miRNAs) are noncoding RNA mole- 
cules that regulate the stability and/or the translational 
efficiency of target mRNAs, and several miRNAs have 
been found to be specifically expressed or highly 
enriched in skeletal muscle. The expression of muscle- 
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specific miR-1, miR-133, miR-206, while miR-208 is 
regulated by muscle transcriptional networks involving 
SRF, MyoD and MEF2. (Flynt & Lai, 2008; Latronico et 
al, 2007; Thum et al, 2008; Van Rooij et al, 2008; Callis 
et al, 2008). Interestingly, non-muscle-specific miRNAs, 
e.g., miR-26a and miR-181, also regulate skeletal muscle 
differentiation (Wong & Tellam, 2008; Naguibneva et al, 
2006). While a great deal of attention has been paid to 
miRNAs involved in control of muscle development, a 
recent study suggests that several other miRNAs, 
including miR-499, miR-208b and miR-23a, also play an 
important role in human muscle growth (Drummond et al, 
2009). Among fish species, Huang et al (2012) previo- 
usly detected differentially expressed miRNA between 2 
strains of Nile tilapia and identified miR-140, miR-192, 
miR-204, miR-218a, miR-218b, miR-301c, miR-460, 
miR-133, miR-152, miR-15a, miR-193a, miR-30b and 
miR-34 as being associated with body growth in tilapia. 
Our present study focused on exploring the potential 
role(s) of miRNAs as a new layer of control in the 
postprandial regulation of the muscle development 
among grass carp. Typically, nutrient availability is 
among the most important environmental variable 
altering muscle growth (Valente et al, 2012). As such, 
starvation and refeeding experiments have served as an 
effective model for studying the regulation of muscle 
growth in fish, including the Atlantic salmon (Salmo 
salar) (Bower et al, 2009), rainbow trout (Oncorhynchus 
mykiss) (Montserrat et al, 2007), and Atlantic halibut 
(Hippoglossus hippoglossus) (Hagen et al, 2009). 
Similarly, in humans MiRNAs turned over quite rapidly 
(i.e. hours) in skeletal muscle following amino acid 
ingestion. However, little information is available 
regarding the early transcriptional changes of miRNA 
during the postprandial period, especially among fish. 
We therefore focused on exploring the postprandial 
regulation of growth-related miRNAs shortly after 
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Figure 1 Relative expression of miRNAs with significant up-regulation following refeeding after 1 week fasting B-actin expression 





was detected as the internal control. All values are presented as mean+SE, п-б. Different letters indicate significant differences 
between columns (P«0.05). 


feeding a single meal in grass carp. Our results showed elevated at 1 or 3 h after refeeding in the fast muscle of 
that miR-1a, miR-133a-3p, miR-133b-3p, miR-146, grass carp. Drummond et al (2009) previously found a 
miR-181a-5p, miR-206 and miR-214 were significantly rapid up-regulation of the miR-1, miR-208b, miR-23a 
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Figure 2 Heat map summary of hierarchical clustering of miRNAs in skeletal muscle during fasting-refeeding periods 


Blue and yellow respectively denote a decrease or increase. The absolute signal intensity ranged from -1.82 to +1.82, with corresponding color changes from 


green to red. 


and miR-499 following the amino acid ingestion in 
humans. These findings suggest that the identified 
miRNAs may be promising candidate miRNAs involved 
in a fast-response signaling system that regulates fish 
skeletal muscle growth. The other finding of significant 
decreased of miR-1a, miR-133b-3p and miR-181a-5p at 
12 h after single meal suggests there may be other 
signaling pathways regulated by the miRNAs that limit 
excessive regulation of muscle growth. 

A further finding of our study was that miR-206 and 
miR-181a-5p showed a dramatic and simultaneous up- 
regulation following feeding by a single meal. MiR-206 
is known to be a muscle-specific miRNA, with its role in 
muscle development having been verified in some 
animal models, including mice, rats and zebrafish 
(Anderson et al, 2006; Kim et al, 2006; Mishima et al, 
2009; Shan et al, 2009). MyoD acts as a transcriptional 
activator of the miR-206 pre-miRNA transcript, and 
subsequently the induced high levels of the mature miR- 
206 result in the down-regulation of specific target 
muscle growth-related genes (Rosenberg et al, 2006). 
MiR-181 is also thought to function partly through 
inhibition of Нох-А 11 expression, а known repressor of 
MyoD, which is required for new muscle growth 
(Naguibneva et al, 2006). Taken together, these different 
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Effects of surround suppression on response adaptation of V1 neurons 
to visual stimuli 


Peng LI, Cai-Hong JIN, San JIANG, Miao-Miao LI, Zi-Lu WANG, Hui ZHU, Cui-Yun CHEN, 
Tian-Miao НОА“ 


College of Life Sciences, Anhui Normal University, Wuhu 241000, China 


Abstract: The influence of intracortical inhibition on the response adaptation of visual cortical neurons remains in debate. To clarify 
this issue, in the present study the influence of surround suppression evoked through the local inhibitory interneurons on the 
adaptation effects of neurons in the primary visual cortex (V1) were observed. Moreover, the adaptations of V1 neurons to both the 
high-contrast visual stimuli presented in the classical receptive field (CRF) and to the costimulation presented in the CRF and the 
surrounding nonclassical receptive field (nCRF) were compared. The intensities of surround suppression were modulated with 
different sized grating stimuli. The results showed that the response adaptation of V1 neurons decreased significantly with the 
increase of surround suppression and this adaptation decrease was due to the reduction of the initial response of V1 neurons to visual 
stimuli. However, the plateau response during adaptation showed no significant changes. These findings indicate that the adaptation 
effects of V1 neurons may not be directly affected by surround suppression, but may be dynamically regulated by a negative 
feedback network and be finely adjusted by its initial spiking response to stimulus. This adaptive regulation is not only energy 
efficient for the central nervous system, but also beneficially acts to maintain the homeostasis of neuronal response to long-presenting 
visual signals. 


Keywords: Surround suppression; V1 neurons; Response adaptation; Cat 


Viewing a long-presenting visual stimulus with 
specific patterns (e.g., orientation, motion direction 
and spatial frequency) often inhibits or perturbs 
perception of a subsequent test stimulus with similar 
attributes. This 
adaptation, and has attracted considerable attention 
since the 1960s (Clifford et al, 2007; Dao et al, 2006; 
Greenlee & Heitger, 1988; Hua et al, 2009; Kohn, 
2007; Maffei et al, 1973; Marlin et al, 1988; Movshon 
& Lennie, 1979; Smith & Hammond, 1985). Since 
visual adaptation shows both evident interocular 


phenomenon is termed visual 


transfer and specificity to adapted stimulus attributes, 
it is generally regarded as a physiological process 
occurred in the cortical level, especially in the primary 
visual cortex (V1) (Duong & Freeman, 2007; Howarth 
et al, 2009), although subcortical neurons also exhibit 
a weak adaptation to visual stimulus (Brown & 
Masland, 2001; DeBruyn & Bonds, 1986; Smirnakis et 
al, 1997). 


Science Press 


The neuronal mechanisms of adaptation to visual 
stimuli are still in debate (Hua et al, 2009; Kohn, 2007; 
Liu et al, 2013). The contrast gain control mechanism, 
which suggests a somatic afterhyperpolarization due to 
an increasing potassium ion current triggered by sodium 
ion influx during prolonged stimulation (Carandini & 
Ferster, 1997; Sanchez-Vives et al, 2000a; Sanchez- 
Vives et al, 2000b), cannot interpret the specificity of 
adaptation to stimulus attributes. Synaptic mechanisms 
can fully account for stimulus-specificity of adaptation 
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but lack of consistent experimental evidences. Some 
studies highlight the roles of excitatory synaptic 
depression in the adaptation process (Chung et al, 2002; 
McLean and Palmer, 1996; Reig et al, 2006; Vidyasagar, 
1990). Some suggest that the changes of inhibitory 
synaptic activities may contribute to the adaptation 
effects (Hua et al, 2009; Yang et al, 2003). Others 
propose that adaptation may be caused by a network 
mechanism concerning a relative weight of recurrent 
excitation and inhibition in local neural circuitry (Teich 
& Qian, 2003). An important factor underlying these 
discrepancies is that previous studies fail to directly 
assesse the correlations of the changes of neuronal 
response adaptation and the changes of local excitation / 
inhibition. Studies on microiontophoresis found that 
administrations of glutamate, gamma-aminobutyric acid 
(GABA) and GABA receptor’s antagonists fail to change 
the adaptation strength of visual cortical neurons 
(DeBruyn & Bonds, 1986; Vidyasagar, 1990). However, 
it is premature to conclude that inhibition is not involved 
in the adaptation process because: 1) the actual effects of 
iontophoretic drug delivery may be challenged if drug 
diffusion time, diffusion range and synaptic spatial 
alignment were concerned; ii) regulatory mechanisms 
from inhibitory synapses other than GABAergic ones 
may exist in adaptation (Ego-Stengel et al, 2002; 
McLean & Palmer, 1996; Waterhouse et al, 1990). 
Moreover, we recently found that relative to young 
adults, the adaptation of V1 neurons in the aged brain 
with compromised intracortical inhibition is actually 
enhanced (Hua et al, 2009). 

The spiking response of a Vl neuron to a high- 
contrast stimulus placed within its classical receptive 
field (CRF) can be suppressed by a simultaneously 
presented stimulus within the surrounding nonclassical 
receptive field (nCRF), especially by the one with the 
similar orientation, motion direction and spatial frequ- 
ency (Cavanaugh et al, 2002b; Haider et al, 2010; Series 
et al, 2003; Webb et al, 2005). This phenomenon, termed 
surround suppression, is induced by the increased activ- 
ation of local inhibitory interneurons that are driven 
chiefly by the lateral horizontal connections and / or the 
feedback from higher visual cortical areas (Bair et al, 
2003; Durand et al, 2007; Haider et al, 2010; Li & 
Freeman, 2011; Series et al, 2003; Smith et al, 2006). 
Therefore, the local inhibition (Akasaki et al, 2002; Fu et 
al, 2010; Walker et al, 2000) on the surround-suppressed 
neurons can be regulated by the varying stimulus size 


Zoological Research 


outside the CRF. 

In this study, grating stimuli of different sizes were 
presented outside the CRF to evaluate the effects of 
intracortical inhibition on the response adaptation of V1 
neurons. 


MATERIALS AND METHORDS 


Animals 

Four healthy young adult cats (2—3 years old) were 
examined ophthalmoscopically prior to experimentation 
to confirm that no optical or retinal problems impaired 
their visual function. All experiment procedures were 
performed strictly in accordance with the guidelines 
published in the NIH Guide for the Care and Use of 
Laboratory Animals. 


Electrophysiological recording preparation 

All cats were prepared for acute in vivo single-unit 
recording using a previously described method (Hua et al, 
2010; Hua et al, 2009; Hua et al, 2006; Meng et al, 2013). 
Briefly, anesthesia was induced by injection of ketamine 
НСІ (40 mg/kg, im) and xylazine (2 mg/kg, im). After 
intubation of intravenous and tracheal cannulae, the cat 
was immobilized in a stereotaxic apparatus with ear, eye 
and bite bars. Glucose (5%)-saline (0.9%) solution cont- 
aining a mixture of urethane (20 mg/h/kg) and gallamine 
triethiodide (10 mg/h per kg of body weight) was 
infused intravenously by a syringe pump to keep the 
animal anesthetized and paralyzed. Pupils were 
maximally dilated with atropine (1%) eye drops, and 
contact lenses (zero power) were used to protect the 
corneas from dryness. Neosynephrine (5%) was applied 
to retract the nictitating membranes. Artificial respiration 
was performed, and expired pCO, was maintained at 
approximately 3.8%. Anesthesia level was closely 
evaluated during the experiment by continuously 
monitoring the animal’s heart rate (180—220 pulses/min) 
and electrocardiogram (ECG) throughout the experiment. 

V] was partly exposed (8 mm posterior to the 
earbar, 4 mm lateral to the midline) by removing the 
skull and dura over V1 (area 17) with the aid of a light 
microscope (77019, Reward, China). The small hole over 
VI was filled with 4% agar saline solution prior to 
electrophysiological recording. The optic discs of the two 
eyes were reflected onto a movable transparent tangent 
screen positioned 57 cm from the animal's eyes and 
overlapped with a CRT monitor (resolution 1024x768, 
refresh rate 85 Hz) for visual stimuli presentation. The 
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area centralis of each eye was precisely located 
according to the position of the optic discs reflected onto 
the tangent screen (Bishop et al, 1962). After all the 
preparations were completed, single-unit recordings were 
performed using a glass-coated tungsten microelectrode 
(with an impedance of 3—5 MQ) which was advanced by 
a hydraulic micromanipulator (Narishige, Japan). When 
the experiment was finished, the distance of each 
recorded cell’s receptive field from the retinal central 
area was measured and calculated as visual acuity 
(1°/cm). 


Visual stimuli and recording procedures 

Visual stimuli were drifting sinusoidal gratings, 
which were generated in MATLAB with the aid of 
extensions provided by the high-level Psychophysics 
Toolbox (Brainard, 1997) and low-level Video Toolbox 
(Pelli, 1997). Once a cell’s visually-evoked response was 
detected, the cell’s receptive field center was prelim- 
inarily determined using bars of light emitted from a 
hand pantoscope and then precisely located by consec- 
utively presenting a series of computer-generated flashi- 
ng bars of light on the CRT. The cell’s preferred stimulus 
attributes, including orientation, motion direction, spatial 
and temporal frequency were determined by comparing 
the cell’s response to a series of grating stimulus pack- 
ages. Then, the cell’s responses to grating stimuli with 
optimal attributes but different sizes were recorded to 
build the response-stimulus size tuning curve (Figure 
1A). We fitted the size tuning curve with a function 
described in previous papers (Cavanaugh et al, 2002a; 
Tailby et al, 2007): 


2 
хе kerf (x/w,) 


= 1 
1+k,erf (x/w, )? E 


Where, x is stimulus size, R(x) is the neuronal response 
to a stimulus with size х, kc and w, are the gain and 
spatial extent of the center mechanism, k, and w, are the 
gain and spatial extent of the surround mechanism, erf is 
the error function. 

From the fitting curve, we acquired three test 
stimulus sizes (a, b, c; Figure 1A), at which the cell’s 
response reached maximum, half of the maximum and 
minimum value on the right side of the fitting curve, 
respectively. Size a is the optimal stimulus size that only 
stimulates the cell’s CRF, but not induces surround 
suppression. Stimulus with size b and c can co-stimulate 
both CRF and nCRF, but may also evoke medium and 
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maximum surround suppression, respectively. 

The contrast of each stimulus was set at 10096. The 
mean luminance of the display was 19 cd/m?, and the 
environmental ambient luminance on the cornea was 0.1 lux. 


Data acquisition and analysis 

Action potentials of the recorded cells were 
amplified with a  microelectrode amplifier (Nihon 
Kohden, Japan) and a differential amplifier (Dagan 
2400A, USA), and then fed into a window discriminator 
with an audio monitor. The original voltage traces 
(Figure 1C, E, G) were digitized by an acquisition board 
(National Instruments, USA) controlled by IGOR 
software (WaveMetrics, USA), and saved for on- or off- 
line analysis. A cell's response to a grating stimulus was 
defined as the mean firing rate (spontaneous response 
subtracted) corresponding to the time of stimulus 
presentation, which was used to acquire the curves of 
tuning response to stimulus orientations, temporal and 
spatial frequencies. The optimal orientation of each cell 
was obtained as previously described. The optimal 
temporal and spatial frequency were determined 
respectively by comparing the cell's response to high 
contrast (10096) grating stimuli with different temporal 
and spatial frequencies, and selecting the temporal and 
spatial frequency with the maximum response. 

The adaptation index (AI) was defined as the ratio 
of the cell's mean response during plateau period of 
adaptation to visual stimulation, a period when the cell's 
response reached a stable minimum value, to the mean 
initial response of the cell (Figure 1D, F, H). The change 
of AI with different stimulus sizes was plotted for each 
studied cell (Figure 1B). The smaller the AI is, the 
stronger the adaptation of the cell becomes. In order to 
assess the impact of surround suppression on the 
response adaptation, several neurons that did not exhibit 
surround suppression to visual stimuli presented in its 
nCRF were excluded from our data analysis. All studied 
neurons had a receptive field within 8? from the central 
area of the dominant eye. 





All values were expressed as mean+S£. Variations 
between different stimulus sizes and subjects were 
assessed using analysis of variance (ANOVA) or t-test. 


RESULTS 


A total of 61 V1 cells from four young male adult 
cats were analyzed in this study (Table 1). АП cells 
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Figure 1 Response adaptation of a sample V1 cell to visual stimuli 

A: The response-size tuning curve fitted with equation 1. B: AI changes with stimulus sizes. The AI at stimulus sizes a, b and c for this 
sample cell was 0.44, 0.57 and 0.73 respectively. Subsequently, the cell's response to prolonged stimulation (90 stimulus cycles) were 
recorded with three stimulus sizes, respectively, which were used to assess the cell's response adaptation changes with the magnitude of 
surround suppression. Each stimulus was presented monocularly to the dominant eye and repeated 4—6 times with a 3-minute interval 
between adjacent trials for the cell's functional recovery. Before each stimulus was presented, spontaneous activity was acquired during a 10 
s period while a mean luminance was shown on the CRT. C, E and G: Voltage traces of the sample cell’s response to 90 cycles of preferred 
visual stimuli with size a, b and c respectively, which were employed to evaluate the cell's response adaptation strength. Spontaneous activity 
was obtained during the first 10 s period while mean luminance was shown on the screen. The dashed horizontal line in each voltage trace 
indicated the threshold for action potential counting. D, F and H: PSTHs show the cell's average response (counted across each 5 stimulus 
cycles, with spontaneous activity subtracted) changes as a function of time. Spikes in the first bar were defined as the average initial response, and 
the mean spikes from the 7th to the 17th bar as an average response, a period when the cell's response decreased to a stable minimum level. 


showed an evident adaptation to prolonged visual stimuli 
(90 stimulus cycles) as indicated by the AI value ranged 
from 0.104 to 0.760. 


Changes of neuronal response adaptation with the 
stimulus size outside the CRF 


The comparison of mean AI of studied neurons with 
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three stimulus sizes (a, b and c) showed that the surround 
suppression effects on the adaptation strength of 
neuronal response to visual stimuli from weak to strong 
were a, b and c, respectively. The ANOVA analysis 
showed significant differences in the averaged AI value 
of all the studied neurons with three different stimulus 
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Table 1 Mean adaptation index of V1 neurons at different 
stimulus sizes in each cat 








Subject Cell number (n) Adaptation index at different stimulus size 
a b с 

Catl 19 0.33+0.021 0.41+0.028 — 0.45+0.028 

Cat2 17 0.33+0.021 0.40+0.020 0.48+0.025 

Cat3 13 0.280.034 0.36+0.036 — 0.46::0.033 

Cat4 12 0.18+0.018 0.22+0.025 0.290.032 


а, b and с: the stimulus size at which the cell’s response reached maximum, 
half of the maximum and minimum value on the right side of the response- 


stimulus size fitting curve, respectively. 


sizes (Fo, 183)=25.7, P<0.0001). These differences were 
independent of subjects (Кв, 183)=0.38, Р>0.5), although 
the mean AI exhibited a significant variance from cat to 
cat (Fe, 183)=22.4, Р<0.0001) (Figure 2A). The mean AI 
of each individual cat was found significantly different 
with different stimulus sizes (сай: Fo, 5776.12, Р<0.01; 
cat2: Ео,51)-9.472, Р<0.001; cat3: Fe, 396.927, P<0.01; 
сай: Fo, 364.935, Р<0.05). The mean AI at stimulus 
size b was significantly less than that at stimulus size c 
(t-test, P«0.0001), whereas, was significantly larger than 
that at stimulus size a (t-test, P«0.0001) (Figure 2B, C), 
indicating that the neuronal response adaptation decr- 
eased with the increase of surround suppression. These 
results suggest that the response adaptation of V1 neur- 
ons to visual stimuli is negatively correlated with the 
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Figure 2 Average AI changes of studied neurons with different 
stimulus size 

A: mean AI at different stimulus sizes of individual cat; B and C: AI at 

stimulus size a v.s. size b and size c vs. size b of all the cells from all cats, 


respectively. 
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surround suppression that was modulated by the stimulus 
size. 


Changes of neuronal response with the stimulus size 
The increase or decrease of AI could result from a 
change of the initial response (IR) of neurons to visual 
stimuli, a change of the plateau response (PR), the 
response during the plateau period of adaptation, or any 
combination thereof. As such, we compared the IR (the 
mean response to the first five cycles of visual stimuli) 
and PR (the mean response of visual stimuli cycles from 
the 36" to the 85", which represents a minimal and stable 
response after adaptation) of V1 neurons to prolonged 
visual stimuli with different stimulus sizes, respectively. 
ANOVA analysis showed significant differences in 
the averaged IR of the studied neurons with different 
stimulus sizes (Fo, 183)=111.207, P«0.0001 =. These 
differences were independent of subjects (s 183= 0.536, 
P>0.5), although the mean IR varied significantly from 
cat to cat (Fo, 183)=14.633, Р<0.0001) (Figure ЗА). The 
mean IR of each individual cat also showed significant 
differences at different stimulus sizes (catl: Fo, 57)= 
21.810, Р<0.0001; cat2: Fo, 51)=75.240, Р<0.0001; cat3: 
Ёозоу= 23.203, P<0.0001; cat4: Fo жу-20.203, P<0.0001). 
The mean IR at stimulus size b was significantly 
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Figure 3 Average IR of neurons to prolonged visual stimuli 
with different sizes 
A: The mean IR at different stimulus sizes of individual cat; B and C: The 
IR at stimulus size a vs. size b and size c vs. size b of all the cells from all 


cats. 
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larger than that at stimulus size c (t-test, P<0.0001), 
whereas, was significantly less than that at stimulus size 
a (t-test, P«0.000001) (Figure 3B, C), indicating that the 
IR of neurons to prolonged visual stimuli reduced greatly 
with the increase of surround suppression. 

However, although the mean PR varied significantly 
from cat to cat (Fo, 1830.667, P«0.0001), no significant 
differences were found in either all the studied neurons 
from all the cats (Ро, 10.403, Р>0.1) or individual cat 
(сай: Fo, 57у=0.466, Р>0.5; саб: Fo. 5)=1.173, Р>0.1; 
Cat3: Fo, з9)=0.561, Р>0.5; cat4: Ер, 360.237, Р>0.5) with 
different stimulus sizes (Figure 4). These results indicate 
that the responses of neurons during the plateau period of 
adaptation to visual stimuli are stable and do not change 
significantly with the changes of surround suppression. 
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Figure 4 Average plateau response of neurons to prolonged 
visual stimuli with different sizes 
A: The mean plateau response at different stimulus sizes for individual cat; 
B and C: the plateau response at stimulus size a vs. size b and size c vs. size 


b for all the cells from all cats. 


Therefore, the response adaptation changes of the 
studied neurons with different surround suppression may 
attribute to the changes of their IR to prolonged visual 
stimuli, whereas, the PR maintains relatively stable. 


DISCUSSION 


Adaptation mechanisms 
Visual cortical neurons exhibit a reduction in firing 
rate to prolonged visual stimulation. However, the 
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underlying mechanisms remain in debate, and previous 
studies proposed several hypotheses. For example, the 
response adaptation is caused by the activity fatigue of 
the neuron because the prolonged stimulation may evoke 
sustained firing and the fatigued neurons respond less 
than they normally do (Carandini, 2000; Sekuler & 
Pantle, 1967). The contrast gain control mechanism 
suggests that adaptation leads to a strong somatic 
afterhyperpolarization due primarily to the activation of 
voltage-gated potassium channels, triggered by the 
sodium influx during generation of action potentials 
(Carandini & Ferster, 1997; Sanchez-Vives et al, 2000a; 
Sanchez-Vives et al, 2000b). Although the above two 
mechanisms can interpret the neuronal response 
reduction during visual adaptation, they are unfortunately 
unable to account for the specificity of adaptation to the 
adapted stimulus attributes, such as orientation. 

To date, more studies support the synaptic mecha- 
nism due to its advantage in explaining the stimulus- 
specificity of visual adaptation. However, debates 
concerning the contributions of excitation and inhibition 
in local circuitry to the adaptation still exist. Some 
studies emphasize the importance of excitatory synaptic 
depression in mediating the adaptation process (Chung et 
al, 2002; Nowak et al, 2005; Reig et al, 2006). Some 
suggest an involvement of inhibitory synaptic activation 
in the adaptation effect, and others propose a network 
mechanism based on recurrent excitation and inhibition 
models (Kohn, 2007). Current evidences on the role of 
local inhibition in the adaptation process are mutually 
inconsistent. An in vivo study reported that the iontop- 
horetic delivery of GABA, receptor antagonists could 
significantly improve the adaptation strength of relay 
cells in the dorsal geniculate nucleus (LGNd) and the 
administration of baclofen, a GABAg receptor agonist, 
could decrease the adaptation strength (Yang et al, 2003). 
However, similar manipulation of GABA inhibition 
failed to alter the amplitude of visual cortical cells to 
visual stimuli (DeBruyn & Bonds, 1986; Vidyasagar, 
1990). Interestingly, V1 neurons of aged cats showed 
stronger adaptation to visual stimuli than that of young 
adult cats (Hua et al, 2009). This enhanced adaptation of 
V] neurons during aging may indirectly suggest that 
neuronal response adaptation is correlated with the 
changes of intracortical inhibition (Hua et al, 2008; Hua 
et al, 2006; Leventhal et al, 2003). 

In the present study, we evaluated the effects of 
surround suppression on the adaptation strength of V1 
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neurons. By manipulating the levels of surround 
suppression using different stimulus sizes outside the 
CRF, we found that the amplitude of response adaptation 
of V1 neuron decreased significantly with the increase of 
surround suppression due to the decrease of the neuron's 
IR to the adapted stimulus. These results indicate that the 
surround suppression might only modify a neuron's IR 
but not the adaptation process and the adaptation strength 
depends closely on the neuron’s IR to the adapted 
stimulus. IR decreases with the increase of surround 
suppression and the decrease of response adaptation, vice 
versa. Therefore, the response adaptation of V1 neurons 
may under the dynamic regulation of a negative feedback 
mechanism. Our results, together with several recent 
findings (Benucci et al, 2013; Compte & Wang, 2006; 
Levy et al, 2013; Liu et al, 2013) suggest that visual 
adaptation may depend on a network mechanism that 
involves an interplay between inhibitory and excitatory 
neurons in the local neural circuitry. 


Benefits of visual adaptation 

The functional benefits of adaptation remain unclear 
due to the inconsistent evidences suggest that adaptation 
sometimes decreases sensitivity for the adapting stimuli, 
and sometimes it changes sensitivity for stimuli very 
different from the adapting ones (Gepshtein et al, 2013). 
Some studies claimed that adaptation could improve the 
detectability of the adapting stimuli (Abbonizio et al, 
2002; Greenlee & Heitger, 1988; Määttänen 4 Koend- 
erink, 1991). Others reported that adaptation increased 


References 


Abbonizio G, Langley K, Clifford CW. 2002. Contrast adaptation may 
enhance contrast discrimination. Spatial Vision, 16(1): 45-58. 


Akasaki T, Sato H, Yoshimura Y, Ozeki H, Shimegi S. 2002. 
Suppressive effects of receptive field surround on neuronal activity in 
the cat primary visual cortex. Neuroscience Research, 43(3): 207-220. 


Bair W, Cavanaugh JR, Movshon JA. 2003. Time course and time- 
distance relationships for surround suppression in macaque V1 neurons. 
Journal of Neuroscience, 23(20): 7690-7701. 


Benucci A, Saleem AB, Carandini M. 2013. Adaptation maintains 
population homeostasis in primary visual cortex. Nature Neuroscience, 
16(6): 724-729. 


Bishop PO, Kozak W, Vakkur GJ. 1962. Some quantitative aspects of 
the cat's eye: axis and plane of reference, visual field co-ordinates and 
optics. Journal of Physiology, 163(3): 466-502. 


Brainard DH. 1997. The psychophysics toolbox. Spatial Vision, 10(4): 
433-436. 


Kunming Institute of Zoology (CAS), China Zoological Society 


perception of novel stimuli in the environment while 
suppressing the perception of adapted stimuli (Dragoi et 
al, 2002; Hosoya et al, 2005; Sharpee et al, 2006). 
Benucci et al (2013) measured adaptation in the response 
of populations of V1 neurons to stimulus ensembles with 
markedly different statistics of stimulus orientation, and 
found that adaptation might act as a mechanism of 
homeostasis by maintaining time-averaged response 
quality and orientation selectivity independence across 
the population of neurons. 

In this study, we determined the response adaptation 
changes of V1 neurons with different degree of surround 
suppression. We found that the response adaptation of V1 
neurons decreased significantly with the increase of 
surround suppression due to the reduction of its IR to the 
adapted stimulus, whereas, the response of neurons 
during the plateau period of adaptation remained stable. 
These results are consistent with previous studies 
(Cavanaugh et al, 2002a) and suggest that V1 neurons 
may dynamically adjust its adaptation strength 
according to its initial spiking activities evoked by the 
adapted stimulus: adaptation enhances if initial 
activities are high or otherwise weakens if initial 
activities are low. The response of the neuron can 
eventually be reduced to the similar level, which is 
independent of the amplitude of initial response. This 
adaptation strategy may be critical in maintaining the 
homeostasis of neuronal response to long-lasting visual 
signals and aiding the energy efficiency/frugality of 


brain activities (Hua et al, 2009). 


Brown SP, Masland RH. 2001. Spatial scale and cellular substrate of con- 
trast adaptation by retinal ganglion cells. Nature Neuroscience, 4(1): 44-51. 


Carandini M. 2000. Visual cortex: Fatigue and adaptation. Current 
Biology, 10(16): R605-607. 


Carandini M, Ferster D. 1997. A tonic hyperpolarization underlying 
contrast adaptation in cat visual cortex. Science, 276(5314): 949-952. 


Cavanaugh JR, Bair W, Movshon JA. 2002a. Nature and interaction of 
signals from the receptive field center and surround in macaque V1 
neurons. Journal of Neurophysiology, 88(5): 2530-2546. 


Cavanaugh JR, Bair W, Movshon JA. 2002b. Selectivity and spatial 
distribution of signals from the receptive field surround in macaque V1 
neurons. Journal of Neurophysiology, 88(5): 2547-2556. 


Chung S, Li X, Nelson SB. 2002. Short-Term depression at 
thalamocortical synapses contributes to rapid adaptation of cortical 


sensory responses in Vivo. Neuron, 34(3): 437-446. 


Volume 35 Issue 5 


418 LI, et al. 


Clifford CW, Webster MA, Stanley GB, Stocker AA, Kohn A, Sharpee 
TO, Schwartz O. 2007. Visual adaptation: neural, psychological and 
computational aspects. Vision Research, 47(25): 3125-3131. 


Compte A, Wang XJ. 2006. Tuning curve shift by attention modulation 
in cortical neurons: a computational study of its mechanisms. Cerebral 
Cortex, 16(6): 761-778. 


Dao DY, Lu ZL, Dosher BA. 2006. Adaptation to sine-wave gratings 
selectively reduces the contrast gain of the adapted stimuli. Journal of 
Vision, 6(7): 739-759. 


DeBruyn EJ, Bonds AB. 1986. Contrast adaptation in cat visual cortex 
is not mediated by GABA. Brain Research, 383(1-2): 339-342. 


Dragoi V, Sharma J, Miller EK, Sur M. 2002. Dynamics of neuronal 
sensitivity in visual cortex and local feature discrimination. Nature 
Neuroscience, 5(9): 883-891. 


Duong T, Freeman RD. 2007. Spatial frequency-specific contrast 
adaptation originates in the primary visual cortex. Journal of 
Neurophysiology, 98(1): 187-195. 


Durand S, Freeman TC, Carandini M. 2007. Temporal properties of 
surround suppression in cat primary visual cortex. Vision Neuroscience, 
24(5): 679-690. 


Ego-Stengel V, Bringuier V, Shulz DE. 2002. Noradrenergic 
modulation of functional selectivity in the cat visual cortex: an in vivo 
extracellular and intracellular study. Neuroscience, 111(2): 275-289. 


Fu Y, Wang XS, Wang YC, Zhang J, Liang Z, Zhou YF, Ma YY. 2010. 
The effects of aging on the strength of surround suppression of 
receptive field of V1 cells in monkeys. Neuroscience, 169(2): 874-881. 


Gepshtein $, Lesmes LA, Albright TD. 2013. Sensory adaptation as 
optimal resource allocation. Proceedings of the National Academy of 
Sciences of the United States of America, 110(11): 4368-4373. 


Greenlee MW, Heitger F. 1988. The functional role of contrast 
adaptation. Vision Research, 28(7): 791-797. 


Haider B, Krause MR, Duque A, Yu Y, Touryan J, Mazer JA, 
McCormick DA. 2010. Synaptic and network mechanisms of sparse 
and reliable visual cortical activity during nonclassical receptive field 
stimulation. Neuron, 65(1): 107-121. 


Hosoya T, Baccus SA, Meister M. 2005. Dynamic predictive coding by 
the retina. Nature, 436(7047): 71-77. 


Howarth CM, Vorobyov V, Sengpiel F. 2009. Interocular transfer of 
adaptation in the primary visual cortex. Cerebral Cortex, 19(8): 1835- 
1843. 


Hua Т, Bao Р, Huang CB, Wang 7, Xu J, Zhou Y, Lu ZL. 2010. 
Perceptual learning improves contrast sensitivity of V1 neurons in cats. 
Current Biology, 20(10): 887-894. 


Hua TM, Kao CC, Sun QY, Li XR, Zhou YF. 2008. Decreased 
proportion of GABA neurons accompanies age-related degradation of 
neuronal function in cat striate cortex. Brain Research Bulletin, 75(1): 
119-125. 


Hua TM, Li GZ, Tang CH, Wang ZH, Chang S. 2009. Enhanced 
adaptation of visual cortical cells to visual stimulation in aged cats. 
Neuroscience Letters, 451(1): 25-28. 


Zoological Research 


Hua TM, Li XR, He LH, Zhou YF, Wang YC, Leventhal AG. 2006. 
Functional degradation of visual cortical cells in old cats. Neurobiology 
of Aging, 27(1): 155-162. 


Kohn A. 2007. Visual adaptation: physiology, mechanisms, and 
functional benefits. Journal of Neurophysiology, 97(5): 3155-3164. 


Leventhal AG, Wang Y, Pu M, Zhou Y, Ma Y. 2003. GABA and its 
agonists improved visual cortical function in senescent monkeys. 
Science, 300(5620): 812-815. 


Levy M, Fournier J, Fregnac Y. 2013. The role of delayed suppression 
in slow and fast contrast adaptation in V1 simple cells. Journal of 
Neuroscience, 33(15): 6388-6400. 


Li B, Freeman RD. 2011. Neurometabolic coupling differs for 
suppression within and beyond the classical receptive field in visual 
cortex. Journal of Physiology, 589(Pt 13): 3175-3190. 


Liu RL, Wang K, Meng JJ, Hua TM, Liang Z, Xi MM. 2013. 
Adaptation to visual stimulation modifies the burst firing property of 
V1 neurons. Zoological Research, 34(3): E101-E108. 


Маанапеп LM, Koenderink JJ. 1991. Contrast adaptation and contrast 
gain control. Experimental Brain Research, 87(1): 205-212. 


Maffei L, Fiorentini A, Bisti S. 1973. Neural correlate of perceptual 
adaptation to gratings. Science, 182(4116): 1036-1038. 


Marlin SG, Hasan SJ, Cynader MS. 1988. Direction-selective 
adaptation in simple and complex cells in cat striate cortex. Journal of 
Neurophysiology, 59(4): 1314-1330. 


McLean J, Palmer LA. 1996. Contrast adaptation and excitatory amino 
acid receptors in cat striate cortex. Visual Neuroscience, 13(6): 1069- 
1087. 


Meng JJ, Liu RL, Wang K, Hua TM, Lu ZL, Xi MM. 2013. Neural 
correlates of stimulus spatial frequency-dependent contrast detection. 
Experimental Brain Research, 225(3): 377-385. 


Movshon JA, Lennie P. 1979. Pattern-selective adaptation in visual 
cortical neurones. Nature, 278(5707): 850-852. 


Nowak LG, Sanchez-Vives MV, McCormick DA. 2005. Role of 
synaptic and intrinsic membrane properties in short-term receptive 
field dynamics in cat area 17. Journal of Neuroscience, 25(7): 
1866-1880. 


Pelli DG. 1997. The VideoToolbox software for visual psychophysics: 
transforming numbers into movies. Spatial Vision, 10(4): 437-442. 


Reig R, Gallego R, Nowak LG, Sanchez-Vives MV. 2006. Impact of 
cortical network activity on short-term synaptic depression. Cerebral 
Cortex, 16(5): 688-695. 


Sanchez-Vives MV, Nowak LG, McCormick DA. 2000a. Cellular 
mechanisms of long-lasting adaptation in visual cortical neurons in 
vitro. Journal of Neuroscience, 20(11): 4286-4299. 


Sanchez-Vives MV, Nowak LG, McCormick DA. 2000b. Membrane 
mechanisms underlying contrast adaptation in cat area 17 in vivo. 
Journal of Neuroscience, 20(11): 4267-4285. 


Sekuler R, Pantle A. 1967. A model for after-effects of seen movement. 
Vision Research, (5): 427-439. 


WWNW.ZOOres.ac.cn 


Effects of surround suppression on response adaptation of V1 neurons to visual stimuli 419 


Series P, Lorenceau J, Fregnac Y. 2003. The "silent" surround of V1 
receptive fields: theory and experiments. Journal of Physiology, 97(4- 
6): 453-474. 


Sharpee TO, Sugihara H, Kurgansky AV, Rebrik SP, Stryker MP, Miller 
KD. 2006. Adaptive filtering enhances information transmission in 
visual cortex. Nature, 439(7079): 936-942. 


Smirnakis SM, Berry MJ, Warland DK, Bialek W, Meister M. 1997. 
Adaptation of retinal processing to image contrast and spatial scale. 
Nature, 386(6620): 69-73. 


Smith AT, Hammond P. 1985. The pattern specificity of velocity 
aftereffects. Experimental Brain Research, 60(1): 71-78. 


Smith MA, Bair W, Movshon JA. 2006. Dynamics of suppression in 
macaque primary visual cortex. Journal of Neuroscience, 26(18): 4826- 
4834. 


Tailby C, Solomon SG, Peirce JW, Metha AB. 2007. Two expressions 
of "surround suppression" in V1 that arise independent of cortical 
mechanisms of suppression. Visual Neuroscience, 24(1): 99-109. 


Kunming Institute of Zoology (CAS), China Zoological Society 


Teich AF, Qian N. 2003. Learning and adaptation in a recurrent model of V1 
orientation selectivity. Journal of Neurophysiology, 89(4): 2086-2100. 


Vidyasagar TR. 1990. Pattern adaptation in cat visual cortex is a co- 
operative phenomenon. Neuroscience, 36(1): 175-179. 


Walker GA, Ohzawa I, Freeman RD. 2000. Suppression outside the 
classical cortical receptive field. Visual Neuroscience, 17(3): 369-379. 


Waterhouse BD, Azizi SA, Burne RA, Woodward DJ. 1990. 
Modulation of rat cortical area 17 neuronal responses to moving visual 
stimuli during norepinephrine and serotonin microiontophoresis. Brain 
Research, 514(2): 276-292. 


Webb BS, Dhruv NT, Solomon SG, Tailby C, Lennie P. 2005. Early and 
late mechanisms of surround suppression in striate cortex of macaque. 
Journal of Neuroscience, 25(50): 11666-11675. 


Yang Y, Jin J, Zhou Y, Shou T. 2003. GABA(A) and GABA(B) 
receptors mediated inhibition affect the pattern adaptation of relay cells 
in the dorsal lateral geniculate nucleus (LGNd) of cats. Brain Research, 
959(2): 295-303. 


Volume 35 Issue 5 


Zoological Research 35 (5): 420—425 рог10.13918/ 155п.2095-8137.2014.5.420 





Protease-activated receptor (PAR)1, PAR2 and PARA expressions in 
esophageal squamous cell carcinoma 


Si-Man LI', Ping JIANG’, Yang XIANG’, Wei-Wei WANG”, Yue-Chun ZHU’, Wei-Yang FENG’, 
Shu-De LI', Guo-Yu YU" 


1. Kunming Medical University, Kunming 650500, China 

2. Key Laboratory of Animal Models and Human Disease Mechanisms of the Chinese Academy of Sciences & Yunnan Province, Kunming Institute of Zoology, 
Kunming 650223, China 

3. Third Affiliated Hospital of Kunming Medical University, Department of Thoracic Surgery, Kunming 650000, China 


Abstract: Here, we used reverse transcription-PCR (RT-PCR) and western blot to detect protease-activated receptor (PAR) 1, PAR 2 
and PAR 4 expression in cancer tissues and cell lines of esophageal squamous cell carcinoma, and investigated the co-relationship 
between PAR expression and clinic-pathological data for esophageal cancer. The methylation of PAR4 gene promoter involved in 
esophageal carcinoma was also analyzed. By comparing the mRNA expressions of normal esophageal tissue and human esophageal 
epithelial cells (HEEpiC), we found that among the 28 cases of esophageal squamous cell carcinoma, PARI (60%) and PAR2 (71%) 
were elevated in 17 and 20 cases, respectively, and PAR4 (68%) expression was lowered in 19 cases. Whereas, in human esophageal 
squamous cells (TE-1 and TE-10), PARI and PAR2 expression was increased but PAR4 was decreased. Combined with clinical data, 
the expression of PARI in poorly differentiated (P=0.016) and middle and lower parts of the esophagus (P=0.016) was higher; 
expression of PAR4 in poorly differentiated carcinoma was lower (P=0.049). Regarding TE-1 and TE-10 protein expression, we 
found that in randomized esophageal carcinoma, PARI (Р-0.027) and PAR2 (P=0.039) expressions were increased, but lowered for 
PARA (P=0.0001). In HEEpiC, TE-1, TE-10, esophageal and normal esophagus tissue samples (case No. 7), the frequency of 
methylation at the 19 CpG loci of PAR4 was 35.4%, 95.2%, 83.8%, 62.6% and 48.2%, respectively. Our results indicate that the 
expression of PARI and PAR2 in esophageal squamous cell carcinoma is increased but РАВА is decreased. Hypermethylation of the 
promoter of the PAR4 gene may contribute to reduced expression of PAR4 in esophageal squamous cell carcinoma. 


Keywords: PARI; PAR2; PARA; Esophageal squamous cell carcinoma; PCR; Western Blot; Methylation 


The development of esophageal squamous cell carc- 
inoma (ESCC) is a complicated process with multiple 
pathological stages. Among various regulatory factors, 
proteases play critical roles in activating signal trans- 
duction pathways and regulating gene expression (Ikeda 
et al, 1999). Protease-activated receptors (PARs) are a 
subfamily of the single seven-transmembrane G-protein- 
coupled receptors and include PARI, PAR2, PAR3 and 
РАКА (Macfarlane et al, 2001). PARI, PAR3 and РАКА 
are thrombin receptors, and PAR2 15 trypsinase/tryptase 
receptor (Xu et al, 1998). Ribeiro et al (2009) found high 
expression of PARI but low expression of PAR2 in the 
tissues of ESCC; however, Wang et al (2010) found high 
expression of PAR2 in ESCC. PARA is a recently 
discovered novel subtype expressed in high amounts in 
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colon cancer and hepatocarcinoma and induces the 
proliferation and migration of cancer cells (Gratio et al, 
2009; Kaufmann et al, 2007). РАКА expression in stom- 
ach cancer tissues (Zhang et al, 2011) and adenocarc- 
inoma of the lung are low (Jiang et al, 2013), but its 
expression in ESCC remains unclear. 

Here, we determined the expression of PARI, PAR2 
and PARA in tissues of ESCC and analyzed correlations 
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between PAR4 promoter hypermethylation and the 
development of ESCC. Our aim was to provide 
theoretical evidence for clinical diagnosis, treatment and 
prognosis of ESCC and provide a basis for drug research. 


MATERIALS AND METHODS 


Experimental materials 

Human esophageal epithelial cells (HEEpiC) and 
human esophageal squamous cells (TE-1, TE-10) were 
obtained from the Cell Bank of Kunming Institute of 
Zoology, Chinese Academy of Sciences. Tissues of 28 cases 
(male=21, female=5; 51-81 years old) of diagnosed ESCC 
were from hospitals affiliated with the Kunming Medical 
University, Yunnan, China. All patients were clean of any 
chemotherapy or radiation treatment prior to surgery. The 
carcinoma tissues and corresponding normal control tissues 
(at least 5 cm away from carcinoma tissue) obtained during 
surgery were quick frozen in liquid nitrogen and stored at — 
80 ?C. AII the experimental protocols were approved by the 
Ethics Committee of Kunming Medical University. 


Experimental procedures Cell culture 

ТЕ-1 and TE-10 cells were cultured with RPM- 
11640 (Takara, Beijing, China) completed culture 
medium containing 10% FBS and incubated at 37 °С 
with 5% СО; (relative humidity=95%). Culture medium 
was changed 1—2 days later. Cells were subcultured 3—5 
days later and digested with 0.25% trypsinase: 0.03% 
EDTA (1:1) HEEpiC were cultured with HEEpiC- 
specific (Takara, Beijing, China) culture medium and the 
culture procedure was the same as that for TE-1 cells. 


Reverse transcription-PCR (RT-PCR) 

Total RNA of ESCC tissues, TE-1, TE-10 and HEE- 
piC were extracted using a RNA extraction kit (Tiangen 
Biotech, Beijing, China). The purity and integrity of total 
RNA were tested. cDNAs were reverse transcribed from 


total RNAs of tissues and cells (2 ng-2 ир) using a 
reverse transcription kit (Takara, Beijing, China) and 
then stored at 20 °C. Gene amplifications were perfor- 
med by taking cDNAs as templates and GAPDH as 
internal reference. Primers for GAPDH, PARI, PAR2 and 
PAR4 are shown in Table 1. PCR products were run by 
2% agarose gel electrophoresis and observed under 
ultraviolet light and photographed. 


Western blot 

Carcinoma tissues and corresponding normal control 
tissues were randomly selected from three cases. 
Cultured cells in logarithmic growth were rinsed twice 
with pre-chilled PBS (4 °C) and then lysed with 6х SDS 
(0.6 mL) and water-bathed (95 °C) for 15 min. After the 
polymerization of resolving gel (12%) and staking gel 
(5%), samples were loaded. Tris-glycine electrophoresis 
buffer (x1) was poured into the running chamber and 
loaded samples were run for 1 h at a constant current of 
160 V. Then the gel was transferred onto the PVDF 
membrane for 2 h at a constant current of 250 mA and 
4 °С. The blots were removed from the transfer unit and 
blocked by placing in 3% BSA-TBST blotting solution 
for 2 h with shaking at room temperature. After washing, 
the blot was incubated with the primary antibodies for 
PARI, PAR2, РАКА and f-actin over night at 4 °C. After 
washing, the blot was incubated with secondary 
antibodies. After another three washes, the blot was 
revealed via HRP-ECL chemiluminescence detection and 
scanned using the SynGene scanning system. The 
illumination densities of the detected protein bands were 
transferred into optical density values (OD values). The 
developing strength of each target protein band and its 
internal reference (f-actin) were determined by OD 
values. The ratio of the OD values of each target protein 
band to its internal reference (f.-actin) was taken as the 
result of RT-PCR. 


Table 1 Primers and reaction conditions used in RT-PCR 





Gene Primer sequences NCBI Accession Number Annealing temperature (?C) Product (bp) 

Е: 5-АТССССААССТСААССТСС-3' 
GAPDH NM 001101.3 60 308 

R: 5-GGGGTCATTGATGGCAACAATA-3' 
F: 5-GCCGCCTGCTTCAGTCTGTGC-3' 

PARI NM_001992.3 67 648 
R: 5-GGCCAGACAAGTGAAGGAAGC-3' 
Е: 5-CCATCCAAGGAACCAATAGATC-3' 

PAR2 NM_005242.3 60 643 
R: S-ATGTCTCCCACCAAGAGCTGCTCA-3' 
F: 5-GGCAACCTCTATGGTGCCTA-3' 

PAR4 NM_003950.2 58 244 


К: 5'-TTCGACCCAGTACAGCCTTC-3’ 
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Bisulfite genomic sequencing PCR 

Genomic DNAs of TE-1, TE-10, HEEpiC and the 
No. 7 tissue sample were extracted using a genome DNA 
extraction kit (Takara, Beijing, China) and their concent- 
rations were determined. The CpG sites in PAR4 were 
amplified with the MethylCode Bisulfite Conversion Kit 
(Invitrogen, US). PCR products were eluted and purified 
with the AxyPrep DNA elution kit (Axygen). Primer 
sequences were: F-5'-TTTAAGGGTGATTTTAGGAAA 
GGTTTAGAG-3'and R-5'-ACTATAACCTCAAACTTC 
CTACCTC-3’. The products of ligation were transformed 
with DH5a component cells. Transformation products 
were spread on LB plates (with Ampicillin) and grown 
overnight. Clones were selected and sent for sequencing. 


Statistical analyses 

Data were analyzed using SPSS 17.0 (SPSS Inc., 
Chicago, USA). Correlations between PAR expression 
and clinical pathologies were tested using Fisher's exact 
tests. The illumination strength of protein expression is 





expressed as теап+50 and comparisons between 
expressions were conducted using t-tests. Statistical 
significance was set at P«0.05. 


RESULTS 


mRNA expression of PARI, PAR2 and PAR4 and 
correlation with clinical manifestation 

When the expression of GAPDH of each sample 
was at a similar level, the brighter the mRNA band was, 
the higher the expression level of PAR was, and vice 
versa. Among the 28 cases of ESCC, 17 (60%) and 20 
(7196) cases were found with increased expression of 
PARI and PAR2, respectively. Totally, 19 (68%) cases 
were found with decreased PAR4. In TE-1 and TE-10, the 
expression of PAR/and PAR2 was high. The expression 
of PAR4 mRNA in TE-1 and TE-10 was lower than in 
HEEpiC. The expression level of PARI, PAR2 and PAR4 
mRNA in the tissues and cells from cases 7, 8 and 9 are 
shown in Figure 1. Combined with clinical pathological 
data, the expression of PAR/ in central and low located 
ESCC was significantly higher than in upper ESCC 
(Р=0.007 and P=0.008, respectively). The expression of 
PAR2 in phase III+IV ESCC was significantly higher 
than in phase I-II (P=0.004). Significantly lower 
expression of PAR4 was found in lower ESCC, compared 
with central and upper ESCC (P=0.036). No significant 
correlation for the expression of PAR/, PAR2 and PAR4 
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and patients’ gender, sex and distant or lymph node 
metastases were found (Table 2). 


Patient NO. _7 8 
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о У 
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PAR4 





PAR2 


PARI 


GAPDH 


Figure 1 PARI, PAR2 and PAR4 mRNA expression in 
esophageal squamous cell carcinoma and normal control tissues 


Protein expression of PARI, PAR2 and PAR4 in 
tissues and cells of ESCC 

Protein expression was determined by taking B-actin 
as the internal reference (Figure 2). The expression of 
PARI and PAR2 proteins were increased (Р<0.05), 
whereas, those of РАКА were decreased (P<0.05). The 
light degrees of protein expression of PARI, PAR2 and 
PARA in sample tissues are shown in Table 3. The 
expression of PARI and PAR2 were high and that of 
PAR4 was low in TE-1 and TE-10 cells. 


Patient NO. _7_ 8 9. 
PAR4 
PAR2 


PARI 


Actin 





Figure 2 Protein expression of PARI, PAR2 and РАКА in 
tissues and cells of esophageal squamous cell carcinoma 


PAR4 promoter hypermethylation in tissues and cells 
of ESCC 

The BSP results of the PCR products of CpG sites 
in PAR4 are shown in Table 4. In HEEpiC, TE-1 and TE- 
10, esophageal and normal esophagus tissue from case 
No. 7, the methylation frequency of the 19 CpG sites in 
PAR4 was 35.4%, 95.2%, 83.8%, 62.6% and 48.2%, 
respectively. 


DISCUSSION 


PAR2 over-expression is commonly found in 
malignant tumors. In this study, among 28 cases of 
ESCC, 20 (71%) cases were found with high PAR2 
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Table 2 Correlation between PARI, PAR2 and PAR4 mRNA levels and clinical manifestation of 
esophageal squamous cell carcinoma 











PARI PAR2 PAR4 
Clinical data Case numbers (л) 
Decreased Increased P-value Decreased Increased P-value Decreased Increased P-value 
Age (year) 
<65 13 5 8 3 10 8 5 
1.000 0.686 0.689 
>65 15 6 9 5 10 11 4 
Gender 
Male 21 9 12 4 17 15 6 
0.668 0.142 0.646 
Female 7 2 5 4 3 4 3 
Clinical stage 
I-II 16 8 8 7 9 11 5 
0.253 0.088 1.000 
I+IV 12 3 9 1 11 8 4 
Location of the tumor 
Upper 4 4 0 У 0 4 2 2 
0.016 0.549 0.574 
Central and lower 24 7 17 5 16 17 7 
Differentiation 
Well and moderated 17 10 7 7 10 9 8 қ 
0.016 0.099 0.049 
Роог 11 1 10 1 10 10 1 
Distant metastasis 
Positive 1 0 1 0 1 1 0 
1.000 1.000 1.000 
Negative 27 11 16 8 19 18 9 
Lymph node metastasis 
Positive 11 4 7 2 9 8 3 
1.000 0.419 0.704 
Negative 17 4 10 6 11 11 6 





^: P<0.05. 


Table 3 Gray values of PAR expression in sample tissues 





Gene Cancer tissue Normal tissue t P 
PAR4 0.18+0.05 0.42+0.12 9.77 0.0001“ 
PAR2 0.84+0.28 0.69+0.25 7.97 0.039" 
PARI 0.74+0.31 0.58+0.21 2.26 0.027" 


*: P<0.05; 7: P«0.01. 


mRNA expression. The expression of PAR2 mRNA was 
also increased in TE-1 and TE-10 cells. The gray value 
indicates that PAR2 protein expression is increased in 
ESCC (P=0.039). PAR2 can be activated by trypsinase, 
human mast cell tryptase (MCT), coagulation factor УПа 
and tissue factor complex (Caruso et al, 2006; Uusitalo- 
Jarvinen et al, 2007). The in vitro synthesized SLIGKV 
by the degradation of PAR2 can also activate PAR2 
(Déry et al, 1998). Activated PAR2 may release vascular 
endothelial growth factor (VEGF), Interleukin-6 (IL-6) 
and IL-8 of tumor cells and thereafter promote the 
generation and invasion of novel vessels by malignant 
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tumors (Knecht et al, 2007; Matej et al, 2007). Activated 
PAR2 may also increase the reverse activity of epidermal 
growth factor (EGF) and the release of transformation 
growth factor a (TGFa), and thus promote the prolife- 
ration of gastrointestinal cancer cells in stomach cancer, 
colon cancer, pancreatic cancer and ESCC (Darmoul et al, 
2004; Fujimoto et al, 2006; Yada et al, 2005). 

PARI and PARA are both thrombin receptors. 
Among the 28 cases of ESCC examined in the present 
study, 17 (60%) cases were found with high expression 
of PARI mRNA; however, 19 (68%) cases were found 
with low expression of PAR4 mRNA. The expression of 
PARI mRNA was increased in TE-1 and TE-10 cells, 
whereas, that of PAR4 was decreased. The gray value 
indicates that РАКІ protein expression was increased in 
ESCC (P=0.027) and that of PAR4 was significantly 
decreased (P=0.0001). High PARI expression was found 
in ESCC, particularly in central and lower (P=0.016) and 
poorly differentiated ESCC (P=0.049), whereas, the low 
expression of РАКА was found in ESCC, particularly in 
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poorly differentiated ESCC. 


Table 4 Methylation frequency of the 19 CpG sites in PAR4 of cells and tissue of esophageal squamous cell carcinoma 





CpG position Gene 3 12 41 68 95 124 195 215 227 259 275 292 327 330 332 341 345 350 377 Total 
HEEPIC 18 73 0 64 27 55 73 9 46 0 46 55 27 0 0 0 73 55 55 354 

TE-1 100 100 90.9 100 100 81.8 100 30 100 100 82 91 100 100 91 91 100 100 100 952 

TET TE-10 82 100 100 100 100 100 100 91 100 100 46 73 55 55 100 100 100 55 30 38338 
Cancer tissue 64 46 30 46 64 9 64 100 91 100 100 100 91 91 9 4 46 30 64 626 

Normal tissue 91 64 9 55 64 30 55 9 55 100 100 9 0 55 55 30 64 30 46 482 


То investigate the underlying mechanisms of PAR4 
and ESCC, methylation of the PAR4 promoter in tissues 
and cells of ESCC was evaluated. The results show that 
the methylation frequencies of PAR4 promoters in 
sample tissues of ESCC (62.6%) and TE-1 and TE-10 
(83.8% and 95.2%, respectively) were both high. 
Together with the fact that the expression of PAR4 mRNA 
and protein is lower in tissues and cells of ESCC, these 
findings indicate that the methylation frequency of PAR4 
promoters may play a role in its expression in ESCC. 

Kawabata et al (1999) found the relaxation induced 
by РАКА aggravates duodenal-gastric-esophageal reflux, 
and trypsinase within the reflux liquid may activate 
PARA, thus promoting the incidence of ESCC. Han et al 
(2011) found that the promoting effects of PARA on anti- 
angiogenesis factors, including endostatin, thrombo- 
stondin 1, 02 macroglobulin, plasminogen activator, ang- 
lostatin and enzyme inhibitors may remarkably inhibit 
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Differences in cocaine-induced place preference persistence, loco- 
motion and social behaviors between C57BL/6J and BALB/cJ mice 
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Abstract: C57BL/6J and BALB/cJ mice display significant differences in sociability and response to drugs, but the phenotypic variability of 
their susceptibility to cocaine is still not well known. In this study, the differences between these two mice strains in the persistence of 
cocaine-induced conditioned place preference (CPP), as well as the locomotion and social behaviors after the 24-hour withdrawal from a 
four-day cocaine (20 mg/kg/day) administration were investigated. The results showed that the cocaine-induced CPP persisted over two 
weeks in C57BL/6J mice, while it diminished within one week among BALB/cJ mice. After 24-hours of cocaine withdrawal, high levels of 
locomotion as well as low levels of social interaction and aggressive behavior were found in C57BL/6J mice, but no significant changes 
were found in BALB/cJ mice, indicating that cocaine-induced CPP persistence, locomotion and social behavior are not consistent between 


these two strains, and that overall C57BL/6J mice are more susceptible to cocaine than BALB/cJ mice at the tested doses. 


Keywords: Cocaine; Conditioned place preference; Locomotion; Social behavior; Withdrawal 


C57BL/6J and BALB/cJ mice are two widely used 
strains in biomedical research, especially in behavioral 
neurosciences. Because these two strains have different 
neurochemical and endocrinological substrates (Agmo et 
al, 1999; Scislowska-Czarnecka et al, 2004; Bach et al, 
2011; Kundakovic et al, 2013), they display different 
responding patterns on many behavioral tasks. For 
example, BALB/cJ mice are with less sociability (the 
tendency to seek social interaction) (Sankoorikal et al, 
2006; Brodkin, 2007; Moy et al, 2007) but higher levels 
of anxiety-like behaviors (Bouwknecht & Paylor, 2002; 
Priebe et al, 2005; Verleye et al, 2011). CS7BL/6J mice 
performs better in learning and memory tasks (Crawley 
et al, 1997; Van Dam et al, 2006; Shi et al, 2008). 
C57BL/6J mice exhibit approaching responses toward a 
novel BALB/cJ mice exhibit 
avoidance (Belzung & Berton, 1997; Belzung & Barreau, 
2000). Bardo et al (1996) previously claimed that high 
levels of novelty seeking are associated with an 


environment, while 


increased risk of drug abusing. Studies also indicated that 
C57BL/6J апа BALB/cJ mice differ in several aspects of 
drug abuse. For example, C57BL/6J mice are prone to 
cocaine self-administration (Deroche et al, 1997; Thom- 


Science Press 


sen & Caine, 2011) and display a cocaine-induced condi- 
tioned place preference (CPP), however, BALB/cJ mice 
fail to demonstrate place preference to cocaine at the 
previously tested doses (Belzung & Barreau, 2000). 
Additionally, morphine preferences of C57BL/6J and 
BALB/cJ mice are differentially affected by social group 
and isolation during the CPP test (Kennedy et al, 2012). 
Drug abuse is often associated with sociability, 
emotion and memory (Curtis & Wang, 2007; Perrine et al, 
2008; Niigaki et al, 2010). Given the differences in these 
aspects between C57BL/6J and BALB/cJ mice, the 
comparisons of drug abuse between the two strains may 
provide more background information. CPP is a widely 
used paradigm in studying the rewarding effects of drugs 
and modeling some aspects of long term drug-seeking 
and relapse (Schechter & Calcagnetti, 1993; Sakoori & 
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Murphy, 2005). The drug paired cues play a critical role in 
the reinstatement of drug-seeking behaviors after a period 
of abstinence (Crombag et al, 2008; Su et al, 2013). 
Although differences in the susceptibility to the reinfor- 
cing properties of cocaine, morphine, and ethanol have 
been described among inbred mice strains (Cunningham et 
al, 1992; Eisener-Dorman et al, 2011), the persistence of 
CPP during cocaine withdrawal in C57BL/6J and BALB/ 
cJ mice remains unclear. Zhang et al previously (2002) 
proposed that the persistence of cocaine-induced CPP is 
study 
investigated: (1) the phenotypic variability between these 


strain dependent. Accordingly, the present 
two mice strains in cocaine-seeking behaviors reflected by 
the persistence of cocaine-induced CPP; (2) the differe- 
nces in locomotion and social behaviors after the 24-hour 


cocaine withdrawal. 


MATERIALS AND METHODS 


Animals 

Male C57BL/6J and BALB/cJ mice at 8-week of 
age were obtained from the Laboratory Animal Center of 
ХГап Jiaotong University (Shaanxi, China). The animals 
were housed in groups of four in standard transparent 
Makrolon cages (42 cmx26 cmx20 cm, lengthxwidthx 
height). The colony room was illuminated on a 12:12 
light-dark cycle (lights on 2000h) and the temperature 
was maintained at 2342 °C. Food and water were 





available ad libitum. Mice were allowed to adapt to 
housing conditions for one week and were handled daily 
by the same technician for three days prior to testing. AII 
protocols and procedures were approved by the Animal 
Care and Use Committee of Shaanxi Normal University. 


Chemicals and injection 
Cocaine-hydrochloride (Northwest Pharmaceutical 
Co., Ltd. Sinopharm, Xi’an, China) was diluted in saline 


(0.9% NaCl) and was administrated subcutaneously (s.c. 


20 mg/kg) (Zhang et al, 2002; Eisener-Dorman et al, 2011). 


Conditioned place preference test 

The place preference apparatus consisted of two 
large compartments (34 cmx25 cmx 32 cm, lengthx 
widthxheight) with different visual cues (one had gray 
walls and the other had white-black striped walls) separ- 
ated by a small middle compartment (11 cmx25 cmx32 
cm, lengthxwidthxheight). The middle compartment was 
an acclimation chamber with a door (7 cmx9 cm, 
height* width) in the center of the base. 
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Pre-test: On the day prior to conditioning, all 
animals were tested to determine any individual innate 
preference to either of the large lateral chambers. The 
mice (C57BL/6J, n-12; BALB/cJ, n-12) were given free 
access to each cue-decorated chamber when received a 
subcutaneous injection of physiological saline. 
Following 10-minute acclimation, the time spent in two 
lateral chambers was recorded for 15 minutes by a 
camera (Sony, HDR-XR260E) mounted 70 cm above the 
arena. After each trial, the chamber was thoroughly 
cleaned with 70% ethanol solution. 

Conditioning with cocaine: Each strain of mice 
(C57BL/6J, n=8; BALB/cJ, n=9) were conditioned with 
cocaine. Both cocaine and equivalent volume of 
physiological saline injections were given on the same 
day for four consecutive days. Specifically, in the 
morning, subjects were placed in one of the outer 
chambers with cocaine injections and were placed in the 
opposite chamber with saline injections at afternoon. 
Two injections per day in an alternating counter balanced 
sequence for four days, thus, providing four associative 
pairings for cocaine and saline. Mice were conditioned 
for 2 h after injections. The morning session and the 
afternoon session were at least 6 hr apart to allow time 
for cocaine clearance (Thiel et al, 2008). 

Post-test: CPP testing was conducted 24 hour after 
the last conditioning trial. Mice in a drug-free and a 
saline-free state were allowed to free access to each 
compartment. The time spent in two lateral chambers 
was recorded for 15 minutes. 

The persistence of CPP: Mice were housed in their 
home cages after the CPP testing. The day of post-test 
was taken as intermission day zero and then the place 
preference testing was conducted on intermission day 7, 
15 in BALB/cJ and on day 7, 15, 25 in C57BL/6J mice, 


respectively. 


Open-field test 

Another group of C57BL/6J and BALB/cJ mice was 
used in open-field test and social interaction test. Mice 
were randomly assigned into cocaine-treated groups (77-8) 
and control groups (7-8). The cocaine-treated groups (CC: 
C57BL/6J mice treated with cocaine, BC: BALB/cJ mice 
treated with cocaine) were administrated with cocaine (20 
mg/kg) at 0900h for four consecutive days. Control groups 
(CS: C57BL/6J mice treated with saline, BS: BALB/cJ 
mice treated with saline) were administrated with saline 
instead. Locomotion and anxiety-like behaviors were asse- 
ssed in an open field chamber 24 hr after the last injection. 
The open field chamber (50 cmx50 cmx25 cm, lengthx 
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widthxheight) made of white glacial polyvinyl chloride 
was brightly and evenly illuminated by four 60 W lamps 
mounted 1.5 m above (400 lux in the center of arena). The 
square arena was divided into 16 quadrants (four central 
and 12 peripheral) (Fiore & Ratti, 2007). Mouse was 
placed in the center of the arena and allowed to explore for 
5 minute. The time spent in the central and peripheral 
zones and numbers of crossings between quadrants were 
recorded. The anxiety-like behaviors were assessed by the 
time spent in the center of the arena and the locomotion 
was determined by the numbers of crossings. 


Same-sex social interaction test 

The social interaction test was conducted between 
1500h and 1700h. To eliminate possible influences from 
sexually motivated behaviors, only male-male dyads 
were used. The stimulus mouse was an unfamiliar, 
sexually naive individual that was approximately of the 
same age and size as the tested mouse. Testing were 
conducted in a neutral plastic cage (44 cmx22 cmx16 cm, 
lengthxwidthxheight) with wood shavings bedding (2 cm) 
and a removable opaque divider in the middle. The 
stimulus and tested mouse was confined in each side of 
the cage for 3 minute, then the divider was removed and 
the activities of the mice were recorded for 15 min by a 
video-recorder mounted 70 cm above. 

Mice behaviors were classified as investigatory behav- 
lor (sniffing face, body or anogenital area), aggressive beha- 
vior (pouncing, i.e. jumps or lunges; fighting, 1.e. tumbling 
or biting; chasing), body contact (staying together with 
another mouse or amicable grooming); self-grooming (cep- 
halocaudal progression that begins with rhythmic moveme- 
nts of the paws around the mouth and face, ears, descending 
to the ventrum, flank, anogenital area and tail) and other 
behaviors (digging, jumping, climbing the cage and resting). 


Statistical analysis 

All behavioral variables were scored from video 
footage according to established definitions by a naive 
observer using Observer 5.0 (Noldus, Netherlands). 
Statistical analyses were carried out via SPSS 10.0 
(SPSS Inc., Chicago, Illinois, USA). Data were checked 
for normality using the one-sample Kolmogorov- 
Smirnov test. The expression of place conditioning was 
analyzed using repeated-measures, with the intermission 
day as a repeated measure. Paired-samples t-test was 
used to evaluate the differences in time spending during 
the pre-test and intermittent tests. Data from the open- 
field test and the social interaction test were compared 
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using two-way ANOVA with strains and cocaine 
treatment as factors. Group differences were compared 
using post-hoc test except aggressive behavior which 
was abnormally distributed and was compared using 
Mann-Whitney U-test. All data were expressed as 





mean-SE. Statistical significance was taken at P< 
0.05. 


RESULTS 


Conditioned place preference 

The pre-test indicated that neither C57BL/6J 
(ta1y71.76, Р=0.106) nor BALB/cJ mice (ї(1)=0.484, 
P=0.638) showed preferences to either of the chambers 
(Figure 1A). Post-testing and intermittent testing 
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Figure 1 Time spent in saline- or cocaine-paired compartments 
by C57BL/6J and BALB/cJ mice 

A: Pre-test; В and C: Post-test and intermittent test in C57BL/6J and 

BALB/cJ mice following cocaine conditioning, respectively; The day of 


post-test was taken as intermission day 0; *: Р<0.05; **: Р<0.01. 
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indicated that cocaine could induce place preference in 
both mice strains (C57BL/6J, Аз в)=217.381, P<0.001; 
BALB/cJ, Ео,24)-7.101, P<0.001). However, a signific- 
ant interaction between cocaine conditioning and interm- 
ission day was found only іп BALB/cJ mice (Fo, = 
7.239, P=0.003) but not in C57BL/6J mice (Ғҙ2зу-1.043, 
Р-0.384). 

A significant preference to the cocaine-paired 
compartment was shown in both mice strains on the day 
of post-testing (intermission day 0) (C57BL/6J: t= 
3.849, Р-0.006; BALB/cJ: %у-4.092, P=0.03) (Figure 
ІВ, C) and persisted to intermission day 7 (47)-3.304, 
P=0.013) and 15 (¢7)=3.862, P=0.006), but not on day 25 
(1.985, Р-0.094) іп C57BL/6J mice (Figure 1B). 
However, this preference persistence was not found in 
BALB/cJ mice on either day (day 7: (10.779, P=0.458; 
day 15: (0.638, P=0.054) (Figure ІС). 


Open field behavior 

The main effect of strain was significant in the time 
spent in the central area (Fg, 294.736, P=0.038) and 
transitions (Ёз, 28)=28.776, Р < 0.001). However, the 
interactions between strain and cocaine administration 
had no effects on either the total transitions (F6, 28)=2.028, 
P=0.165) or the time spent in the central area (Fo, 
23.461, Р-0.073). 

Although no differences were found in the time 
spent in the central area between С57ВІ/6) and 
BALB/cJ control mice (Mean difference=1.158, P= 
0.996), C57BL/6J control mice showed a higher level of 
locomotor activity (total transitions) (Mean differe- 
nce=46.000, P=0.02). Cocaine-administrated C57BL/6J 


** 


























ЗО гд cs 
| CC 
25} E BS 
ШЕШ BC 





пина 
КА 
MANN 


ОКК] 


Time spend in the center (sec)/5 min 














Mice groups 


mice showed a greater number of total transitions than 
both the C57BL/6J control mice (Mean difference= 
35.375, P=0.04) and the cocaine-administrated BALB/cJ 
mice (Mean difference=79.250, P << 0.001), as well as 
spent more time in the central area (Mean diffe- 
rence=14.799, P=0.008) than the cocaine-administrated 
BALB/cJ mice. However, cocaine-administration did not 
affect either the total transition (Fg, 28)=2.580, Р-0.119) 
or the time spent in the central area (F6, 20.178, 
P=0.677) in BALB/cJ mice (Figure 2). 


Same-sex social interaction 

The male-male interactions indicated that strain and 
cocaine both significantly affected the social investi- 
gations (strain [duration: F5 зу-18.216, P<0.001, freque- 
ney: Ға, 28)=107.113, Р < 0.001]; cocaine [duration: 
Ев28)-8.817, P=0.006, frequency: F55-20.502, Р < 
0.001], the contact behaviors (strain  [duration: 
F(358y-30.856, P<0.001, frequency: Ка) -F(35sj- 356.608, 
P<0.001]; cocaine [duration: Аз 28)=8.891, P=0.006, 
frequency: Fig.9= Fa, 230.888, P=0.354]) and the 
frequencies of self-grooming (strain: Рззу=66.336, P< 
0.001; cocaine: Аз 28)=14.131, Р<0.001). The frequen- 
cies of aggressive behaviors between two mice strains 
were significantly different (F6, 28)=33.618, P<0.001). 
Moreover, significant interactions between strain and 
cocaine were found in the social investigation (duration: 
Ға, 28)=15.386, P=0.001; frequency: F6, 2у-16.084, P 
< 0.001), the self-grooming (Duration: Fe, 914.598, 
P=0.001; Frequency: Fg, 29=25.121, P<0.001), duration of 
contact behavior (Fg, 29=9.675, Р-0.004) and aggressive 
behaviors (Fa, 294.423, Р-0.045). 
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Figure 2 Open-field behavior of C57BL/6J and BALB/cJ mice after the 24 h cocaine withdrawal 
A: The time spent in the central area; B: Total transitions; *: Р<0.05; **: Р<0.01; CS and CC: C57BL/6J mice administrated with saline and cocaine, 


respectively; BS and BC: BALB/cJ mice administrated with saline and cocaine, respectively. 
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Figure 3 Duration and frequency of same-sex social interactions in C57BL/6J and BALB/cJ mice after the 24 h of cocaine 
withdrawal 
A: Investigation; B: Aggression; C: Body contact; D: Self-grooming; *;P<0.05, ““Р<<0.01; CS and CC: C57BL/6J mice administrated with saline and cocaine, 


respectively; BS and BC: BALB/cJ mice administrated with saline and cocaine, respectively. 


C57BL/6J control mice showed high levels of the P «0.001; frequency: mean difference=27.750, P< 
social investigations (duration: mean difference=39.646, 0.001), the aggressive behaviors (duration: U-3; 
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P=0.001; frequency: U-1, P<0.001) and the contact 
behavior (duration: mean difference=93.608, P=0.000; 
difference-30.005, Р < 0.001) 
compared with the BALB/cJ control mice. Moreover, the 


frequency: mean 
levels of the self-grooming behaviors between these two 
control mice groups were also significantly different 
(duration: mean difference- 17.296, P=0.007; frequency: 
mean difference— 1.875, Р< 0.001) (Figure 3). 

Cocaine-administrated C57BL/6J mice exhibited 
attenuations in the social investigations (duration: 
mean difference=—33.360, P «0.001; frequency: mean 
difference=—16.500, P < 0.001), the self-grooming 
behaviors (duration: mean difference——23.524, P «0.001; 
frequency: mean аШегепсе--5.250, Р < 0.001), the 
contact behaviors (duration: mean difference=—65.813, 
P=0.001; frequency: mean difference=—2.875, P=0.230) 
and the aggressive behaviors (duration: U=9, P=0.015; 
frequency: U=22, P=0.328) than those of their control 
mice (Figure 3). 

Cocaine withdrawal had no apparent effect on social 
investigations (duration: mean difference=4.6140, 
P=0.506; frequency: mean difference=—1.000, P=0.717), 
U=31.5, P=0.927; 
contact 


aggressive behaviors (duration: 
U=31.5, Р-0.927), 
(duration: mean difference=1.390, P=0.928; frequency: 
mean difference=0.250, P=0.916) or self-grooming 
behaviors (duration: mean difference=8.401, P=0.166; 
frequency: mean difference=0.750, P=0.383) in BALB/ 
cJ mice (Figure 3). 


frequency: behaviors 


Compared with cocaine-administrated BALB/cJ 
mice, cocaine-administrated C57BL/6J mice were more 
engaged in social investigations (duration: mean 
difference=1.6725, P=0.809; frequency: mean differe- 
nce=12.250, P<0.001), aggressive behaviors (duration: 
U=15, P=0.083; frequency: U=6.5, P=0.005), contact 
behaviors (duration: mean difference=26.495; P=0.095; 
frequency: mean difference=32.625, P<0.001), but less 
in self-grooming behaviors (duration: mean difference=— 
14.630; P=0.02; frequency: Mean difference =0.875; P= 
0.35) (Figure 3). 


DISCUSSION 


Persistence of cocaine-induced place preference 
Although Belzung & Barreau (2000) claimed that 
cocaine may be not able to induce place preference in 
BALB/cJ mice, in this study, after cocaine conditioning, 
both C57BL/6J and BALB/cJ mice showed a significant 
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preference to the cocaine-paired compartment, indicating 
that cocaine induces rewarding effects in both mice 
strains (Miner, 1997; Zhang et al, 2002; Eisener-Dorman 
et al, 2011). Since the rewarding effects of cocaine could 
be influenced by administration patterns (Zhang et al, 
2002), this discrepancy may be due to the fact that 
Belzung & Barreau (2000) chose a dose of 10 mg/kg and 
a conditioning trial of half-hour, instead of 20 mg/kg and 
two-hour, respectively, as in this present study. These 
phenomena indicate that the CPP establishment could be 
affected by both the doses of cocaine and the duration of 
reinforcement. 

In this study, the cocaine-associated preference 
diminished within one week after the withdrawal in 
BALB/cJ mice, but persisted at least two weeks after the 
withdrawal in C57BL/6J mice, which is consistent with 
previous studies indicating that the cocaine-induced CPP 
could be maintained by repeated testing two or four 
weeks after conditioning (Mueller & Stewart, 2000; 
Zhang et al, 2002). Tran-Nguyen et al (1998) found that 
the cocaine-seeking behaviors were getting more intense 
during the course of cocaine withdrawal in rats. Su et al 
(2013) reported that the cocaine-induced CPP could 
remain viable at three weeks of the withdrawal. Mueller 
& Stewart (2000) found that in animals tested only once, 
a spontaneous reduction in CPP was shown at six weeks 
after conditioning. These findings suggest that CPP may 
be variably maintained or extinguished dependent upon 
the timing and frequency of testing following 
conditioning (Sakoori & Murphy, 2005). Intermittent re- 
testing for CPP itself may act as a secondary reinforce 
that strengthens the association between drug experience 
and environment and may actually help maintain 
cocaine-induced CPP (Mueller & Stewart, 2000; Sakoori 
& Murphy, 2005). However, no such reinforcement was 
shown in BALB/cJ mice. The genetic variations in 
pharmacokinetic may not account for the differences in 
cocaine responsiveness observed because no difference 
in the incorporation of [*H]-cocaine has been found 
between C57BL/6J and BALB/cJ mice (Seale, 1991). An 
alternative explanation of the differences ш the 
persistence of CPP expression is the poor capacity of 
learning and memory of BALB/cJ mice (Crawley, 2000; 
Shi et al, 2008) because the CPP task involves learning 
processes (Fleming et al, 1994; Thiel et al, 2008). 
Similarly, Oler & Markus (1998) found that young rats 
showed a stronger retention for conditioning context 
after conditioning than aged rats. Thus, the persistence of 
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CPP expression indicates that the development of 
ordinary memory and addictive memory between 
C57BL/6J and BALB/cJ mice are quite different. 

Studies suggested that place conditioning may be 
used to model some aspects of long term drug-seeking, 
drug relapse, and the hedonic properties of drugs 
(Sakoori & Murphy, 2005). Our results indicate that the 
persistence of cocaine-induced CPP during prolonged 
withdrawal may contribute to the sustained vulnerability 
to “relapse” of cocaine-seeking behavior. Moreover, 
compared with C57BL/6J mice, BALB/cJ mice are 
characterized by higher levels of anxiety-like behaviors, 
lower locomotor activities and are less sociability 
(Crawley et al, 1997; Sankoorikal et al, 2006; Brodkin, 
2007; Moy et al, 2007). The different CPP persistence 
showed in this study between the two mice strains is in 
accordance with the previous study indicating that the 
highly social species have different susceptibility to the 
effects of drug compared with the less social species 
(Curtis & Wang, 2007). 


Cocaine-induced locomotor activities 
In this present study, the 24-hour 
withdrawal had no effect on the anxiety-like behaviors of 


cocaine 


the two mice strains, which is consistent with previous 
reports (Niigaki et al, 2010; Stoker & Markou, 2011). 
However, in rats, cocaine withdrawal is associated with 
increased anxiety-like behaviors in the elevated plus 
maze (Perrine et al, 2008; Hall et al, 2010). Moreover, 
the anxiogenic effects of abstinence from cocaine are 
also correlated with different protocols (e.g. the elevated 
plus maze, the open field or the light-dark box) (Stoker 
& Markou, 2011; de Oliveira Cité Mdo et al, 2012). 

In this study, after a 24-hour cocaine withdrawal, 
the level of locomotor activities of C57BL/6J mice 
increased while that of BALB/cJ mice was maintained at 
the same level. These results are consistent with previous 
studies indicating that a marked increase or dose- 
dependent stimulant effects on locomotor activities in 
response to cocaine have been found in C57BL/6J mice 
(Zhang et al, 2002; Eisener-Dorman et al, 2011; 
Thomsen & Caine, 2011). BALB/cJ mice are less 
sensitive to the cocaine-associated stimulant effects on 
locomotor activities (Ito et al, 2007; Eisener-Dorman et 
al, 2011; Thomsen & Caine, 2011), though some reports 
claimed that BALB/cJ mice show hyperlocomotor 
activities in response to cocaine (Miner, 1997; Kuzmin et 
al, 2000). In the open field test, the increased locomotion 
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activities were found in rats after 24-hour of abstinence 
from cocaine (de Oliveira Cité Mdo et al, 2012). 
However, some studies showed that cocaine withdrawal 
suppresses locomotor activities in rats (Baldo et al, 1999; 
Koeltzow & White, 2003). One potential explanation for 
these discrepancies is that the locomotor activity may be 
influenced by the measurement environment (e.g. the 
open field, the conditioning chamber or the light-dark 
box) or cocaine dosage, administration pattern and 
withdrawal time (Zhang et al, 2002; Niigaki et al, 2010; 
Eisener-Dorman et al, 2011; Stoker & Markou, 2011; de 
Oliveira Cito Mdo et al, 2012). For example, the context 
in which a drug is experienced can significantly 
influence both acute and sensitized responses to the drug 
(Badiani & Robinson, 2004; Eisener-Dorman et al, 2011). 
A specific interaction between the cocaine and the 
environment may result in context-dependent sensitiz- 
ation in BALB/c mice (Eisener-Dorman et al, 2011). The 
present results indicate that cocaine withdrawal induces 
locomotion changes in C57BL6/J mice. 


Cocaine-induced social behaviors 

In this present study, when interacting with same- 
sex individuals, C57BL/6J mice were more active in 
social investigation, body contact and aggression than 
those of BALB/cJ mice, which are consistent with 
previous reports (Sankoorikal et al, 2006; Brodkin, 2007; 
An et al, 2011). Fairless et al (2008) demonstrated that 
the size of the corpus callosum relative to brain weight is 
associated with sociability among these two mice strains. 
Moreover, compared with the saline control, cocaine 
withdrawal induced remarkable decreases in social 
investigation, contact behavior and aggressive behavior 
in C57BL/6J mice as described in previous reports 
(Rademacher et al, 2002; Estelles et al, 2007). Although 
cocaine induces complicated changes in social behaviors, 
no agreement has been reached on its specific effects on 
aggression (Moeller et al, 1997; Dhossche, 1999). Some 
studies also indicated that other than strains or species, 
the aggressive behaviors may also be affected by the 
patterns of drug administration (single or binge 
administration), the dosing regimens and the specific 
temporal window assessed (Estelles et al, 2004, 2007; 
Wang et al, 2012). 

Interestingly, no significant effects of cocaine 
withdrawal on either the social behaviors or the locom- 
otor activities were found in BALB/cJ mice, indicating 
that there may be a dissociation between cocaine-induced 
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CPP and locomotion or social behavior in BALB/cJ mice. 
Compared with C57BL/6J mice, different mechanisms 
may underlie such effects in BALB/cJ mice. Combined 
the persistence of CPP test, these results indicate that 
С57В1/6] mice are more 
withdrawal than BALB/cJ mice. However, because 


susceptible to cocaine 


BALB/cJ mice are characterized with low levels of 
social investigation and contact behavior, it is possible 
that what we have observed in this study is only a ‘floor 
effect’ of cocaine on the two behaviors. 

Mesolimbic dopamine (DA) is responsible for 
cocaine-induced behavior and locomotor activation 
(Sarnyai, 1993; Tran-Nguyen et al, 1998). Hyperloc- 
omotion induced by psychostimulants is mediated by the 
mesolimbic dopaminergic system, whereas stereotyped 
behaviors are mediated by the nigrostriatal dopaminergic 
system (Ito et al, 2007). The balances of the activation of 
dopaminergic neurons between mesolimbic апа 
nigrostriatal systems may play an important role to 
engender corresponding behavioral outcomes (Ito et al, 
2007). The differences between C57BL/6J and BALB/cJ 
mice in dopaminergic function within the prefrontal 
cortex and the striatum have been reported (Hervé et al, 


1979; Helmeste & Seeman, 1982). Thus, the differences 
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Abstract: Transmissible spongiform encephalopathies (TSEs), also known as prion diseases, are a group of fatal neurodegenerative 
diseases detected in a wide range of mammalian species. The “protein-only” hypothesis of TSE suggests that prions are transmissible 
particles devoid of nucleic acid and the primary pathogenic event is thought to be the conversion of cellular prion protein (PrP^) into 
the disease-associated isoform (PrP?^). According to susceptibility to TSEs, animals can be classified into susceptible species and low 
susceptibility species. In this review we focus on several species with low susceptibility to TSEs: dogs, rabbits, horses and buffaloes. 
We summarize recent studies into the characteristics of low susceptibility regarding protein structure, and biochemical and genetic 


properties. 
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Transmissible spongiform encephalopathy (TSE), 
or prion disease, is an invariably fatal neurodegener- 
ative disease detected in a wide range of mammalian 
species, including Scrapie in goats (Capra hircus) and 
sheep (Ovis aries); bovine spongiform encephalopathy 
(BSE) in cattle (Bos taurus); chronic wasting disease 
(CWD) in elaphure (Elaphurus davidianus) and moose 
(Alces americanus); feline spongiform encephalopathy 
(FSE) in cats (Felis catus); transmissible mink 
encephalopathy (TME) in minks (Mustla vison); and 
Creutzfeldt-Jakob disease (CJD), variant Creutzfeldt- 
Jakob disease (vCJD), fatal familial insomnia (FFI), 
Gerstmann-Straussler-Scheinker syndrome (GSS) and 
Kuru in humans (Homo sapiens) (Collins et al, 2004; 
Prusiner, 1982). Humans and other animals infected 
with TSE are clinically and pathologically charac- 
terized with neuronal progressive vacuolation, stellate 
cell gliosis, spongiform lesions in gray matter, 
amyloid deposition and eventually disastrous degener- 
ation and death (Prusiner, 1998). No effective treatm- 
ents have been found and the World Health Organ- 
ization has named TSE and AIDS as two major health 
problems of the 215 century. 
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The “protein-only” hypothesis of TSE suggests that 
the pathogenic factors of TSE are not bacteria or a virus, 
but a protein devoid of nucleic acid, which has been 
named prion protein (PrP). PrP is encoded by the prion 
protein gene (PRNP) (Prusiner, 1982). Normal cellular PrP 
(PrP^) expresses in the cells of mammalian species and 
the number of amino acid residues varies from 253 to 
264 across species (Wopfner et al, 1999) and are all 
highly conserved (Figure 1). PrP has two signal peptide 
sequences, a N-terminal and a C-terminal. Mature PrP 
has an intra-molecular disulfide bond and two 
glycosylation sites, and 15 anchored on the cell membr- 
ane surface via glycosylphosphatidylinositol (ОРГ) at the 
C-terminal (Aguzzi et al, 2008). According to the 
protein-only hypothesis TSE is a conformational disease 
and under certain circumstances cellular PrP“ mistakenly 
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Figure 1 Amino acid sequences of PrP in 16 mammals (data from GenBank) 
Human (NM 000311.3), chimpanzee (NM 001009093.3), Rhesus (NM 001047152.1), deer (AY330343.1), elk (EU082291.1), mouse (NM 011170.3), rat 


(NM 0126312), pig (NM. 001008687.1), 


sheep (NM 001009481.1), goat(JF729302.1), rabbit (NM 001082021.1), dog (NM 001013423.1), cat 


(EU341499.1), horse(NM_001143798.1), cattle (NM. 181015.2), buffalo ( KC.1). 137634. 
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converts into the disease-associated isoform (PrP*). 
Although the primary structures of PrP* and PrP* are the 
same, their secondary structures are quite different. PrP“ is 
enriched with a-helix (42% are a-helix, 3% are p-fold), 
whereas PrP is enriched with f-fold (43% are f-fold, 
30% are a-helix) and is protease resistant (McKinley et al, 
1983; Pan et al, 1993; Prusiner, 1982). The massive 
intracellular accumulation of PrP® induces formations of 
oligomer and amyloid fibrils, and eventually neuronal 
degeneration (Barron et al, 2007; Caughey et al, 2009). 
PrP plays vital roles in the pathological process of TSE. 
Knockout and low expression of PrP effectively abolishes 
or reduces susceptibility to TSEs, respectively (Brandner 
et al, 1996; Büeler et al, 1993), whereas high expression is 
associated with susceptibility and a shortened incubation 
time for disease development (Manson et al, 1994). 

TSE susceptibility is species-specific. Previous stu- 
dies show high susceptibility in hamsters (Mesocricetus 
auratus) to TSE, as they can be infected by various PrP*° 
virus strains isolated from human, cattle, goats, mice 
(Mus musculus) and minks (Bessen & Marsh, 1992; 
Gibbs & Gajdusek, 1973; Kimberlin & Walker, 1977; 
Thomzig et al, 2006). Similar high susceptibility is also 
found in mice (Chandler, 1961; Gibbs & Gajdusek, 1973; 
Hill et al, 2000; Lasmézas et al, 1997; Thomzig et al, 
2006). However, rabbits could not be infected by PrP“ 
strains isolated from human, goats and mice (Barlow & 
Rennie, 1976; Gibbs & Gajdusek, 1973). During the 
outbreak of BSE in the UK, infections in humans and 
several species of feline were reported, but no infection 
was found in dogs (Canis familiaris) or horses (Equus 
caballus) (Aldhous, 1990; Kirkwood & Cunningham, 
1994). Collectively, species with confirmed susceptibility 
to TSE include humans, rhesus monkeys (Macaca 
mulatta), hamsters, mice, minks, elaphures, moose, goats, 
sheep, cattle and raccoons (Procyon lotor) (Imran & 
Mahmood, 2011). Only a few species, such as dogs 
(Canis familiaris), rabbits (Oryctolagus cuniculus) and 
horses (Equus caballus) have been recognized as TSE 
resistant (Fernandez-Funez et al, 2011; Yuan et al, 2013; 
Zhang, 2011а). Interestingly, although more than 
190,000 cattle were infected by BSE, and buffaloes 
(Bubalus bubalis) and cattle are closely related, no 
buffalo has been reported with BSE 
(http://www.oie.int) and are of low susceptibility to BSE 
(Zhao et al, 2012). In this review, based on TSE 
susceptibility, animals have been classified into TSE 


infection 


susceptible animals and TSE low susceptible animals. 
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The pathological mechanisms of TSE are yet to be 
clarified. Although highly susceptible animals are 
important to understanding this disease, studies on anim. 
als with low susceptibility provide a new angle from 
which to examine TSE. Here, we review recent research 
developments on protein structures, biochemical charac- 
teristics and genetic features of four animals (dogs, rabb- 
its, horses and buffaloes) with low susceptibly to TSE. 


Dogs 

During the outbreak of BSE in the UK, several 
species of feline were reportedly infected, including 
cheetahs (Acinonyx jubatus), pumas (Pumaconcolor) and 
cats (Kirkwood & Cunningham, 1994). Since 1990, 
about 100 cats and 29 captive felines, including 15 
cheetahs, four lions (Panthera leo), three leopard cats 
(Prionailurus bengalensis), three pumas, three tigers 
(Panthera tigris) and one Asian golden cat (Catopuma 
temminckii), have been diagnosed with FSE (Imran & 
Mahmood, 2011). The presumed infection source was 
PrP**-contaminated food; however, dogs and cats аге 
provided similar food and no dogs were reported with 
TSE (Imran & Mahmood, 2011; Kirkwood & Cunni- 
ngham, 1994; Wopfner et al, 1999). With further laborat- 
ory cell experiments, dogs have been recognized as a 
species with low TSE susceptibility. For example, when 
Madin-Daby canine kidney cells (MDCK) were infected 
with brain tissue homogenates from CJD patients or RML 
prion strain isolated from scrapie animals, although the 
biosynthesis and processing of PrP^ in MDCK are similar 
with those in М2аРКІ cells of murine neuroblastoma, 
which are highly susceptible to TSE, no PrP* was found 
in MDCK. When infected MDCK were used to infect 
N2aPK 1 cells, no PrP% was found in N2aPK1 cells either 
(Ploymenidou et al, 2008; Zhang & Liu, 2011). 

The gene polymorphism of PRNP is correlated with 
TSE susceptibility (Westaway et al, 1994). In humans, at 
least 30 mutations of PRNP are intertwined with TSE 
susceptibility (Lloyd et al, 2011). In dogs, the amino acid 
residue 187 and 229 of the PrP sequence are histidine 
and glycine, respectively, whereas, they both are arginine 
in cats (Wopfner et al, 1999). No FSE-related polym- 
orphic site was found by screening encoding sequences 
of PRNP in 609 animals (including 15 FSE infected 
cases) and 29 species from 22 genera of the Order 
Carnivora, but Stewart et al (2012) did notice that amino 
acid residue 163 in all canines is either aspartate or 
glutamic acid, indicating this locus may have some 
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connection to TSE susceptibility. 

The three-dimensional structure of PrP“ may be 
another tool in resolving the puzzle of TSE susceptibility 
(Lin & Wen, 2011). To understand the structural 
differences of PrP^ in animals with low and high 
susceptibility to TSE, Lysek et al (2005) carried out a 
study on nuclear magnetic resonance (NMR) structures 
of PrP“ in dogs (canine PrP, cPrP), cats (feline PrP, 
fPrP), pigs (sus scrofa PrP, scPrP) and goats (ovine PrP, 
ovPrP). Their overall three-dimensional structures are quite 
close, consisting of a N-terminal (constituted of about 100 
amino acid residues in random coil) and a globular domain 
in the C-terminal (including three o-helixes and a pair of 
short, reverse paralleled B-folds constituting about 100 
amino acid residues). The globular domain in the C- 
terminal is species-specific, e.g., four amino acids 
(Aspl59Asn, Arg177His, Lys185 Arg and Gly229Arg) аге 
different between cPrP^ and fPrP^; Asp-159 and Arg-177 in 
cPrP^ make it unique in potential distribution; in fPrP^, 
scPrP^ and ovPrP^ same positive potential distribution 
patterns are observed in their C-terminals. 

The conversion of PrP“ into PrP* is critical in TSE 
pathogenesis. Besides PrP“ and Ргр“, the existence of 
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the third state-B intermediate state, consisting mainly of 
B-fold in circular dichroism (CD) (Hornemann & 
Glockshuber, 1998), has raised attention due to its 
capability in introducing TSE (Collinge & Clarke, 2007). 
Khan et al (2010) claimed that in different species, the 
propensities in forming the B intermediate state are one 
of the vital factors influencing TSE susceptibility. Using 
dual wavelength CD, Khan et al (2010) compared the 
structures of the globular domain in high-susceptible 
(hamsters and mice) and low-susceptible species to TSE 
(rabbits, horses and dogs) and found that under inducing 
conditions with different pH values and urea concentr- 
ations, the propensities of forming the B intermediate 
state vary in different species. At pH 7.0, urea 
concentrations have no effect on PrP and no В interm- 
ediate state is observed in all five species; at pH 5.0, the 
structure of hamster PrP“ (hPrP‘) is most unstable and 
easily forms the В intermediate state; at pH 4.0, the PrP* 
in all five species is unstable and easily forms the f 
intermediate state, the lowest concentration of p inter- 
mediate state occurs in dogs (Figure 2). The propensities 
in forming the В intermediate state (hamsters > mice > 
rabbits > horses > dogs) can also be adopted in evaluating 
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Figure 2 Propensities of conversions of PrP“ into the p intermediate state in different species at pH 7.0 (A), 5.0 (B), 4.5 (C) and 4.0 
(D) and different concentrations of urea (modified from Khan et al, 2010) 
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species’ susceptibilities to TSE (Fernandez-Funez et al, 
2011). Using molecular kinetic methods, Zhang & Liu 
(2011) found stable molecular structures of wild cPrP^ 
under both neutral and acidic pH conditions, and in the 
neutral condition the salt bridge between D177 and R163 
improves structural stability. These studies help to 
explain the mechanism of low TSE susceptibility in dogs 
and general TSE pathogenesis. 


Rabbits 


No cases of spontaneous TSE infection in rabbits 
have been reported to date. In 1973, Gibbs & Gajdusek 
failed to infect rabbits with either brain tissues from 
human CJD or Kuru patients or brain tissues from 
scrapie animals and minks with TME. In 1976, Barlow & 
Rennie failed to infect rabbits with ME7 strains isolated 
from animals with scrapie. By constructing rabbit PrP* 
(RaPrP‘) over-expressed murine neuroblastoma tumor 
cell lines, Vorberg et al (2003) confirmed that RaPrP* can 
neither be infected by RML strains nor convert into the 
РгР5°, The in vivo experiments conducted by Fernandez- 
Funez et al (2010) support the view that rabbits are TSE 
resistant species. Fernandez-Funez et al (2010) expressed 
full length PrP of hamsters, mice and rabbits in 
drosophilae voided of endogenous PrP, and found that 
cavernous transformation and isomers similar with PrP* 
can only be found in the brains of transgenic drosophilae 
expressing shPrP^ and mouse PrP“ (moPrP^), but not in 
drosophilae expressing RaPrP©. Moreover, Bellotti & 
Chiti (2008) reported that TSE correlates with the 
deposition of amyloid fibrils. Zhou et al (2011) found 
that Ficoll 70 and dextran 70 significantly accelerate the 
fibration of hPrP^ and bPrP*, but prevent fibration of 
RaPrP^; however, different from hPrP^ and bPrP^, 
RaPrP^ does not have fragments resisting protease K 
digestion. 

Initially it was presumed that one of the possible 
reasons rabbits are resistant to TSE infection is that 
certain transforming factors are lacking in their cellular 
environment or some inhibitory factors are expressed and 
the conversion of PrP“ to PrP^ is prevented. This 
presumption was later proved wrong. When PK13 cells, 
expressing ЅҺРІРС, moPrP© and vole (Microtus spp.) 
PrP^, respectively, were infected with shPrP*, moPrP* 
and vole PrP*, massive replications of Ргр“ were 
observed, indicating that factors necessary for PrP 
conversion exist in rabbit cells (Courageot et al, 2008; 
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Vilette et al, 2001). Rabbits and mice share 8796 
similarity in amino acid sequences (33 amino acids are 
different, including 22 in mature peptides). When amino 
acid residues 99, 108, 173 and 214 in moPrP^ were 
mutated into the corresponding amino acid residues in 
RaPrP (Asn99Gly, Leu108Met, Asn173Ser and Val214lle, 
respectively) 
neuroblastoma (MNB), and then RML prion strains were 
used to infect MNB, the mutants of moPrP* could not 
convert into PrP™, indicating that several amino acid 


and were overexpressed in mouse 


residues in RaPrP^ can prevent the replication of isomers 
of РгР5° (Vorberg et al, 2003). About 33% of the different 
amino acids in moPrP* and RaPrP^ locate around the 
attachment site of the GPI-anchor. Nisbet et al (2010) 
stably transfected RK13 cells voided of endogenous PrP* 
with a constructed double-mutant (Ser230Gly and 
Ser231Val) model of moPrP^, MoPrP-RbGPI, and then 
infected RK13 cells expressing MoPrP-RbGPI using 
human prion strains M1000 or MU-02 isolated from 
mice brain homogenates, The results showed that 
MoPrP-RbGPI could change neither the attachment of 
the GPI-anchor nor the location of PrP© on cells, but no 
PrP* produce either, indicating that rabbit-specific 
amino acids may interfere with PrP* and PrP“ contact 
and eventually prevent conversion of PrP^ to PrP* 
(Nisbet et al, 2010). These findings indicate that the 
resistance of rabbits to TSE may be attributable to their 
unique RaPrP^ structure (Lin & Wen, 2011). Wen et al 
(2010a) adopted NMR techniques to study the solution 
structure of RaPrP* and found that compared to hPrP^, 
mPrP© and bPrP®, RaPrP© features a unique charge 
distribution pattern. A large consecutive positive 
potential area exists on the protein surface of RaPrP© 
which may interfere to interact with molecules such as 
chaperon protein X, prevents the proliferation of PrP*° 
(Wen et al, 2010a). Khan et al (2010) found a critical 
helix-capping motif interacting with the third а-һейх and 
regulating the f-intermediate state by exploring the 
crystal structure of RaPrP*. As we mentioned earlier, the 
complexities of PrP“ forming the B-intermediate state in 
different species are correlated with their susceptibility to 
TSE. Compared with ShPrP* and moPrP^, RaPrP^ is 
difficult to transform into the B-intermediate state (Figure 
2) (Khan et al, 2010). The irregular curling fragment, 
called an a2-p2 loop, locates between the second p-fold 
and the second a-helix (165-172). The epitope consisting 
of the «2-82 loop and the C-terminal of the third a-helix 
is considered capable of recognizing protein X and 
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regulating progression of TSE (Kaneko et al, 1997). 
Protein dynamics analysis shows that RaPrP has a 
constructively highly ordered 82-02 loop (Wen et al, 
20102) which may function as a species barrier for TSE 
dissemination (Lin & Wen, 2011). However, compared 
with wild RaPrP^, S173N and I124V mutations affect the 
interactions of the B2-a2 loop with the third a-helix, and 
thereafter decrease the stability of the entire construct 
(Wen et al, 2010a, b). In addition, when the salt bridges 
between D202-R156 and D178-R164 were removed in 
hPrP^ and moPrP^, although secondary protein structures 
remained intact, the helix structures of RaPrP© were 
destroyed, indicating that salt bridges are important to 
the stability of RaPrP^ (Zhang, 2009, 2010, 201 1a). 

Recently, Joaquin Castilla’s research group has 
raised questions about the view that rabbits are resistant 
to TSE. They amplified rabbit brain homogenates using 
serial automated protein misfolding cyclic amplification 
(saPMCA) and then inoculated this in vitro novel PrP 
into the brains of three other rabbits. One rabbit was 
found with TSE symptoms 766 days after inoculation 
even although no exogenous PrP* was involved. Then 
the brain homogenates from this infected rabbit could 
100% infect RaPrP© over-expressed transgenic mice. 
Therefore, Chianini et al (2012) claims that rabbits are 
not TSE resistant. Furthermore, using saPMCA, when 
the amplified proteins from mixtures of rabbit brain 
homogenates and BSE prion strains were inoculated into 
the brains of RaPrP^ over-expressed transgenic mice, the 
resultant strains similar to BSE prion strains were 
discovered (Vidal et al, 2013). Fernandez-Borges et al 
(2012) claims that in vivo infective experiments are 
imperfect when forming the conclusion that rabbits are 
resistant to TSE, especially when supported only by the 
observation that rabbits can not be infected with TSE 
naturally. 


Horses 


As there is no reports of horses being naturally 
infected with TSE, horses are recognized as low 
susceptibility species (Zhang, 2011a). Relative to dogs 
and rabbits, fewer studies have looked at low suscept- 
ibility in horses. Studies on the conversion of PrP“ into 
the В intermediate state show that under unstable 
conditions at pH 4, the PrP“ of hamsters, mice, rabbits, 
dogs and horses can convert into the B intermediate state, 
but the lowest level of B intermediate state is found in 
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horses (Figure 2), indicating that equus caballus PrP* 
(ecPrP^) is relatively stable (Khan et al, 2010). Structural 
NMR on ecPrP© found two horse-specific amino acid 
alterations in its В2—02 loop (Ser-167 and Lys-173, 
respectively), among which, S167 affects the highly 
ordered solution structure of the В2—02 loop and may 
influence the low susceptibility of horses to TSE (Pérez 
et al, 2010). However, when amino acid residue 167 in 
moPrP^ was mutated from asparagine into ecPrP‘- 
specific serine (MoPrPP!^5), although МорРгРР!®? and 
ecPrP^ share similar NMR structures and their 82-02 
loops in solutions are both relatively highly ordered, 
spongiform lesions were found in mice expressing 
MoPrPP!65 and neural diseases can be induced with the 
accumulation of PrP* in the brain (Sigurdson et al, 2011). 
Therefore, the ordered state of the В2—02 loop in solution 
alone does not fully explain different susceptibilities to 
TSE (Lin & Wen, 2011). Moreover, as the salt bridge in 
RaPrP* stabilizes protein structures, similar salt bridges 
consisting of GLU196-ARG156-HIS187, ARGI56- 
ASP202 and GLU211-HIS177 are also found in ecPrP* 
(Zhang, 2011b). The structures of ecPrP^ and cPrP^ are 
stable under both neutral and acidic conditions (Zhang, 
2011a). A common phenomenon found among RaPrP^, 
ecPrP^ and cPrP^, is the salt bridge ASP177-ARG163 
connects with the В2—02 loop of PrP and is probably 
correlated with TSE susceptibility (Zhang, 2011a). 
Nevertheless, the low susceptibility of horses to TSE 
requires further work. 


Buffalo 


BSE was initially found in the UK in 1986, rapidly 
spread to over 25 countries, and caused major economic 
losses (Harman & Silva, 2009; Wells et al, 1987). BSE 
can also infect humans via the food chain and cause 
human vCJD (Collinge et al, 1996; Hill et al, 1997). The 
multiple pathogenic pathways of TSE, including spont- 
aneous mutant, inheritance and infection (Nicholson et al, 
2008), may explain why even after meat and bone meal 
was strictly forbidden, more than 15 000 BSE infected 
cattle were found in the UK (http://www.oie.int). Worl- 
dwide, there were over 190 000 taurus cattle, 1 Bos 
indicus, and 1 Bos indicus х Bos taurus cross reported 
with BSE infections (data of OIE, Novakofski et al, 2005; 
Seuberlich et al, 2006). Although buffaloes and cattle are 
quite close phylogenetically, no case of BSE infected 
buffalo was ever reported, suggesting that genetic factors 
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are crucial to BSE susceptibility (Zhao et al, 2012). The 
expression level of PrP is closely correlated with BSE 
susceptibility. Studies show that the PRNP gene of cattle has 
two indel (insertion and deletion) polymorphisms (a 23-bp 
indel in putative promoter, and a 12-bp indel in intron 1). 
These Indel polymorphisms affect gene expression (Msalya 
et al, 2011; Sander et al, 2005) and eventually BSE 
susceptibility (Haase et al, 2007; Juling et al, 2006; Sander 
et al, 2004). Studies on polymorphisms in buffalo in 
Anatolia (Oztabak et al, 2009), Pakistan (Imran et al, 2012), 
Indonesia and Thailand (Uchida et al, 2014) show signify- 
cant differences in frequency distributions between buffa- 
loes and cattle. Recently, genotyping analysis on Chinese 
buffalo showed that the distribution frequencies of BSE 
susceptibility related to genotypes and alleles, including the 
23-bp deletion allele (D5;) and 12-bp deletion allele (О), 
were significantly lower than those of healthy cattle and 
BSE infected cattle, indicating that the low PrP expressed in 
buffalo may influence BSE susceptibility. Our later experi- 
ments proved that in tissue of the cerebellum, brain stem, 
mesenteric lymph nodes and bronchial lymph nodes, the 
expression of PrP is lower in buffalo than in cattle 
(submitted data). 

Although PRNP play a vital role in the pathogenesis 
of TSE, the underlying pathological mechanisms of TSE 
remain unclear. Some propose that other than prions, 
there may be other factors or proteins regulating the 
pathogenesis and pathological process of TSE (Daude & 
Westaway, 2011; Watts et al, 2007). The SPRN (shadow 
of prion protein) gene and its encoded protein Shadoo 
(Sho) have drawn lots of attention due to their roles in 
the pathogenesis of TSE. Comparative genomics analysis 
indicates that Sho is a newly discovered member of the 
prion protein family. Sho has been found in mammals 
such as mice and humans and is highly conserved from 
fish to mammals (Premz et al, 2003). Sho and PrP“ have 
a lot in common regarding structure and expression 
(Wang et al, 2014). In PrP* infected animal brains or 
nervous cell, with increasing PrP* expression, the level 
of Sho decreases dramatically (Watts et al, 2007, 2011; 
Westaway et al, 2011). Beck et al (2008) reported that 
the insertion of a base (heterozygous) within the 
encoding area of SPRN induces a frame-shift mutation 
which is correlated with vCJD. So, it is highly possible 
that Sho regulates the process of TSE by functioning as 
an inhibitory factor (Daude & Westaway, 2011). Our 
analysis of differences in the genetics and expression of 
SPRN between buffalo and cattle show that in the 
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hydrophobic domain (HD) within the encoding area, 
cattle have a 12-bp indel polymorphism which induces 
insertion/deletion of four amino acids. However, this 
phenomenon was not observed in buffalo (Zhao et al, 
2012). The HD of Sho not only protects against 
physiological stressors in nervous cell, but also helps Sho 
to interconnect with PrP^ (Wang et al, 2010). This 
interconnection is a prerequisite of Sho regulating the 
pathogenesis of disease (Wang et al, 2010). The 
exploration of the indel polymorphism within the Sho 
HD structural area is critical in fully understanding 
underlying mechanisms of TSE. Our luciferase reporter 
and immuno-blotting experiments confirm that compared 
to cattle, buffaloes have higher promoter activity and 
higher Sho expression, consistent with our prediction 
that buffalos have more transcription factor binding sites 
than cattle (Zhao et al, 2012). These findings suggest that 
the low susceptibility of buffalos to BSE is probably 
attributable to significant genetic differences in SPRN. 


Further Research 

Joaquin Castilla’s research group denies there are 
TSE resistant mammals, and believes that with impro- 
vements in detection any species can be found to be at 
risk of TSE infection (Fernandez-Borges et al, 2012). 
However, from available data, TSE susceptibility does 
vary between species and we can classify animals as high 
susceptibility species and low susceptibility species. 
Scientists worldwide have applied various techniques to 
the study of TSE pathogenesis and have mainly focused 
on the genetic polymorphism of PRNP, expression levels 
of PRNP/PrP, three-dimensional structure and stability of 
PrP“ and dynamics of PRNP. However, the pathogenesis 
of TSE remains unclear. Lin & Wen (2011) claim that the 
Ргр and the 
physiological function of PrP^ are keys to resolving this 


three-dimensional structure of 
puzzle but these two research directions have proved 
extremely difficult, even after two decades of attention. 
With breakthroughs in novel technologies and methods 
however, progress is likely. Studies on TSE low 
susceptibility molecules and newly discovered SPRN 
(Wang et al, 2014) provide important clues about the 
formation of PrP% and our understanding of TSE 
pathogenesis. Due to similarities in Sho and PrP© 
regarding structure and function, especially their 
important roles in the pathogenesis and development of 
TSE, exploration of the biological functions of Sho and 


its regulatory effect on TSE will be vital. 
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